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Design of a Hydraulic Dexterous Manipulator for 
Minimally Invasive Surgery 


by Devin Rodney Berg 
ABSTRACT 


The research described here identifies the limitations of existing robotic surgi¬ 
cal platforms, which include the balance between the scale of the robot and its 
manipulability in terms of range of motion, load capacity, and tool capability, 
then develops a means of overcoming them by taking advantage of fluid power as 
an enabling technology with its inherent power density and controllability. The 
approach described here differs significantly from conventional surgical robots in 
that the robot is embedded within the surgical device itself, whereas in the con¬ 
ventional system, a general-purpose robot is used to manipulate various surgical 
tools. This is done in order to demonstrate that fluid power can be used advanta¬ 
geously for the design of embedded surgical robotic systems for minimally invasive 
surgery. 

To enable the design of a fluid powered surgical robot, it was first necessary to 
identify the design requirements for a robot of this nature as well as the consid¬ 
erations unique to this approach. To this end. a quantification of the necessary 
load capacity for natural orifice robots was conducted. Further, through a review 
of the literature in the fields of surgery and robotics, considerations of necessary 
workspace and limitations for the prevention of tissue damage were explored. The 
results of these analyses are presented. 
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The technologies that comprise thus novel surgical robotic system include a hy¬ 
draulic control valve, actuation units, and an enabling structure. The intended 
application of these technologies introduced numerous limitations and challenges 
to the design process. The most stringent of these limitations was that of overall 
size, due to the realities of patient anatomy, which prevented the use of commer¬ 
cially available hydraulic components. An assemblage of components to achieve 
the aforementioned design requirements is described including the design of a 
novel hydraulic control valve to enable manipulation of three actuators using a 
single valve sized to fit within the working channel of a surgical endoscope. 

The advantage of the described approach is that the device enables greater minia¬ 
turization, improves cost effectiveness, and has better ease of mobility. The mobil¬ 
ity and the relaxed requirements for operating room cleanliness can be potentially 
useful for mobile clinics, out-patient clinical settings, and on the battlefield. Being 
more cost effective and having a small overall size, the robotic assisted surgical 
devices can be widely deployed, even in rural or other less technology intensive en¬ 
vironments. Through careful review of the literature and analytical evaluation of 
the various proposed concepts, it was possible to arrive at a design that meets the 
needs of modern surgical interventions while addressing the perceived limitations 
of existing surgical robotics. 

Through the efforts described in this dissertation, much new information was 
produced and developments resulted. The considerations of hydraulic power for 
surgical robots were evaluated and are applicable to other surgical tasks where 
hydraulic power may be used advantageously. A quantification of the load require¬ 
ments for surgical robots performing abdominal procedures was produced which 
will provide a guide for other researchers developing surgical robots. These values 
are difficult to find in the literature and are a valuable resource for the field. An al¬ 
ternative, simplified model for predicting the behavior of continuum beams under 
load was developed to provide an inverse formulation for computing beam shape 
and end loads. This is useful as continuum beams are widely used for minimally 
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invasive surgical manipulators as well as in a wide variety of other applications. 
Finally, a novel valve concept and two possible designs realizing this concept were 
developed. These valve designs fac ilitate control over the three actuators in an 
antagonistic arrangement. Further, the valve designs enable proportional control 
of the three actuators at a size scale not commercially available. In summary, 
the design of a novel hydraulic surgical manipulator as a summation of its parts 
has been performed. This design demonstrates the' feasibility of the fluid power 
approach to embedded minimally invasive surgical robotics. The pursuit of this 
research has provided many unique challenges and the work presented here has 
addressed many of them, as well as laid the foundation for future developments 
in the application of hydraulic power to the growing field of surgical robotics for 
minimally invasive surgery. 
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Chapter 1 


Introduction 


As the field of surgery developed through the centuries from early experimen¬ 
tation to the advanced procedures common today, so too have the design and 
implementation of novel surgical tools. This introductory chapter first briefly lays 
out the changes observed and the changes that ore yet in progress in the surgical 
practice and the way that procedures are performed. Next, some background and 
discussion of the current state regarding the use of robotics in the operating room 
is presented. Finally, the topic of this dissertation is defined t hrough discussion 
of the motivation, n statement of the problem, and a vision for the future. 


1.1 Background and Related work 

Modern surgical intervention can trace its history to the work of Ambrose Pare in 
16th century France JlJ: however, evidence of more primitive procedures and ex¬ 
periments have been found dating hack several t housands of years earlier |2). Since 
that time, the field has undergone extensive transformation as surgical knowledge, 
as well as technology, cont inues to improve. The majority of these advancements 


1 


2 

have ultimately led to improved efficiency, safety, irncl overall effectiveness while 
continuing to progress towards increasingly less invasive surgical techniques. 

In terms of surgical procedures, there are typically three levels of invasiveness 
that describe a particular intervention. These consist of non-invasive, minimally 
invasive, and invasive arranged in order of escalating disruption to normal body 
function. The term non-invasive refers to procedures during which no perforation 
of tissue occurs and are typically of a diagnostic or therapeutic nature. This 
subset of surgical procedures will not be discussed in detail at this time. 


1.1.1 Open Surgery 

Open surgery is commonly considered an invasive procedure due to the use of a 
large incision to gain access to the body. Despite this fact, a significant num¬ 
ber of surgical procedures ;ire still performed using open surgery methods. The 
reasons for this are likely attributed to the difficulties associated with adapting 
complex procedures currently performed using open surgery. Translating these? 
procedures to the methods of more desirable minimally invasive techniques such 
ns laparoscopic surgery, single port access surgery, or natural orifice surgery re¬ 
quires significant effort due to the long learning curve necessary for surgeons to 
develop the required skills in these areas PJ. However, advancements in these ar¬ 
eas are quickly allowing an ever increasing variety of procedures to be performed 
using one or more of these techniques [1]. 


1.1.2 Laparoscopic Surgery 

Laparoscopic surgery is a form of minimally invasive surgery that can be traced 
back to Hans Christian .locobaens who performed the first such procedure on 
humans in the year 1910 [SJ. This technique can be characterized through the 
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use of small incisions in the abdomen of the patient in order !o gain access to the 
peritoneal cavity, as shown in Fig. |1. I| These smaller incisions have the effects of 
reduced post-operative pain, shorter recovery time, and less visible scarring Jfi|. 
However, due to the use of small incisions, it becomes necessary to introduce some 
method of visualization which usually occurs through the use of a charge-coupled 
device (C’(’D) and fiber optics which allow an image to be collected at the distal 
end of the laparoscope and displayed on a monitor for the surgeon. Visualization 
in this manner limits the surgeon's perception and ultimately results in longer 
surgery times and potentially increased risk of surgical error [71 ® ID]- These 
effects can be compounded by t he fact that the surgeon also has limited dexterity 
when using laparoscopic tools |T0]. The typical laparoscopic tool has four degrees 
of freedom (DOF) and requires the surgeon to operate using inverse motions due 
to the fulcrum effect at the location of tool insertion. Laparoscopic surgery would 
greatly benefit from the introduction of tools with greater freedom of motion and 
int nit ive manipulnt ion. 



Figure 1.1: Example of laparoscopic surgery showing access ports and a visual¬ 
ization port PH- 
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Despite these perceived shortcomings, laparoscopic surgery has become a popu¬ 
lar option for several routine surgical interventions. In particular, laparoscopic 
cholecystectomy, or removal of the gallbladder, has become one of the most wide 
spread laparoscopic procedures after being first presented by Dr. Reddick in 1989 
|12| . Further, a panel representing the National Institutes of Health ruled that 
laparoscopic cholecystectomy was a suitable alternative to open cholecystectomy, 
while at the same time warning that the procedure necessitates surgeons of a 
high level of expertise |T3J. This again is likely due to the counterintuitive input 
method. The use of proper surgical robotics would be expected to increase the 
speed at which a surgeon becomes proficient with this technique. 


1.1.3 Single Port Access Surgery 

Single port access (SPA) surgery, also known more recently as lapnro-endoscopic 
single-site surgery (LESS), is similar in form to laparoscopic surgery with the dis- 
t inction of using a single incision to gain access to the abdominal cavity (see Fig. 
1.2[ ), as the name suggests. Recent clinical trials suggest that SPA procedures are 
comparable to their laparoscopic counterparts in terms of operation time and rate 
of complications (Q]. Further. SPA surgery carries the benefits of reduced scar¬ 
ring and possibly shorter hospital stays due to fewer incisions |15[ . One reported 
downside to SPA procedures is a greater difficulty in performing complicated tasks 
due to limited triangulation; however, articulating instruments have proven bene¬ 
ficial and greater tool advancement would likely enhance the range of procedures 
that could he performed using a single port |1-1| . Further, transitioning to on 
endoscopic or serpentine style surgical tool may also provide better performance 
in SPA surgery. 
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Figure 1.2: Example of single port access surgery (Tf>| . 

1.1.4 Natural Orifice Surgery 

Natural orifice surgery is a phrase used to describe surgery in which t he location of 
interest within the patient is accessed through one of the body’s natural orifices, 
such as the urethra, vagina, anus, or esophagus. More specifically, the phrase 
natural orifice translumenal endoscopic surgery (NOTES) is used to describe pro¬ 
cedures during which an endoscope or similar tool crosses the gastrointestinal 
lumen to gain access to the peritoneal cavity as shown m Fig. |1.3| The lirst 
demonstration of a NOTES procedure took place in 1901 bv Dr. Von Ott during 
which the female pelvic cavity was examined by inserting an endoscope through 
the posterior fornix of t he vagina [ITj. Since that time, steady progress has been 
made to bring NOTES to the state of development in which it exists today. Most 
notably, the work of Dr. Knlloo et ah, published in full in 2004, Is often cred¬ 
ited as the first demonstration of the feasibility of oral, transgastric access to the 
peritoneal cavity in a porcine model [TS] |T9] . However, it should also he noted 
that Dr. Reddy and Dr. Rao have described performing NOTES procedures on 
humans in 2003 that went without publication (TSJ. Ultimately, the work of Dr. 
Kalloo et al. led to a joint effort between the American Society of Gastrointestinal 
Endoscopy (ASGE) and the Society of American Gastrointestinal and Endoscopic 
Surgeons (SAGES), which has defined the field of NOTES research and established 
a working group called the Natural Orifice Surgery Consortium for Assessment 
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and Research (NOSCAR) in order to provide structure to the rapidly developing 
field of NOTES [RJ- 


a 



Figure 1.3: Example of transgastric natural orifice surgery [30] . 


The primary advantage of NOTES procedures is the reduced impact on patients, 
which includes no visible scarring, less post-operative pain, and decreased time to 
recovery [Til . Further, analysis of a collection of 180 NOTES procedures found 
that the rate of complications was comparable with laparoscopy while there was a 
significant reduction in the use of medication for pain, post-operativelv [22[ . The 
importance behind these outcomes is significant as they have greatly contributed 
to the transition from open surgery, to laparoscopic surgery, to SPA surgery. 
Presently, the majority of NOTES clinical trials have been performed using tra¬ 
ditional endoscopes adapted for the particular task; however, new advancements 
in specialized devices are becoming more common and will be discussed further 
in Section ll.l.~>l 








1.1.5 Robotics in Surgery 
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The development and application of robot ics for the operating room is a field that 
has undergone significant transformat ions in recent years as the push towards min¬ 
imally invasive surgical techniques requires an increasing level of surgical preci¬ 
sion and repeatability. Advancements in minimally invasive surgery have enabled 
robots to replace the camera-person, provide a more ergonomic environment for 
the surgeon, and even provide the means for a surgeon to perform a procedure 
from thousands of miles away |23J. Further, surgical robots can be employed as 
a tool for minimizing surgical errors such as those that may be caused by hand 
tremors and have been shown to improve the efficiency mid speed of even novice 
surgeons [24| . 

The earliest recorded use of robotics in surgery can be found in the work of Shoo 
et aL in 1985 when a robot designed for industrial applications was used to per¬ 
form a neurosurgical biopsy [251This work formed the foundation for a field 
that produced the first FDA approved surgical robot, the ROBODOC' (Integrated 
Surgical Systems Inc.) for use in hip replacement surgery [2711281 . More recent 
technologies in surgical robotics include the AESOP and Zeus robotic systems 
(Computer Motion) !is well as the da Vinci Surgical System (Intuitive Surgical). 
After the merger of Computer Motion and Intuitive Surgical, the Zeus System 
was discontinued in favor of the da Vinci Surgical System. The da Vinci Surgical 
System consists of two primary pieces of equipment shown in Fig. [l.-l[ The first 
is a surgeon console which provides feedback from the robot to the surgeon and 
accepts input from the surgeon to control the robot. The second Is the patient 
cart which includes a bed for the patient and the laparoscopic surgical instruments 
that, perform the procedure. The tools of the da Vinci system are manipulated 
by an extensive arrangement of tensile wires. These wires run the length of the 
robotic arm and move in controlled sequences in order to articulate the tooLs in 
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the desired manner. As can be seen, the system is rather large due to the technolo¬ 
gies it employs and is prohibitively expensive for many surgical facilities. Despite 
this, the da Vinci Surgical System has become the standard in surgical robotics. 
However, ns surgical techniques have advanced towards more difficult procedures 
and towards less invasive access methods, greater sophistication and ability Is re¬ 
quired of surgical robotics. Further, miniaturization of the robot would allow the 
robotic surgical system to be used in smaller surgical suites and possibly mobile 
environments. 



Figure 1.4: da Vinci Surgical System including surgeon console (left) and patient 
cart (right) [20] , 

The field of laparoscopic surgery 1ms been particularly receptive to the use of sur¬ 
gical robotics. This has allowed the da Vinci Surgical System to become the only 
widely distributed, commercially available surgical robot currently on the mar¬ 
ket. Presently, there are more than 2000 da Vinci systems installed at hospitals 
around the world and thus the conversion of traditional laparoscopic procedures 
to robotic laparoscopy is occurring rapidly [IQ]- For example, the prostatectomy 
performed on the da Vinci system Is becoming the standard treatment for prostate 
cancer £TT]. This level of market penetration Is helping to bring patient accep¬ 
tance to the held of surgical robotics. However, as advanced surgical procedures 
transition towards a laparoscopic approach or other less invasive access methods. 
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improvements are necessary both in terms of form factor and tool capability. 


Modular Robots 

The next realm of development for surgical robotics appears to be headed towards 
the even less invasive surgical access techniques such as single-port access or nat¬ 
ural orifice access [3=]. There me several emerging technologies related to surgical 
robotics that fall into broad categories ranging from modular robots to flexible en¬ 
doscopic or serpentine robots (TTJ. Beginning with modular robots (meaning that 
they consist of small standardized units that can be assembled or arranged), the 
work of Slat kin et ah. shown in Fig. |1.5[ describes the development of a robot ic 
endoscope intended to gain access to the small intestine without the use of an ex¬ 
ternal incision, something that was previously not achievable through the use of 
traditional endoscopes |T1]. The device described therein employs pneumatic ac¬ 
tuators to achieve various states of expansion and contraction to produce motion, 
similar to that of a segmented earthworm. Further, alternate segments of the de¬ 
vice provide radial expansion which permits the device to achieve traction against 
the intestinal wall. Thus, by controlling both the radial and longitudinal expan¬ 
sion and contraction of the various segments, locomotion of the device tltrough 
the intestines is possible. Experimental studies showed that the robot performed 
well even in situations involving flexible and slippery terrain. Although this device 
creates a basis for a whole new held of possibilit ies through surgical robotics, it 
provides imaging only with no tool capabilities. 

Further developments in the area of modular robots axe beginning to address 
the issues with positioning, manipulation, and teleoperation |33j . For exmnple, 
significant advances have been made in these areas using a remotely controlled 
single-port surgical robot known as SPRINT [351 |3fi| . The SPRINT robot is 
a master/slave system which consists of two independently manipulated arms 
intended to be introduced to tire abdominal cavity through a single 3 cm access 
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port located in the patient's navel. The joints of each arm are driven by DC 
motors (both external and embedded) and gears for each degree of freedom at the 
"shoulder", "elbow ", and "wrist" joints. This device demonstrates the capabilities 
of these forms of robots in terms of dexterity and maneuverability fay realizing 6 
DOF per manipulator and a load capacity of 5 N [3D]. While the geometric scale of 
this design is larger than the application presented here (maximum outer diameter 
of 18 mm), it does show feasibilit y of using embedded robotic components. Lessons 
learned from the development of robots for SPA procedures can be used to inform 
the design of natural orilice robots. 



Figure 1.5: Endoscopic robot design described by Slatkiu et al. [34] . 

While partially addressed by the work of Slatkiu et al. |33], one constraint on 
novel surgical robots that attracts significant attention is the level to which the 
robot is connected to the outside world while deployed in the abdominal cavity. 
One method for addressing this issue is to use miniature in vivo robots [331I37[. 
The complexity of these devices ranges from stationary or mobile devices with 
dedicated tasks such as visualization or lighting to robotic assistants that are ca¬ 
pable of providing both visualization and tissue manipulation [Rj[33]. Particularly 
interesting are a recent category of devices which make use of magnetic anchor¬ 
ing to physically isolate the device from the external world while maintaining the 
control over positioning and manipulation [3311301 t-TOl |TI| . One passible realiza¬ 
tion of this concept is the use of multiple specialized robots deployed through 
the same access channel that work together to collectively provide visualization, 
lighting, grasping, cutting, and cautery [331 [391 32]. Thus through the strategic 
placement of the various devices, procedure's such as a cholecystectomy cau be 
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performed with no external incisions [3BJ. However, one disadvantage to this ap¬ 
proach is that it introduces an extra layer of complexity to the surgical procedure 
since each of the various components must he managed independently. Another 
possible approach, is to use a more advanced in vivo robot that incorporates vi¬ 
sualization iuid the necessary tools into one unit such as is described by Lehmsui 
et al. |-ll| . This work presents an in vivo robot consisting of a main body housing 
lighting and visualization, as well as two arms located on either end of the body 
that provide grasping and cautery dissection. The device demonstrates promise 
for the future of advanced in vivo robots; however, issues with robustness of the 
device and insufficient tool maneuverability need to be addressed for such a device 
to be successful QTj. In particular, the use of rigid arm-like structures severely 
constrains the tools ability to manipulate safely within the abdominal cavity due 
to t he necessity for large sweeping motions of the arms resulting in unintended 
collisions with abdominal tissue |J3J. Further, the tool has the disadvantages of 
requiring assembly once inserted into the abdomen and using motion mechanisms 
that: make its movements too slow for practical application. 


Serpentine-like Robots: Jointed 

Another approach Ls the use of endoscopic or serpentine-like robots which have 
advantages over modular robots in terms of large flexibility for minimal size [3ti| . 
Many of the clinical trials performed thus far in this direction have employed tra¬ 
ditional endoscopes by passing the necessary tools through the scope's working 
channels. However, standard endoscopes are ill-equipped to provide the necessary 
freedom of tool manipulation and force capacity for many advanced Abdominal 
procedures [311 UTl jJoj. Therefore, novel endoscopic platforms are rapidly being 
developed to address the concerns of these advanced procedures. Robots of this 
variety take their inspiration from nature utilizing the concepts of biomimicry and 
bioinspired design Plti, |T7| . Many of these devices utilize a series of rigid joints 
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driven by robotically manipulated cables to achieve tool positioning and manip¬ 
ulation such as is seen in the work of Phee et al. and others J3S 1391 1501 [5T|. 
Such devices are subject to several constraints that must be considered for them 
to be successful. These include the diameter of the instrument which should be 
approximately 20 mm if adapted to NOTES procedures, a minimum of two to 
three working channels, and flexibility in steering and control of position fXfl [52] . 
These devices operate in a similar manner to the da Vinci robot previously de¬ 
scribed. while many of the same limitations remain such ns tool manipulation and 
triangulation. This is particularly true in the application to transgastrle NOTES, 
which often requires greater instrument flexibility than rigid laparoscopic tools 
can provide |53| . Phee et al. describe the design of jointed manipulators intended 
for use in endoscopicallv accessed abdominal procedures [38| . The robotic manip¬ 
ulators me controlled by externally located DO motors and long tendons. Each 
manipulator has 5 positioning DOF realized through the use of serially arranged 
joints and an overall diameter of approximately 10 mm. Subsequent test ing of 
this device revealed the ability to produce end loads in the range of 2.!) N to 5.2 N 
[491 . This device demonstrates significant capability for an endoscopicallv intro¬ 
duced manipulator but requires the use of externally located power components 
to achieve these results. The work of Kobaynshi et al. [50] presents a surgical 
tool for SPA procedures which consists of an external positioning manipulator 
and internal tool manipulators. The overall tool has a diameter of .'1 cm and the 
individual tools have diameters of (> mm and 8 mm. The positioning of the tools 
axe externally controlled using rotating shafts and screws rather than long ten¬ 
dons thus avoiding the issues associated with long tendon manipulation such as 
sticking, slipping, and jerking [ofl| . Each tool has a bend range of -15 c which may 
be seen as a limitation in terms of the desired surgical workspace requirements. 
Finally. Shang et al. demonstrate the capability of a tendon driven robot to actu¬ 
ate through a large range of motion (±90° vertically and ±225° horizontally) by 
routing the tendons around universal joints to magnify the output motion (oil 153] . 
Tills design makes use of DC' micromotors with a diameter of -I mm to drive the 
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universal joints. This has the benefit of enabling embedded actuation, but comes 
at the cost of load capacity being limited to 0.5 N at the tool tip and an overall 
outer diameter of 12.5 mm. This design shows a lot of promise but currently has 
severe limitations on output load capability. Improvements in load capacity and 
greater miniaturization would make this device an attractive option for a surgical 
robot with embedded components. 


Serpentine-like Robots: Continuum 

Similar to the jointed, hvper-redundant robots presented in the previous section 
are a class of robots which are based on the use of a continuous flexible structure. 
Continuum manipulators have the advantage of improved simplicity relative to 
their jointed counterparts thus enabling greater miniaturization which is critical 
for many surgical applications as well as elimination of backlash that is present in 
jointed structures. On the small scale, the design of continuum manipulators is 
typically coupled with the use of tendons due to the geometric limitations along 
the length of the manipulator. One interesting example of a flexible serpentine 
robot design presented by Webster et al. is based on pre-tuned, telescoping tubes 
fo5| . The concept consists of tubes made of superelastic materials, such as nitinoi, 
that are pre-curved in such a way that ns they are extended from a straight semi¬ 
rigid outer tube the elastic energy stored within the material produces a known 
curvature. Subsequently, through the layering of multiple pre-curved tubes, the 
device can achieve complex shapes through variations in longitudinal and rota¬ 
tional positioning. However, since the device relies upon the strain energy stored 
within the pre-cun*ed tubes in order to achieve manipulation passively, it is sub¬ 
ject to severe limitations in situations where the application of a force is necessary. 
This introduces difficulties when it comes to procedures that require significant 
tissue manipulation. More recently this work has been adapted for use with pneu¬ 
matic actuation in combination with the pre-curved concentric tube manipulators 
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for surgical treatment of epilepsy [561 I57[ . This development represents an inter¬ 
esting direction that may lead to greater capability for the pre-curved tube based 
manipulators and presents an interesting example of a fluid power approach. How¬ 
ever. the current iteration includes large external frame to house the manipulation 
components. Further, the manipulator has a relatively short range and would re¬ 
quire the access point to be located in close proximity to the surgical site thus 
making it difficult to apply this technology to endoscopic access techniques in its 
current form. 

As previously mentioned, many of the upcoming designs in surgical robotics fea¬ 
ture cable or tendon style inputs to drive the end effector. The ViaCath System 
(EndoVia Medical) presented by Abbot et id. [T5] consists of two (i DOF ma¬ 
nipulators which are externally manipulated bv servomotors and introduced to 
the abdominal cavity endoscopically. Each manipulator Is approximately 5 mm 
in diameter. In experimental analysis the tool was found to be capable in terms 
of positioning and range of motion but was limited to 0.5 N of load at tlie distal 
end which was determined to be insufficient for most procedures. A secondary 
tool design was also presented which consisted of a jointed structure of 7.2 mm 
diameter, similar to the robots presented in the previous section. The second 
design demonstrated greater load capacity but was subject to the need for ad¬ 
ditional controlled inputs to manipulate the given joints. More recent work in 
this area has demonstrated that miniaturizat ion of the continuous backbone style 
manipulator is passible and lars been demonstrated down to a diameter of 4.2 mm 
[561159) . The distal manipulators of the presented design each have 4 DOF and 
are manipulated by an externally located actuation unit. Testing of the proto¬ 
type revealed uncertainties due to the use of long tendons which extended from 
the distal manipulators back to the external actuation units. Each manipulator 
was shown to be capable of bending 90° for the proximal segment and 60° for the 
distal segment meaning that the manipulator as a whole could achieve a bending 
range of 150°. However, the tool load capacity is only 1 N at the tip which may 
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not meet the requirements of common surgical procedures (see Section 27T|. Ding 
et al. [HD] and Xu et al. |ilT] both present extensions of this design which use 
the same tendon style control to operate two manipulators and a vision unit on 
the same system which can be deployed endoscopicallv. This design features the 
ability to fold each of the three parts together for insertion and then deploy them 
for the procedure. The manipulators are then controlled via external actuation 
units similar to the previous design. Another feature of this design is an innovative 
method for suturing tissues through the use of pre-curved nitiuol needles based on 
the work of Webster et al. |o5| . This approach may alleviate part of the need for 
higher tool force capability. However, it severely limits the adaptability of the tool 
during a procedure by having a suture needle of specific size that cannot make 
use of existing suture needles already available in the operating room. Further, 
each of these designs requires a large external actuation unit to control the device 
which limits portability and introduces the potential for lag and inconsistency 
when controlling the unit. 


In addition to the issues .surrounding miniaturization, another major challenge 
faced by designers of continuum robots is that of producing accurate models of the 
robot's performance or behavior. These models are necessary both to facilitate a 
better understanding of the device as well as to enable design of a controller for the 
device. Early models were developed under the assumptions of hyper-redundant. 
rigid link based devices p21[?>!?I- While these models are quite effective, they do not 
account for the specific behavior of continuum devices present, due to the existence 
of an infinite number of degrees of freedom |f>3|. Much work has been performed 
to model the behavior of the continuum manipulator using a variety of methods 
[tl?>l Ififil [fi71 IfiS]. While these models provide great insight into the behavior of 
the continuum manipulator, they often require significant computational effort to 
arrive at a solution using numerical simulation. However, in the case of a low 
inertia system it is often possible to simplify the given models in such a way as to 
allow for more real-t ime calculation or simulation efforts for prediction of dynamic 
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behavior. The desire to develop models with greater accuracy and less complexity 
is an ongoing effort and the best solution to this problem continues to evolve 

mm- 

Ultimately, robotic surgery presents several possible benefits over manual sur¬ 
gical techniques. As previously mentioned, robotic control of a surgical device 
introduces the possibility of reducing surgical error through cancellation of hand 
tremors or other methods such as image recognition (23|. Further, robotic sur¬ 
gical instruments have been shown through clinical studies to improve surgical 
efficiency, which translates to greater patient throughput and cost savings when 
compared with traditional procedures [32| . Beyond these benefits, other possibil¬ 
ities for an improved surgical experience include force feedback and telesurgery. 
Force feedback is one component of a surgical procedure that is often missing in 
robotic surgery. However, it has been shown that force feedback provides impor¬ 
tant information to the surgeon that can help to prevent tissue damage or other 
unintentional errors [7fil fTTl I72[ , Therefore, the inclusion of force feedback in any 
robot ic surgical instrument should be given strong consideration. Another benefit 
of robotic surgery is the potential for teleoperation. As discussed in the literature, 
teleoperation permits surgeons to operate from a more ergonomic position, which 
reduces fatigue and allows the surgeon to be located at a great distance from 
the operating room [23|. Tills in turn opens up a number of possibilities, such 
as surgical practice in under-developed countries where surgeons may be in high 
demand. Further, by operat ing the tools robotically, the task can be performed in 
a way that is more intuitive for the surgeon, with more natural kinematics. With 
all of these potential benefits, development in the held of surgical robotics is likely 
to continue at a rapid pace. Through the development of a surgical robot for MIS. 
it is possible to enable these types of advancements to better equip surgeons as 
they perform complicated tasks. The advent of endoscopic robots for perform¬ 
ing surgical procedures enables such developments as natural orifice access tech¬ 
niques through inherent flexibility and maneuverability. However, even though 
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eudoscopic instruments provide greater flexibility, they are still hindered by lim¬ 
ited stability while performing delicate procedures. Swnnstrom ct aL present a 
possible solution to the Issue of scope rigidity |73| . The device uses ShapeLock 
technology which allows the scope to be oriented in the correct surgical position 
and then locked into place through the application of a compressive force within 
the shaft. This platform provides a useful base from which flexible surgical in¬ 
struments may conduct the procedure and presents the possibility of greater tool 
capability. 

The designs discussion in this section have primarily achieved their design re¬ 
quirements either through the use of external actuation units and long tendons 
which introduce uncertainty and a huger operating room space requirement or by 
increasing the size of the manipulators to accommodate embedded components. 
While each of the designs presented here show promise, there remains the need 
for greater capability at a smaller scale to enable the design of robotic surgical 
manipulators which can be embedded within the working channels of existing 
endoscopic instruments. 


1.2 Motivation 

The use of robotics in surgery is becoming increasingly necessary with the advent 
of new minimally invasive surgical techniques such as single port access (SPA) 
surgery, a single incision procedure dining which access to the abdomimd cavity 
is typically gained through the navel, and natural orifice translumenal endoscopic 
surgery (NOTES), an experimental technique that involves performing traditional 
surgic al procedures with a flexible endoscope through a natural orifice such as the 
mouth, rectum, or vagina. These sorts of novel techniques remove the need to 
create incisions ncross the abdominal wall, thus minimizing or even eliminat ing 
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past-operative pain, abdominal scars, and complications such as incisional her¬ 
nias. Other potential benefits include reducing postoperative recovery time, using 
moderate or deep sedation instead of general anesthesia, and the potential to per¬ 
form many surgical procedures outside of traditional operating rooms. Because 
of these advantages, techniques such as SPA surgery and NOTES appear to he 
the next natural step in the progression of minimally invasive surgery for routine 
procedures such as cholecystectomies and exploratory laparoscopies in hospitals 
around the world [RJ. Due to the complexity of these procedures and the nat¬ 
ural limitations of space and site access, the development of advanced robotics 
becomes imperative for the future success of these sort of surgical interventions. 

One approach to the development of surgical robots is fluid power, which can come 
in the form of pneumatic or hydraulic. In particular, hydraulic power advanta¬ 
geously has higher power density (on a per weight or per volume basis) on the 

order of 10:1 or greater when compared with other driving methods [73M75J- Ad¬ 
ditionally, a hydraulic operating fluid such as saline would provide less uncertainty 
due to compressibility than a pneumatic fluid and thus enable greater precision. 
However, the underlying challenges surrounding this .assumption include a ques¬ 
tion of what limitations exist on the miniature scale for which the benefits of a 
hydraulic approach hold true. For example, work is currently in progress explor¬ 
ing the use of fluid power at the meso-scale (mm to cm) for use in prosthetics [7(i| . 
The next logical step is to explore the use of fluid power at the micro-scale (/nil to 
nun) for which the research is ongoing in the field known as microfluidics [771 ITS] , 
However, as pointed out by Love et al.. there exists a gap between conventional 
fluidics and inicrofluidics that includes the type of high pressure-low How rate fluid 
power components that would be necessary for the application of hydraulic power 
to the held of surgical robotics [Tfij. Therefore, it is necessary to conceptualize 
the methods by which these sorts of fluid power components might be realized. 

For surgical robots to adapt to these developments, it is necessary to overcome 
significant challenges in terms of tool manipulability and precision. By embedding 
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the robotic components close to the manipulators, it is possible to reduce the 
uncertainty present within the system while achieving similar tool capability. Most 
systems rely on a large external robot manipulating cable driven tools which 
inherently result in uncertainty due to the use of long tendons piSl 179]. Due to 
advantages in power density and relative incompressibility [ 75] , hydraulics may 
allow surgical robotics to be embedded endoscopically. Yet. the use of hydraulics in 
this field has been under-explored due to difficulties associated with the integration 
of hydraulic components into the distal end of the robotic manipulator due to size 
constraints JSO]. The few attempts at using hydraulics for surgical robotics have 
typically been limit eel to external manipulation of the surgical tools rather than 
an embedded design [8]] IS2J - 


1.3 Problem Description 

There exists a wide range of abdominal procedures for which a robotic surgical 
tool may be employed. Access to the abdominal cavity may lie gained tlirough 
abdominal incisions, the patient's navel, or natural orifices such as the month, 
vagina, or rectum. In the case of oral access, the surgical tool deployed via an 
endoscope which in turn would be passed through an incision in the stomach wall 
and into the abdominal cavity. Prom this trnnsgastric access point, the surgeon 
could perforin the necessary tasks ranging from gallstone removal to gastric by¬ 
pass surgery. The robot as a surgical tool provides the surgeon with the ability to 
perforin internal procedures without direct line of sight or direct external access 
to the surgical site. The greatest difficulties faced when applying robotics to these 
types of procedures come in the form of severe constraints on the design in terms 
of force requirements, range of motion requirements, and limited available space. 
When examining the variety of surgical robots in the literature, it Is apparent that 
there is a trade-off between size ;uid functionality. The goal of this work was to 
demonstrate that hydraulic power can be employed advantageously in the design 
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ol a robotic surgical tool for miniinnlly invasive procedures and enable the use of 
embedded robotic components to achieve a smaller overall footprint within the 
surgical suite. Size constraints for endoscopic surgical devices are dictated by the 
geometric constraints of the patient’s anatomy. The outside shaft diameter for the 
device is typically limited to a range of 18 to 22 millimeters based the average size 
of the patient’s esophagus |"2| . Contained within this diameter, it is necessary to 
include a channel for imaging and at least two tool channels, as well as additional 
ports for insuHiation, suction, and irrigation. Considering these requirements, a 
practical limit on the tool channel diameter of 5 to 6 mm is typical, it is wit hin 
this available space that the mechanisms for both device articulation and tool 
manipulation must be contained. Further, the surgical instrument needs to be 
capable of 180 ° retroflexion in order to navigate the complex abdominal anatomy 
|j-| . This requirement is particularly important for oral access techniques since 
it is not possible to reach the upper quadrants of the abdomen without retroflex¬ 
ion. Looking at the force capabilities required from the tool, it is important to 
consider the forces necessary for both gross tissue manipulation and common pro¬ 
cedural applications such as suturing. Through experimental investigation it was 
determined that a load capability of approximately 5 X would provide sufficient 
tool capability. The lack of a clear solution for a robotic device that meets these 
requirements provided the impetus to pursue a novel design. 

To achieve these goals, it Is necessary to consider each of the various subcom¬ 
ponents that comprise the assembled surgical tool and the interactions between 
them. These components provide the mechanics of tool positioning through the 
support structure, actuation unit, and method of control input. First, the struc¬ 
ture of the tool which provides support and shape needs to be determined and 
modeled. Next, it is necessary to determine the method by which the tool is ma¬ 
nipulated through the use of an actuator or motor. It Is then necessary to specify 
the control method for the device. Finally, the interactions between the individual 
subcomponents and the system as a whole needs to be analyzed and understood 



in older to evaluate how the tool internets with its environment. 
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1.4 Vision 


The work described here intends to identify the limitations of existing robotic 
surgical platforms and develop a means of overcoming them while satisfying the 
tool requirements of force output: and mauipuiability. Currently available surgical 
robots are developing rapidly, but a technological gap exists in the area of small 
scale, high capability platforms. It is hypothesized here that this gap can he filled 
by taking advantage of the power density and other properties of fluid power, 
such as incompressibility and efficiency; thus permitting surgical robots to meet 
the needs of advanced minimally invasive procedures. The approach discussed 
here differs significantly from conventional surgical robots in that the robot, is em¬ 


bedded wit Inn the surgical device itself (see Fig. 1.6ft: whereas in the conventional 


system, a general-purpose robot is used to manipulate various surgical tools. The 
advantages of this approach are that with proper technologies, the robotic ma¬ 
nipulator can exhibit greater freedom of motion, have better controllability, be 
more cost effective, and have superior ease of mobility. The improved mobility 
combined with relaxed requirements for operating room cleanliness can be po¬ 
tentially useful for mobile clinics, out-patient clinical settings, and on the battle 
field. Being more cost effective and having a small overall size, the envisioned 
robotic assisted surgical devices can be widely deployed, even in rural or other 
less technology intensive environments. Hydraulic actuation typically displays 
greater stiffness and faster response time than the comparable electromechanical 
or pneumatic counterparts. Further, due to their ability to operate at low veloc¬ 
ities with high forces and the lack of an upper limit on power input, hydraulic 
systems can achieve greater power density which allows for miniaturization of the 
device components [STIJ. When compared with a pneumatic approach, hydraulic 
power stores less energy in the form of compressed fluid. Thus hydraulic power 
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can be safely operated at a higher pressure than a similar pneumatic system. Fur¬ 
ther. if a compat ible hydraulic fluid is selected, such as saline, leakage is not likely 
to result in contamination of the abdominal cavity. Lastly, the low costs of the 
components used in this design would permit the device to be disposable and thus 
loosen requirements for sterilization. 



Figure 1.6: Envisioned design of the robotic surgical tool for natural orifice surgery. 
Shown with example end effector tools including a grasper and a scissors [84| . 

Through the development of a novel surgical manipulator. I seek to show that the 
use of fluid power in this application is an achievable goal and would allow for 
greater miniaturization while maintaining or improving robotic capability. The 
surgical robotic platform described here consists of three basic units: manipula¬ 
tors. actuators, and control valves. All components are designed to be located in 
the distal end of the endoscope and contained within its workiug channel with a 
diameter of approximately 6 mm. Figure |1.7| shows a simplified version with two 
manipulators. These manipulators would serve ns a platform for the end effec¬ 
tor tools (see example tooLs in Fig. |1.6[ ). which could include graspers, scalpels, 
scissors, or clip applicators, each manipulated by its own actuator [85j. These 
end effectors could be deployed individually with the manipulator or stored dis- 
tally within a tool holder. Each manipulator can be divided into multiple serially 
connected segments with a single flexible backbone and Is actuated by a set of 
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tendons/'wires, similar to other serpentine-like robots [SB] [87] [88|. The bending 
of each manipulator segment is cont rolled by a group of three tendons antag¬ 
onistically operated by hydraulic artificial muscle actuators at the base of the 
platform. The tendons, which are arranged circumferentially, pass through inter¬ 
mediate discs and terminate at the end of the given segment. The tension in each 
tendon controls the bending of the segment in one direction. Unlike cable driven 
surgical tools that are actuated via long cables connected to external motors or 
actuators, in this design the tendons are kept relatively short since the muscle 
actuators are small enough to be stacked serially and placed (listally within the 
lumen of the endoscope. Hydraulic flow/pressure to the muscle actuators are 
modulated bv valves which are also contained within the endoscope behind the 
muscle actuators. Thus the only connections required externally are lines for the 
hydraulic supply (assumed to be n constant pressure supply) and return as well as 
the electronic signal wires. Figure [P| shows the basic structure of both a single- 
segment manipulator (Fig. |1.8a[ i and a multi-segment manipulator (Fig. |1.8h( . 
Each version is shown with only a minimum of included components and is not 
meant to illustrate the full system assembly. 



Figure 1.7: Model of the robotic surgical tool design shown with two manipulators, 
each with two individual segments. 


The usage scenario for this device would be similar to the scenario used by many 
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Figure 1.8: Line drawings showing a minimal version of how the valve, actuator, 
and tendon components are connected for (a) a single-segment manipulator and 
(h) a multi-segment manipulator. 


similar robotic systems discussed in Section [~1.1.5| In practice, this system would 
consist of an easily deployable bedside cart with a hydraulic supply. Surgical 
input would be provided either on-location or remotely from a surgeon input 
console. The device components are inexpensive enough to be disposable but are 
sealed such that sterilization is possible when necessary. Once in the surgical 
suite, preparation of the patient would involve gaining access lo the abdominal 
cavity using the methods currently being developed in the field of NOTES (Section 
1.1.4( |. Once access is achieved, the endoscope would be positioned near the area 
or organ of interest such as the gallbladder or small intestine. The endoscope 
would then be rigidly locked in place using a technology such as ShapeLock |7H| . 
From this platform the robotic instruments presented here would be deployed 
(see Fig. |1.9| . The surgeon as the operator would manipulate the robot from a 
console similar to what is employed by the da Vinci Surgical System (see Section 
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1.1.5[ . This console would allow the surgeon to control multiple segments per 
manipulator in a manner that mimics natural body motion by measuring mult iple 
joint orientations and mapping them each to the actuators for each manipulator 
segment, for a total of three DOF per manipulator segment. A stereoscopic video 
feed would he streamed to the console from the distal end of the robot using 
liber optics. The manipulators allow the surgeon to remotely perform nil of the 
necessary common surgical tasks such as suturing, tissue manipulation, or tissue 
resection. Upon completion of the procedure, the manipulator are retracted from 
the endoscope to allow for closure of the access site. The relatively low cost of 
the manipulator components would allow the manipulators to he discarded after 
use. In the subsequent pages 1 intend to demonstrate that this vision of greater 
miniaturization while maintaining or improving robotic capability can he achieved 
through the application of fluid power as an enabling technology. 
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Figure 1.9: Illustration of the robotic design during a procedure showing access 
to the abdominal cavity gained through an incision in the stomach wall. 

1.5 Dissertation Overview and Contributions 

This introductory chapter presented a brief history of the evolution of surgical 
intervention as well as an overview of recent developments surrounding the ap¬ 
plication of robotics to the field. While the range of approaches and technologies 
being adapted for use in minimally invasive surgery is vast, it is clem that there 
exists a need for a small scale, highly capable surgical manipulator. A vision for 
the design of a surgical manipulator was presented that takes advantage of the 
properties of fluid power to meet the needs of emerging surgical procedures. 

Through the efforts presented in this dissertation, much new information was pro¬ 
duced and developments resulted. The use of hydraulic power for surgical robots 
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lias not been extensively explored and thus the considerations unique to this ap¬ 
proach which have been evaluated are valuable and applicable to other surgical 
tasks where hydraulic power may be used advantageously. A quantification of the 
load requirements for surgical robots performing abdominal procedures was pro¬ 
duced which will provide a guide for other researchers developing surgical robots. 
These values have typically been reported anecdotally in the literature and thus 
are a valuable resource for the field. An alternative model for predicting the 
behavior of continuum beams under load was developed to provide an inverse for¬ 
mulation for comput ing beam shape and end loads. This is useful as continuum 
beams are widely being explored for use as minimally invasive surgical manipula¬ 
tors as well as in a wide variety of other applicat ions. Finally, a novel valve concept 
and two possible designs realizing t bis concept were developed. These valve de¬ 
signs facilitate control over the three actuators in an antagonistic arrangement. 
Further, the valve d(‘signs enable proportional control of the three actuators at 
a size scale not commercially available. The design of a novel hydraulic surgical 
manipulator as a summation of its parts hits been described. This device builds 
upon the variety of approaches found in the literature while adding several useful 
contributions such its enhanced miniaturization and internally contained compo¬ 
nents. 

In Chapter [2J the design requirements for a surgical robot designed for use in nat¬ 
ural orifice transgastrie endoscopic surgeries are described. These requirements 
include considerations for the size of the surgical workspace, safety and biocom- 
patibility. and load capacity of the tool. Due to a lack of tool force requirement 
quantifications in the literature, a set of experiments were carried out for expected 
high force usage scenarios to support robotic manipulator design through an eval¬ 
uation of tool force requirements. The results presented allow for performance 
analysis of existing surgical tools as well as inform future directions for similar 
devices. 


Chapter [3] provides a detailed description and analysis of the selected backbone 
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structure. An alternative model using the principle of minimum potential en¬ 
ergy and inverse kinematics is presented to predict the behavior of the backbone. 
The model is then experimentally validated to demonstrate that it is capable of 
producing accurate beam shape predictions under most circumstances. Recom¬ 
mendations are provided as to how to best utilize the model under the prescribed 
assumptions. 

In Chapter [1] the selection of an actuator to drive the manipulator is presented. 
The performance of the actuator is predicted analytically and confirmed experi¬ 
mentally. Further, an optimization algorithm is used to determine the ideal pa¬ 
rameters necessary for fabrication of an actuator which maximizes range of motion 
within the specified constraints placer! on actuator diameter and force capability. 
This was accomplished through modeling of the stresses present within the actua¬ 
tor using finite element analysis. It is demonstrated that the requirements placed 
upon the actuator can be met using the given design. 

Chapter [^describes a method of providing control input to the system by way of 
a control valve. To meet the needs of the presented device within the constraints 
of the application, a novel valve design is conceived. Design parameters that are 
necessary to achieve desired valve performance were calculated. Two realizations 
of the design were fabricated at a larger scale to provide a means of evaluating 
the design performance and confirming the hydraulic principles used to develop 
the design. The results presented suggest that the design can be scaled to a size 
appropriate for application to t lie field of minimally invasive surgery. 

In Chapter [(TJ a description of the assembly of the novel manipulator is provided. 
A model is shown which combines each of the manipulator components to allow 
for prediction of manipulator performance. Design specifications are given for 
each of the components that make up the assembly such that the design is theo¬ 
retically capable of meeting the requirements set: forth in Chapter[2] Additionally, 
results from testing of a manually operated manipulator assembly are shown to 


demonstrate efficacy of the design. It is also demonstrated that it was possible to 
achieve the desired performance objectives using this design. 



Chapter 2 


Design Requirements 


The first task in any design problem Ls to lay out both tin* constraints and re¬ 
quirements that help direct the design process towards a workable solution. This 
chapter attempts to do that by identifying the important considerations rele¬ 
vant to the design of a robotic manipulator for minimally invasive surgery and 
then quantifying those requirements using both experimental measurement mid 
information from the literature as appropriate. Consideration is given to the load 
capacity required of the tool, the size of the workspace, and Issues involving safety 
and biocompatibility. Within each of these areas, recommendations are provided 
to help guide the design of the surgical manipulator. 


First, in Section |2.I| the force requirements for a surgical manipulator for use in 
abdominal procedures are experimentally quantified. This is followed by a dis¬ 
cussion of the tool requirements in terms of the workspace, safety considerations, 
and biological constraints in Sections |2.2| |2.3[ anil |2.‘1| respectively. Finally, the 
requirements discussed ill these sections are summarized in Section [To 


:$<) 







2.1 Tool Force Requirements 
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2.1.1 Motivation 

Surgical robots have allowed surgeons lo perform n variety of procedures using 
loss invasive techniques. These new techniques make it possible to minimize or 
even eliminate post-operative pain, scarring, or other complications. Further, the 
advent: of robotic surgical tools allow for the inclusion of technologies such as 
force-feedback to better inform the surgeon and enable greater performance [SflJ . 

It has been shown that expert surgeons exert less force on the tissue when com¬ 
pared with novice surgeons. However, it is important to understand the forces 
necessary to perform a particular surgical task for both t he proper deslgu of sur¬ 
gical tools and to minimize the potential for tissue damage plDllfHl |!)2| . Through 
experimental evaluation of these surgical tasks, it is possible to develop a better 
understanding of the force limitations and better deline the tool force require¬ 
ments. In this section, the force requirements for a surgical manipulator operating 
within the abdominal cavity will be quantified. 

Previous work in this area has quantified surgical forces in terms of both su¬ 
ture knot holding strength and suture tightening, as well as in the quantifica¬ 
tion of the forces associated with tissue puncture with a scalpel blade or needle 
(EH dO] [EH (EE 071 SS [501 • Additionally, the abdominal organ tissues have been 
characterized m intro with rectangular samples [ 100| . However, these studies have 
limited applicability to clinical applications and the results of these studies show 
a wide range of variability in the measured forces due to the selected organ, the 
methodology, mid the inherent nonlinear biomechanical properties of soft biolog¬ 
ical tissues. Thus, a better understanding of the forces associated with specific 
surgical tasks Is necessary. Particularly, the forces associated with suturing, in¬ 
cluding suture needle tissue puncture and the force necessary to break the needle 
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away from the suture line, have not been well characterized previously. In addi¬ 
tion, little work has been doue to quantify the forces necessary for manipulation 
of the abdominal organs. The general assertion is that the biological tissue should 
be handled in such a way to minimize the potential for tissue damage Il0l[ . How¬ 
ever, t hese forces need to be quantified to facilitate the design mid development 
of surgical instruments. 

Though there are a variety of tasks common to surgical procedures, manipulation 
of the abdominal organs and suturing were selected for evaluation of their force 
requirements as part of the present study. These tasks were selected due to their 
widespread applicability to a variety of surgical procedures and their likelihood of 
having greater tool force requirements. Further, in the case of suturing, the force 
capability of the surgical robot can affect the quality of the knot and its ability 
to perform as intended PJ3]- The results of this study will help to define the force 
requirements for novel surgical tools for procedures involving suturing tasks and 
the manipulation of soft biological tissues. 


2.1.2 Methods 

In order to determine the tool force requirements, a Nano25 six-axis force/torque 
transducer (ATI Industrial Automation) was mounted between a handle and prob¬ 
ing tip to provide mi intuitive interface for force measurement that could be used 
in the operating room, as shown in Fig. |2.1[ Through a collaboration with The 
Visible Heart Laboratory (University of .Minnesota) |H)2| . access to a cadaveric 
porcine model was possible, which enabled measurement of the forces necessary 
for abdominal tissue manipulation. Abdominal tissues including the liver, spleen, 
gallbladder, stomach, kidney, and bladder were examined by using a force probe 
to manipulate the tissue such as may be necessary during surgical intervention. 
The force was recorder! ms a magnitude applied to the tip of the probe. Then the 
peak value was selected from the data set for each organ mid reported. In each 
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case, the organ was manipulated in order to gain access to the underlying tissue 
or surrounding organs, (’sue was token to avoid damaging the I issue by using the 
minimum force necessary to move the organ, in accordance with the practice of 
minimizing the potential for tissue damage os suggested in the literature |101| . 
As can he seen in Fig. |2.2| the abdomen was open during testing, which would 
he similar to the conditions during surgery when the abdomen is inflated. 



Figure 2.1: Illustration of the mounted six-axis transducer 


Using the same instrumented force probe, results were also collected for common 
suturing tasks such os tissue puncture and knot tightening, seen in Fig. |'2.3[ 
where high tool force requirements would he expected. Both porcine liver and 
gallbladder t issues were used for needle puncture force measurements. The suture 


needle (size 4-0. 1.5 metric) was attached to the force probe shown in Fig. 2.1 and 
was passed through t he tissue os would be done for a normal suturing procediue. 
The peak magnitude of the force on the tool tip was recorded. 


Next, suture knot tightening forces were measured using both size 0 (3.5 metric) 
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Figure 2.2: Force measurements being taken in a porcine abdomen 


Table 2.1: Suture sizes and types evaluated for break force 


U.S.P. Suture Size 

Metric Size 

Material 

Product 

Manufacturer 

1 

4 

polyhut ester 

Novafil 

Davis & Geek 

3-0 

2 

polypropylene 

Surgilene 

Davis \z Geek 

4-0 

2 

chromic gut 


Davis & Geek 

4-0 

1.5 

polypropylene 

Surgilene 

Davis .*c Geek 

5-0 

1 

polypropylene 

Surgilene 

Davis ,*c Geek 

6-0 

0.7 

braided polyester 

Polydek 

Deknatei 

6-0 

0.7 

polypropylene 

Surgilene 

Davis & Geek 

7-0 

0.5 

polypropylene 

Surgilene 

Davis & Geek 


and size 4-0 (1.5 metric) polypropylene (Surgilene, Davis <k* Geek) sutures. The 
suture lines were pulled, and the tool tip force measured, until a satisfactory knot 
was formed, ns shown in Fig. |2.3| This method introduces inherent variability 
between the individual surgeons tying the knot. However, by maintaining con¬ 
sistency within this study it: is possible to compare the knot tightening forces 
between suture sizes. 


Lastly, the forces necessary to break the suture line were cpiantified using a variety 
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Figure 2.3: Evaluation of suture knot tightening force requirements 

of suture sizes and materials ns shown in Table [271] The suture needle was rigidly 
damped as shown in Fig. |2.4| which was itself mounted to the crosshead of a 
tension testing machine. The free end of the suture was then secured to the 
bottom frame of the tension testing machine using a looped knot. Mount ing the 
suture in this way made it possible to evaluate the suture break force regardless of 
the location of the failure, whether the break occurred in the suture line or at: the 
junction between the line and the needle. The break force was measured using a 
250 N, uniaxial load cell (Fart Number 4501018/B) moiuited on a Qtest QT/10 
(MTS Systems Corporation) tension testing machine. Additionally, the crosshead 
of the tension testing machine was moved vertically at a rate of 0.1 meteis per 
minute until breakage occurred. Displacement rates of 0.05 and 0.2 meters per 
minute were also tested but the dependence of the break force on displacement 
rate was found to be insignificant. 
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Figure 2.4: Rigidly mounted suture for measuring break force 


2.1.3 Results 


As can be seen in Table |2 2| the greatest forces recorded during the abdominal 
force probing experiments were for manipulation of the kidney and liver tissues, 
which each demonstrated a measured load of approximately 3.5 X. This result 
would be expected due to the rigidity of the kidney and the mass of the liver, each 
of which contribute to the increased force necessary to manipulate the tissues. The 
remaining abdominal tissues tested fell within n range of 1.3 X to 2.5 X. 


Table 2.2: Abdominal probing force measurements 


Tissue 

Force (X) 

Liver 

3.45 

Gallbladder 

1.32 

Spleen 

1.47 

Stomach 

2.53 

Kidney 

3.57 

Bladder 

1.53 


Experimental data collected to evaluate the force requirements for the suturing 
tasks of tissue puncture and knot tightening is shown in Table |!P| where n denotes 
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the number of tests performed for each tiuslc. We can see that the tool force 
necessary for tissue puncture ranged front approximately 1.5 to 3.3 N, while the 
force necessary’ for suture knot tightening ranged from 2.3 to 7.2 N. 

Table 2.3: Sut ure task force measurements 


Task 

n 

Force (N) 

SD 

Liver puncture 

1 

3.26 


Gallbladder puncture 

3 

1.49 

0.72 

Size 4-0 knot tightening 

3 

2.28 

0.10 

Size 0 knot tightening 

2 

7.19 

0.96 


The forces necessary to achieve suture breakage are shown in Table |2.4[ where 
the suture material tested was polypropylene (Surgilenc. Davis «Vr (leek) unless 
otherwise noted. It is apparent that the force requirement to tighten the suture 
knot is dependent on the suture size. 


Table 2.4: Suture break force measurements. 


L.S.P. Suture Size 

u 

Force (N) 

SD 

1 (Polyhut. ester) 

4 

38.67 

5.70 

3-0 

6 

11.93 

0.15 

4-0 (Chromic Gut) 

2 

8.73 

1.71 

4-0 

8 

6.75 

0.80 

5-0 

2 

4.84 

0.51 

6-0 (Braided Polyester) 

1 

5.10 

— 

6-0 

3 

3.12 

0.56 

7-0 

2 

2.20 

0.34 


2.1.4 Discussion 


The abdominal tissue manipulation results presented in Table [iO| show that the 
greatest forces were measured for manipulatiou of the kidney and liver tissues. 










38 

This remit would be expected due to the rigidity of the kidney and the inass of 
the liver, each of which contribute to the increased force necessary to manipulate 
the tissues. Alt hough not directly comparable, tissue handling in terms of grasping 
force have been presented in the literature to have a mean value of 8.5 N [103| . 
Tims, the forces recording <luriug the abdominal organ manipulation experiment 
are expected to be within the safe range regarding the avoidance of tissue damage. 

Suturing task forces, presented in Table |?!3j were found to range between 1.49 and 
7.19 N. The lowest force measured was for needle puncture of gallbladder tissue. 
However, a large standard deviation of 0.72 N over three trials was found and was 
likely caused by variable tissue stiffness within the organ. The greatest force was 
meastued for suture knot tightening using sut ures of size 0. This result Is expected 
due to the increase in bending stiffness of a suture line with a linger diameter. 
Finally, it was possible to compare the knot tightening force of size 4-0 suture with 
the literature. The results demonstrate that there is a significant dependence on 
the suture size when considering the force requirement to tighten the suture knot. 
This corresponds well with the work of Kitagawa et al. in which the average 
tension applied to the suture by ;u» attending surgeon was slightly under 2 N for 
size 4-0 polypropylene JOT). Further, it was also shown that this force decreased 
with decreasing suture size. However, a definitive comparison between the two 
studies cannot be made due to dependence on the individual surgeon performing 
the operation. Lastly, the results for suture break force, seen in Table [27T| show a 
range of 2.20 to 38.67 N with a direct correlation with suture size. 

Based on these results, it is evident that there is a wide rauge of force requirements 
depending upon the task to be performed. Considering some examples of emerging 
surgical robots from the literature [351 C16) . which were shown to have a load 
capacity of 5 N, it can be seen that robots of this nature would be capable of 
performing many of the necessary tasks but would be limited to suture sizes of 
U.S.P. 5-0 or smaller. When comparing the results with the literature for size 4-0 
polypropylene, it was found that the force of approximately 6.75 X presented here 
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was much lower than the force of 13.70 N presented by Trail ct al. 1104) . One 
possible explanation for this result: is the location of breakage. In the experiment 
conducted by Trail et al.. the suture was secured by winding the line mound 
two metal rods before taping it down. This would dictate that the suture failure 
must occur in the line. However for this study, the needle end of the suture was 
secured in a vice and it was ohserved that the majority of the failures occurred 
at the junction between the needle and the line. This would suggest: that this 
junction is likely weaker than the tensile strength of the line. Since the focus of 
this study was to evaluate the surgical tool force requirements, and in particular 
for this experiment the force necessary to break the needle away from the rest 
of the suture, these results hold value. If it were determined that greater force 
was necessary to achieve a satisfactory knot, the surgeon would have the option 
of grasping the suture line in order to apply greater tension. 


2.1.5 Conclusions 

The results presented in Section |2.i.3| demonstrate that it is necessary for the 
designers of novel surgical robotics to consider the specific applications that their 
tool is targeting when determining the limits of capability. In a scenario such as 
this where increasing the force capacity of the tool may be size or cost prohibitive, 
it may be necessary to make use of specialized tools such as cutting instruments to 
achieve the forces that may be required for suturing tasks. The results presented 
here were developed under the assumption that the greatest forces required of 
a surgical robot would be found while performing tissue manipulation, suturing 
tasks, or suture needle separation. Under these conditions, the range of surgical 
tool forces was found to be e [1.32.V, 38.G7 jV]. It is apparent that even for this 
relatively small sub-set of common surgical tasks, there is much variability in the 
necessary tool load capacities. Therefore, based on the results presented here and 
the work of others found in the literature, a load capacity of 5 N is suggested as 
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a target goal as an end load on a manipulator for use in the design of endoscopic 
surgical tooLs used for abdominal procedures for which access is gained through a 
natural orifice as it provides a good trade-off between tool capability and the size 
of the tool. 


2.2 Workspace 

The workspace at the procedure location within the patient that the surgical ma¬ 
nipulator intLst be able to maneuver around is highly dependent on the procedure 
being performed. Further, the skill level of the surgeon will have a significant 
effect on the amount of space required to perform the necessary manipulations 
for a given procedure. Data collected using a set of instnunented laparoscopic 
surgical tools show’s that an expert surgeon can perform an example task within 
the confines of a cube of side lengt h 5 cm where a novice may require a cube with 
a side length of 8 or 0 cm to perform the same task | 1 U.~i) . Similar requirements 
were found for a surgical instrument for surgery on the throat |58j . 

Another requirement previously mentioned is the necessity for the surgical instru¬ 
ment to achieve 180° retroflexion. This requirement allows lor tlie surgeon to 
access surgical sites in t he upper abdomen when entering through t he esoplnigus. 
However, the burden of meeting this requirement does not necessarily have to 
fall solely on the individual surgical manipulator, but instead on the endoscopic 
platform as a whole. This means that if retroflexion can be achieved by the endo¬ 
scope that is delivering the manipulators to the surgical site, then the individual 
manipulators do not need to have the ability to retroflex. Assessment of common 
surgical practice has found that when using laparoscopic tools, the surgeon will 
spend 95% of the surgery within a cone with a vertex .uigle of 60° |10(>| . Tints, 
it is possible to trade bending range of motion for tool size or load capacity if 
necessary for a surgical manipulator design to he successful. 
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Putting these two constraints together, a workspace with an edge length of 5 cin 
is a reconmieiuled target for the positional range of the manipulator spanning one 
end of the range limit to the other. Further, the manipulator should be capable of 
operating through an angular range of 60° at a minimum for the majority of surgi¬ 
cal manipulations. Under these conditions it is desirable for the manipulator to be 
capable of spanning the workspace within a 60° bend/orientation range. However, 
under these conditions, a combination of manipulator bending and longitudinal 
positioning of the endoscope would need to be used to reach the proximal and 
distal limits of the workspace. Additional functionality in terms of the orientation 
range of the manipulator would reduce the necessity for longitudinal repositioning 
ol the endoscope. The constraints suggested here preclude the use of excessively 
short length manipulators as well as manipulators with a small angular range of 
motion and long length. 


2.3 Safety Considerations 

The use of hydraulics in surgical applicat ions fact’s some obstacles that need to be 
considered and overcome for this modality to be successful. One common theme 
in the lield of hydraulics is that leakage of the hydraulic medium is difficult to 
prevent and must be expected. However, in the surgical application, there is little 
room for uncertainty when it comes to possible contamination of sensitive areas 
such as the peritoneal cavity IHI7 |. This means that the hydraulic medium, such 
as saline, would need to be subject to the same sterility recpiirements as the rest 
of the device. 

An additional concern stemming from the leakage of hydraulic fluid is any possible 
damage to nearby biological tissue due to the high operating pressures. Due to this 
high pressure, the hydraulic fluid would have the potential to escape the tool at 
high velocity and possibly cause damage by compression of sensitive tissues found 
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within the abdominal cavity, as shown in Fig. |2.5| Indeed, the use of water jets 
have been proposed in the past as a means of tissue resection |10S| . Therefore, it: 
would be necessary to take this into consideration when designing the flow control 
in order to include safeguards against this sort of undesired behavior. 



Figure 2.5: Example of tissue damage due to applied pressure [Iflflj . 

The literature Inis shown t hat it can be difficult to predict whether or not: a given 
load will result in damage to the tissue and further that the magnitude of the load 
that may cause damage varies significant lv depending on what kind of tissue it is 
[H)fllllllHllll| . However, using these studies as general guidelines, it is possible to 
say that tissue damage begins to occur in the range of 160 to 2S0 kPa of applied 
pressure [1001 (110) . 

To find the pressure exerted on a tissue by a fluid leak it is assumed that the 
force applied to the tissue by the fluid results from a complete loss of forward 
momentum when the fluid impacts the tissue such that the force Is calculated as 

F = ni Av (2.1) 

where m' is the mass How rate of the fluid and Ac is the change in the forwiird 
velocity of the fluid. Then assuming that the velocity in the direction normal to 
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the tissue surface is zero after impact, the equation can he written as 

F = pAv 2 (2.2) 

where p is the fluid density. A is the area of the hole, and v is the velocity of 
the fluid as it exits through the hole. This relationship assumes that the fluid jet 
from the leak impacts the tissue normal to its surface and that the tissue is close 
enough to the site of the leak that the jet does uot expand or lose velocity. If 
the velocity of the fluid is approximated using the orifice equation, the resulting 
pressure applied to the tissue at the point of impact is calculated as 

l'=- x = P'' 2 = 2 C 2 d P, (2.3) 

where C'j is the discharge coefficient of the leak orifice and /’, is the pressure 
internal to the system at the leak site. Using an approximate discharge coefficient 
of O.G. this suggests that the pressure applied to the tissue could be as much 
as 72% of the supply pressure in the worst case scenario. It is also important, 
however, to consider the time over which the tissue is exposed to this pressure. 
Rodrigues et ni. )111| showed that tissue subjected to sustained load from a suture 
began to rupture within approximately 15 to 20 seconds when the magnitude of 
the load was sufficiently high. Therefore, one possible method of reducing the 
risk of tissue damage would be to actively monitor for leakage and turn off the 
supply if leakage is detected, thus minimizing the amount of time that the tissue 
is exposed to unintended loads and lowering the risk of permanent tissue damage. 

This safety concern could also be addressed through the introduction of a How 
resistance within the supply line such that the downstream pressure drops signif¬ 
icantly at high flow rates. For example, with the inclusion of an upstream How 
restriction, the pressure drop across the restriction could be calculated ns 



(2.4) 
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where the area of the How orifice. A, would be sized such that the downstream 
pressure drops to a safer value when the flow rate. Q, through the orifice signifi¬ 
cantly exceeds the nominal flow rate of the system under normal conditions. 


2.4 Tool Size Constraints 


Another consideration for a robotic instrument to work in a surgical setting is 
that it must conform to the size constraints placed upon it by the anatomy of 
the patient. For natural orifice surgery, these constraints come primarily from the 
access pathways as they are significantly more restrictive than the available space 
within flte abdominal cavity. In the case of transgastric access, the overall endo¬ 
scope is constrained to a diameter in the range of IS to 22 mm Within this 
cross-section must, be contained vision channels (preferably stereoscopic), suction, 
irrigation, and at least two tool channels. An example layout of the endoscope 
cross-section Ls shown in Fig. 27fi This layout permits the tool channels to have a 
diameter of approximately (i mm as compared against a 20 mm endoscope. Within 
these tools channels, the components of the manipulator must be contained. For 
a two-segment manipulator, this would include two valves and six actuators. The 
actuators can be arranged serially as necessary to meet the overall size constraint. 
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Figure 2.G: Example layout of endoscope cross-section with necessary channels. 
The manipulators described here would be contained within the tool channels. 

2.5 Summary of Design Requirements 


The considerations necessary for the design of a surgical manipulator are summa¬ 
rized here and separated into four primary areas. These considerations include 
the load capacity of the tool, the size of the workspace, safety of operation, and 
biological constraints on the overall size of the device. 


First, considering the tool load requirements as laid out in Section 27T it is ap¬ 
parent that the tasks that a surgical tool may he expected to perform can vary 
widely. Further, these tasks and the associated load requirements me likely to 
dilfer greatly between the types of possible procedures. However, based on the 
conditions specific to abdominal procedures and taking into consideration the lim¬ 
itations of scale, it was suggested that a force load capacity of 5 Ts he the target 
load capacity for a surgical manipulator. 


Next, the size of the workspace for t he manipulator was examined in Sect ion 12,2 


through a review' of the data from the literature collected for laparoscopic tooLs. 
From this information it was gathered that the tool workspace should he on the 
order of 50 mm in all three Cartesian directions to accommodate common surgical 
tasks such as suturing, cutting, or manipulation. Additionally, it: was found that 
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the majority of procedure manipulations take place within a cone rotated around 
the long axis of tlie tool with a vertex angle of 60°. This realization introduces 
some flexibility into the workspace requirements as it is now possible to transfer 
some of the manipulation range to the endoscope, which is used to deliver the 
manipulators to the surgical site, since it would not be unreasonable to re-position 
the endoscope between surgical tasks. 


In considering the safe operation of a hydraulically powered manipulator within 
the abdominal cavity, there are several potential unique issues that must be in¬ 
cluded as discussed in Section 2.3[ Of primary importance is the possibility of 
leakage which inherently necessitates the usage of a sterile hydraulic fluid such 
as standard saline. Additionally, in the event of leakage, it must be ensured that 
the leaking fluid will not result in unintended damage to the surrounding tissues. 
There are a few options for mitigating the risk of tissue damage. First, the How 
rate of tire fluid could be monitored and if a leak is indicated by a sharp unex¬ 
pected increase in flow rate then the supply pump is switched off. Second, mi 
upstream flow resistance could be used to induce a pressure drop if the flow rate 
should drastically increase due to a leak, thus reduciug the risk of tissue due to 
exposure to high fluid pressure. 


Finally, the tool size must be designed such that the instrument Is biologically 
compatible with the patient and the intended access pathway. In the application 
to an esophageal access pathway, tins limitation Is on the order of 18-22 nun in 
diameter for the endoscope as a whole. Thus, considering the space available to the 
tool channels which would contain the surgical manipulators, the size limitation 
is likely to be constrained to 6 mm in diameter which must house all of the 
components necessary to operate the manipulator. 





Chapter 3 


Manipulator Backbone 


This chapter is the first of four that describe the design of the manipulator, each 
focusing on one particular aspect of the overall design. The focus of this chapter is 
the backbone (support) structure for the manipulator. First. Section |Q| presents 
a discussion of the design choices including the selection of a continuum beam 
structure as the best fit for this design. In Sections and 373 the methods 
of modeling a continuum beam are discussed to provide a reasonable simulation 
of the backbone component that will be revisited when considering the system 
assembly. Experimental results ore then presented in Section pTl] to provide a 
measure of the quality of the model under the prescribed assumptions. 


3.1 Motivation and Related Work 

There are a variety of technologies available to serve as the structure for a surgical 
robot. The most familiar robot form factor seen in industry Is comprised of a 
small number of rigid links arranged in series. However, this style of robot is 
not suited for a natural orifice access technique as small motions can require large 
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articulations and thus safety is obtained by providing the robot with room to move 
|47| . The tight space requirements found in minimally invasive surgery discussed 


in Section 2^2| necessitate the use of a robotic structure with compliance and a 
high ratio of length to diameter. In this format, the primary options are hyper- 
redundant robots with a high number of discrete links or robots with a continuum 
struc ture. Discrete link robots have the advantage of being more precise due to 
hitting a finite number of degrees of freedom, each of which can be actuated. 
However, for this style of robot to be feasible on a small scale, the geometry 
of the links and their connections must be relatively simple such as one degree 
of freedom pin joints. For a pin joint it would then he necessary to use two 
tendons/actuated inputs to control each degree of freedom thus increasing the 
robot complexity. Alternatively, continuum robots offer safe operation through 
natural compliance and ability to adapt to an unpredictable environment [TTIRTI] . 
This Is due to the presence of an infinite number of degrees of freedom for the 
continuum robot while only a few of those degrees of freedom can be actuated. The 
downside to this Is that the lmnctuated degrees of freedom passively determine 
the resulting beam shape under a given load which can make precise positioning 
difficult. As discussed in more detail in Section G.2.2[ for this application three 
actuators provide control over three degrees of freedom per manipulator segment. 
Thus in tin* example case of achieving six total degrees of freedom, six actuators 
could be used compared with twelve actuators for the example pin-jointed rigid 
link case. It is then the remaining, unactuatcd degrees of freedom that passively 
determine the beam shape. Further, the inherent compliance of the continuum 
robot mokes il an attractive option for medical applications since the robot will 
conform to the unpredictable environment |112j . For this design, the continuum 
robot solution was selected to take advantage of its compliance and greater ease 
of miniaturization through simplicity. 


The use of a continuous backbone for cable driven robots has been well docu¬ 
mented in the literature [TTfl IS71 (5611%71 [88l 1113[ . Continuum manipulators have 
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many advantages that make them an ideal selection for this application including 
simplicity and dexterity |8ti| . By using a beam material with desirable high strain 
behavior such us nickel titanium alloys, it is possible to reach a wide range of po¬ 
sitions as predicted using beam bending theory. This style of joint would include 
the flexible beam at the center with a variable number of disks attached to allow 
the tensile bin's to be routed appropriately as shown in Fig. |3. 1[ This approach 
Is similar to the hvper-reduudant robots consisting of a multitude of rigid links, 
previously discussed, providing comparable yet distinct operating characteristics 


W1 



Figure 3.1: Graphic of beam bending based design. 


Unfortunately, the prediction of performance for the continuous backbone manip¬ 
ulator is noil-trivial and there have been numerous models developed to predict 
this behavior as discussed in Section |1.1.5| [ 551 [6411671 [65] . However, many of 
them require significant computational effort and thus make it difficult to allow 
for practical application or they do not capture enough information to enable 
integration of the model into the model of a robotic system as a whole. There¬ 
fore. a discussion of the modeling of a continuous backbone for this application is 
presented here. 


3.2 Simplified Euler-Bernoulli Beam 


To enable modeling of beam behavior within the system, it Is first necessary to 
decide on a method of describing the mechanics of the continuous beam. The task 







here is to allow for rapid prediction of manipulator behavior under the conditions 
of a cable driven continuum manipulator. If a model can be specified with reason¬ 
able accuracy, it can then be incorporated into a system model of the manipulator 
us a whole. One of the more common ways to capture the large-deflection behav¬ 
ior of the continuous backbone is to take the Euler-Bernoulli Theory and apply 
it to the special case of bending in a beam subjected to pure moment loading. 
The approximation of a tendon-driven robot as a beam subject to pure moment 
loading has beeu experimentally validated assuming proper selection of spacing 
and height: of the tendon support disks [S71II15) . By neglecting arbitrary external 
loading, the solution to the Eider-Bernoulli Theory becomes trivial. While the as¬ 
sumptions of this model are overly constraining, this solution method is included 
here to provide insight into the predicted beam behavior under the strict condi¬ 
tions of pure moment load. Therefore, under the assumptions of uniform material 
properties, pure moment tendon driven loading, negligible shear deformation, and 
negligible external force loading, we start with the total potential energy function 
for the Euler-Bernoulli beam | llt>) 




(Is 


(3-1) 


where the first: term represents the distortion energy in the beam and the re¬ 
maining terms represent the external work of the applied moment load. Here it 
is assumed that the external loading is known and thus the energy contributed 
by the external load can he accounted for. Further, 0(s) represents the angle of 
the tangent to the curve of the beam at any location along its length measured 
relative to the undefonned configuration, mj is the moment being applied at the 
free end of the beam (see Fig. |3.2t l. and E is the modulus of elasticity for the 
backbone material. In these equations and throughout the remainder of this dis¬ 
sertation except where otherwise noted, the prime denotes the first derivative of 
the variable with respect to the position along the beam, s 6 [0, L ] where L is the 
length of the beam. 
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Now to minimize the total potential energy by minimizing the integrand, we choose 
0(.s) such that 0(11) - 0„ arid ff(s) = AC minimizes - m L 0'. This yields AC = 
which then provides the unique uiinimixer 0 (.s) that satisfies the boundary 
conditions as 


l »-5r* + 


< 3 - 2 ) 


Using this solution, the resulting shape of the beam con be expressed as a circle 
with a radius of 

= — (3.3) 

niL 

as depicted in Fig. |3.2| which is consistent with the basic constitutive relationship 
for the Euler-Bernoulli beam [TT71 ITTSj . 



Figure 3.2: Resulting beam geometry from an end applied moment. 


The resulting radius of curvature as calculated in Eq. |3.3| can be used to find the 
tool tip location by expressing the coordinates of t his point in both the longitudinal 


(z) and transverse (x) axes, as shown in Fig. 3.3 These equations are found to 
be, 

' ' ' (3.4) 


-mi/ = '7,ir,, sin 


r 6m, 


' criif ' 6mm 


1 — COS 


-±~) 

' beam / _ 


(3.5) 


where / is the length of the beam segment. This result suggests that it is possible 
to achieve a smooth motion of t he tool tip when using a structure based ou beam 
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bending principles and a controlling force applied as a moment. This result has 
been confirmed in the literature where it has been shown experimentally that a 
cantilever beam subject to moment loading assumes the shape of a circular arc in 
the absence of external loading JS7J. 





Figure 3.3: Graphic showing location of the beam tip. 

The limitation of this forward kinematic model is that it only accounts for mo¬ 
ment loading at the tip. Any forces or other external loading applied cannot: be 
accounted for in the same simple manner. For this reason it is necessary to ap¬ 
proach the problem from a different direction that allows for greater flexibility in 
terms of loading scenarios and more realistic usage cases. 
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3.3 Model Prediction of Beam Shape and End 
Loads 

An alternative approach to predicting the behavior of a continuous backbone will 
be presented in this section. This model makes use of the principles of inverse 
kinematics and potential energy minimization to predict beam shape and end load 
conditions for a specified end position and orientation of a given beam segment in 
a manner similar to the classical Cosserat rod model |115) . The model begins by 
assuming that the backbone of each manipulator segment is modeled as an Euler- 
Bernoulli beam which Ls cantilevered at the base, similar to what was previously 
discussed. However, in this case the externally applied loading Ls not known 
and thus is not factored in to the total energy. Assuming that the position and 
orientation of the distal end is given, we wish to determine the beam shape and 
the beam reaction force and moment acting on the end-effector at the distal end 
of the beam. 

The method of modeling the deformation of a robot with a continuous backbone 
which will be presented here is iui energy-based formulation |110| . The basis of 
this method is that of all the beam shapes that satisfy the specified end point 
conditions, the continuous backbone with negligible mass will assume the shape 
with the least stored elastic deformation energy. Knowing this shape for the 
specified locat ion and orientation of the distal end of the beam, it is then possible 
to calculate the resulting load applied on the terminating disk located at. the 
distal end of the beam which can include a force vector and a moment vector of 
arbitrary directions. This inverse kinematics, quasi-static model, which neglects 
the effects of the mass of the beam itself [Sfill?$S) . is here referred to as the Inverse 
Formulation Model (IFM). 


3.3.1 IFM in Two-Dimensions 


In the two-dimensional case, the problem simplifies to a cantilevered beam with 
arbitrary position (i and x) and orientation (0) end conditions as shown in Fig. 
|3.4[ Under these end conditions, the beam will have some shape which, if known, 
provides additional informat ion about what conditions (load and/or moment) are 
present at the distal end. 



Figure 3.4: Graphic showing the condit ions of a cantilevered beam with arbitrary 
end coordinates. 

The IFM in two-dimensions assumes planar bending only and thus can only be 
used to describe in-plane external loading. Let * 6 [0, L\ denote the distance 
along the beam of length L. and 0{s) be the angular orientation at location .s. 
Thus, 0(a) is used to describe the shape of the beam for n given set of position 
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iuui orientation coordinates for the distnl end. In other words. I he beam shape 
is specified as a n-degree parameterization of O(-) with respect to s, assuming 
that the location and orientation of the distal end of the beam ore specified as 
inputs such that x cn j, z in j. and 0 rn ,i are known. The parameterization of 0(-) can 
take many forms; however, the use of a polynomial function and a spline function 
iire evaluated here. For a selected parameterization method, the beam shape is 
computed using the potential energy stored within the beam. 


Potential energy is stored within the beam of cross-sectional area A. elastic mod¬ 
ulus E. second moment of area /. and shear modulus G\ as elastic deformation 
consisting of axial, bending, and shear components with stiffness terms of EA. 
EL and (7.4, respectively [ITTf] . This formulation differs from what was presented 
in Eq. iTT in that it includes terms for axial and shear strain and does not include 
terms for external work since the external loading is not known a priori and thus 
only quantifies stored energy. Using these definitions, the potential energy within 
the beam can be expressed as 


PEW- ] ) = \J [iB(s) T £3MCQ(») T B(- 


ElO'l 


r/.s 


(3.6) 


where the prime denotes the first derivative of the variable with respect to the 
position along the beam, s ns before and <V(.s), C. •?(,«). and Jr(.s) are defined as 


Q(«) = 


c = 


cos(0(s)) -Sin(0(*)) 
siu<0(s)) cos(0(s)) 

EA 0 

(» a a 


*(.,)'-cos(0(a)) 
x{sy - sin(0(s)) 

*(*) ' 

*(*) . 


(3.7) 

(3-8) 

(3.9) 

(3.10) 
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By these definitions, (^(.s) represents the crass-sectional orientation as a function 
of s, C is used to collect the axial and shear stiffness terms, and v(s) measures 
axial iuhI shear stretch of the beam. The vector x{s) represents the position vector 
of the beam center line as a function of .v 

Using this definition of the potential energy, the beam shape is predicted using a 
constrained optimization to minimize the potential energy stored within the beam 
subject to the boundary conditions at a* = L. 

x(L)=x tnd . z(L) = z rnd , 0{L)=0 rnd (3.11) 

and at a — 0. the cantilevered beam boundary conditions 

i(0) = 0. 3(0) = 0, 0(0) = 0. (3.12) 


Evaluating the energy function through the application of Hamilton's principle 
j!20| as discussed in the literature results in the following partial differential equa¬ 
tions governing the beam, adjusted for the static case |11P| 

[g(s)C^ r (s)P(s)]' = 0 (3.13) 

£70"(.s) + ^(*)SQ(s)CQ T (s)t7(s) = 0 (3.14) 


where S is the skew-symmetric matrix 



0 1 

-1 0 


(3.15) 
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Those differential equations are further subject to the following, boundary condi¬ 
tions. 


E10'{L) = Mb 

(3-16) 

QCQ T v\ L = K 

(3.17) 

0(0) = 0 

(3*18) 

m =o. 

(3.19) 


These boundary conditions designate that the proximal end of the beam is located 
at the origin, is clamped at s = 0. and that the end loads are defined as fj, and 
Mb with respect to the global coordinate system at .s = L. 


Examination of this model reveals that the majority of the deformation of the 
beam will result from bending deformations due to lower energy requirements. For 
example, for the beam used experimentally in Set tion [iO| the stiffness quantities 
are El = 5.6361 • 10" 5 N’r/r. EA = 4453.2 N. and GA = 1590.4 N. Therefore, it is 
possible to save significant computational effort by reducing the potential energy 
function to the classical model form which accounts for only energy from bending 
deformations as ^ 

pm-)) = \ f EItr{sfds (3.20) 

A justification of this decision is presented in Section |3.4[ With a given solution 
based on the minimization of the strain energy, the beam shape can be specified 
using the positional coordinates of the centerline, x(.s-) and ;(.s), which, since the 
relative position of the end point of a differential element of the beam at s is 
[srn(0(s)); cos(0(s))] fits, can be given as functions of 9(s) such that 


: 


{,.)= jTsin (B[a))da , s€(0.L| (3.21) 

(*) = f o caa(B(tT))<I<T , s e [0, L ] (3.22) 
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An optimization algorithm is used to minimize the stored strain energy within the 
beam in order to find the lowest energy beam shape corresponding to a set of end 
conditions. This is accomplished by parameterizing 0(a) as a high order polyno¬ 
mial and using a constrained optimization to minimize the potential energy (Eq. 


■1.20ft . subject to the constraints of Eqs. |3.11 and 3.12 within the beam as cal¬ 


culated above over all possible combinations of the parameterization coefficients. 
An example of the output from this model is shown in Fig. 13.5| for which the 


bending stiffness. El. was 5.0351 - 10~“ N-wt 2 . Noting, however, that the beam 
shape prediction is independent of El assuming that both are constant along the 
beam. Inputs to the model were 30 mm for both the transverse and longitudinal 
tip positions and mi angle of 10° for the tip orientation. 



Figure 3.5: Example backbone model output for arbitrary input parameters. 


If boundary conditions corresponding to the constant curvature load case are 
applied, it is found that the 2D IFM output using a 0th order polynomial pa¬ 
rameterization for 0(a) matches the beam shape predicted by the basic constant 
curvature model. This is demonstrated in Fig. |3.0[ where the beam shape from 














59 

each model is plotted for a variety of constant curvature input conditions with 
0(L) 6 15°. 150°]. As con he seen, the curves from eacli model perfectly over¬ 
lap for nil inputs. Further, comparing the model predicted end loads from the 
1FM with that of the constant curvature model it was found that they compare 
well. Moment loads matched the loads specified by Eq. and force loads were 
zero, noting however that small deviations in position and orientation occasionally 
result: in huge changes in t he predicted end load. 



Figure 3.6: Comparison of 1FM output with that of coils taut curvature model for 
coiLstaut curvature input conditions. 


Once the potential energy of the beam has been minimized using the optimization 
algorithm, the beam shape is known and the end loads can be calculated. This is 
done using the principle of virtual work which states that the load corresponding 
to a given displacement can be computed if the associated potential energy is 
known [118| . Under this principle, the three relevant end loads. (Fft_p.Fi,-, A/fc), 









can be calculated as 
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dPE' OPE' OPE' 

f'b.l - -• I'b.; -- -*'(> — "TTTJ- 

OXend OZend d 

where PE * is the potential energy stored within the beam as calculated by the op¬ 
timisation. These values can be computed directly using the Lagrange multipliers 
associat ed with each of the end position and orientation constraints. 



3.3.2 IFM in Three-Dimensions 


Similar to the two-dimensional model, the IFM in three-dimensions is based on the 
description of the beam as a smooth parameterized curve. However, for the three- 
dimensional model, the beam shape must be described by three successive rotation 
angles. 0(.v), and &{*). which are Euler angles representing the orientation 
of the beam cross-section with respect to the x, i/. mid axes (the rotated 
local coordinate axes), respectively, as shown in Fig. |3.7| By this definition, the 
reference configuration of the beam originates at the origin and is tangent to the 
Z-axis of the reference coordinate system. 


The rotation matrix. /?(.*>) 6 £0(3), used to represent this series of rotations about 
the local coordinate system contains the intrinsic Euler rotations about the r-axis 
by $(.*»), then about the (/-axis (new v-axis) by 0(a). and finally about the ;"-axis 
(new z-axis) bv t’(s) (see Fig. |3.7[ and is realized as shown in Eq. 3.24 The 
prime notation on the axis labels indicates that the axis has changed orientation 
as a result of previous rotations. 








til 


COS 0 COS lj> 

cos o sin ip + sin 0sin If cos t ■ 
sin 0sin r — cos 0sin If cos tj> 


— cos 0 sin ip 

cos 0 cos ip — sin dsin If sin ip 
sin ocas v + cos ©sin 0 sin i' 


sin 0 

— sin ip cos If 
cos o cos I) 

(3-24) 



Figure 3.7: Three dimensional Cartesian coordinate system showing orientation 
angles for a differential beam segment, At. The y'-axis represents the new, rotated 
v-nxis and the ;"-nxis represents the new rotated z-axis. 


Using the rotation matrix, it is possible to find the skew symmetric matrix which 
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contains the two bending components as well as the torsion component that to¬ 
gether describe the beam shape mid orientation as 


,,(s) = R(s) T R(s)' = 


0 My 

U. 0 -Mr 

a r 0 


(3-25) 


where u, and u y relate to bending about the local r ami y axes, respectively, mid 
u : relates to torsion about the local z-axis [ 1 l r >[ . Evaluating this expression mid 
pulling out the functions for u r , and u. results in 





sin tf>(s) 0 (a)' + cos cos B(s)4>(sy 

n{s) = 

"v 

= 

cos i’(.s)0(.s)' - 0060 (a) sin i'{.s) 0 (.s)' 


Ms 


sin 0 (s)<f>(sy +• 


(3.2G) 


Similarly, the differential expression of x(,v), y(s), and :(.v). which define the 
coordinates of the beam centerline, can be specified as 


&(*) _ 
ds 


*'00 
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(3-27) 


Further, the resultant beam curvature, AC(.v). Ls calculated as the magnitude of the 
arc length derivative of the tangent vector where the tangent vector of the curve 
at any location Ls found as r(.s) = .r'(.s)' 112111122| . Evaluating these expressions 
for the rotation matrix from Eq. |3.2-1| yields the following expression for beam 
curvature 

IC(s) = y]«l + uj - v/u-W 1 + *-(») W#(») (3.23) 

with the values (f(s) and o'(s) defined as the differential change in 0 (s) and o(s) 
over the segment, respectively. Similarly, the torsion. r(.s), comes directly from 
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r(.s) =an0(s)0(a)' + v>(«O' 


( 329 ) 


Thus, for tlie* deformed beam which traces some curved path defined by the above 
parameterization, t he potential energy stored within the beam can now be ex¬ 
pressed as a function of the beam curvature and torsion as 

PEM-).B(-)M-)) =^j L [EIK(s) 1 + GJt{,)*\ ,1, (3.3D) 


where El is the bending stiffness and GJ is the torsional stiffness noting that 
./ = 21 for a differential segment of a circular cross-section. This function for 
the potential energy is of the same form as was shown in Eq. |3.20| with the 
addition of torsional energy, but still neglecting axial extension and direct shear 
which are typically negligible |11!)| . This formulation of potential energy can 
be equivalently represented using the components of the skew symmetric matrix 
directly. t7(s) = [u r , u M , u ; | r . when the expression ii(s ) 1 ■ Kbt * t?(/») is substituted 
for the integrand in Eq. |3.3l i| Here A'bt is the stiffness matrix represented as 



El () () 

0 El 0 
0 () GJ 


(3-31) 


Again, assuming that the location and orientation of the distal end of the beam 
are specified as inputs such that x rlld . y rnd . z rnd , <t> rnd . O c ,ui- and i£> rnd are known, 
and integrating x(s ) 1 to get the vector components delining the beam centerline, 




*(»). .»/(«)• and :(.s) as 
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*W = j 

f fun{0{o))(lo , 

.s € (0. L] 

(3.32) 

y(*) = t 

[ - COS(0(fT))sin(0(lT))f/j7 , 

*6[0 ,L) 

(3.33) 

=(.,) =. 

f cos(0|<r))cos(c>(«7))(fcr , 

s 6 [0, L], 

(3-34) 


0 


the three-dimensional beam shape is then given by the constrained optimization 
of the potential energy of the beam subject to boundary conditions at s = L. 

j(L) = x rnd , y{L) = i/ rnd , z(L) = z rMd . (3.35) 

<P(L) = <S> rnd , 0(L) = 0 cnd , ii>{L) = i> cnd (3.36) 

and at s = 0 

2(0) = 0, 0(0) = 0, 0(0) = 0, 0(0) = 0 (3.37) 


As in the two-dimensional case, once the potential energy of the beam lias been 
minimized, the beam shape is known and the end load can be calculated using 
the principle of virtual work, realized as 


n, = 


i)PE- 


dx rnd ' 




MUL) = 


c)PE‘ 

<M„,i 


MuU) = 


dPE- 

% rtd 

dPE * 


n, = 


OPE' 


dO 


dz tnd 

MUD - 9PE ‘ 


mil 


ck 


mil 


(3.33) 

(3.39) 


where again PE' is the potential energy stored within the beam as calculated by 


the optimization. The moments calculated in Eq. |3.39| are with respect to the 
local rotated coordinate axes. 
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An example of the output from this model Ls shown in Fig. )3.8| where 




20 mm 

Vend 


10 mm 

-end 


35mm 

•Pend 


5° 

Mend 


10 ° 

tend 


15° 


(3.4(1) 


For this example, the orientation angles were parameterized using 6th order poly¬ 


nomials as in Section |3.3.11 and the material and geometric properties of the 
beam were E = 28.0 CPa, G = lO.SGPa. L = 50mm, / = 0.002mm 1 , and 
,/ = 0.004mm 4 . Then using Eqs. |3.38| and [3.39[ the calculated end loads are 


F, 


0.43jV 

n 


0.22jV 

F; 


0.91 A* 

M 0 


O.OOlGiV. m 

Mo 


0.0061 N • m 

M v 


-0.0007 N • m 


(3.41, 

















66 



Figure 3.8: Example of three-dimensional 1FM output for arbitrary input param¬ 
eters. 


3.4 Experimental Model Validation 

To test the effectiveness of the inverse formulation model at predicting the beam 
shape for a given loading, an experimental setup was prepared consisting of a 
flexible nickel titanium alloy beam (E - 28.0 GPa,G = 10 .SGPa) 11231 . support 
structure, and n six-axis load cell (Nano25, ATI Industrial Automation). The 
length of the beam was 50 mm for all test conditions with the exception of Sub- 
Fig. reference (b) (Table |3T{ |, which had a length of 55 mm, and a circular cross 
section with a diameter of 0.45 mm. Position and orientation of the distal end 
of the beam as well as the beam shape was recorded photographically (see Fig. 
|3.9[ | and a grid backdrop was used to provide scale. A variety of combinations of 
position and orientation were selected randomly, the values for which are shown 
in Table |3,1| For each combination, the beam was allowed to assume a stable 
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configuration and the beam was not observed to pass through a bifurcation info 
an unstable configuration. 



Figure 3.9: Example photograph of experimental setup used to record beam shape 
for arbitrary position and orientation of distal end. 


Table 3.1: Distal end location mid orientation conditions arranged according to 
corresponding sub-figure references. 


Sub- 

Fig. 

-end 

(mm) 

•t’rnrf 

(mm) 

&cnd 

(°) 

(a) 

39.7 ± 0.7 

27.3 ± 0.7 

0 ± 2 

(b) 

28.0 ± 1.0 

35.0 ± 1.0 

115 ± 2 

(c) 

33.0 ± 0.4 

25.0 ± 0.4 

45 ± 2 

<‘U 

-10.1 ± 0.4 

27.0 ± 0.4 

60 ± 2 

(e) 

34.9 ± 0.5 

30.0 ± 0.5 

-2.5 ± 2 

(f) 

46.6 ± 0.4 

7.4 ± 0.4 

-60 ± 2 


For each combination of position and orientation, the beam shape was recorded 
using image pixels to measure beam location and transferred to a data file. Next, 
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the same set of position and orientation conditions were used in the inverse for¬ 


mulation model described in Section 13.31 

Table 3.2: Model predicted end loads using both the full potential energy function 
(Eq. |3.0fr and the siinpliiied energy function (E<j. |3.20| l. Prediction of ston'd po¬ 
tential energy, PE' . also included for comparison. Forces are given in N. moments 
and potential energy in N - m. 


Sul>- 

Fig- 

Simplified PE Function 

F : F r My PE' 

F : 

Full PE Fund ion 

F, My PE' 

(n> 

0.13 

-0.11 

0.0039 

0.0025 

0.14 

0.13 

-0.0021 

0.0022 

<b) 

0.02 

0.07 

-0.0020 

0.0022 

0.02 

0.05 

-0.0017 

0.0019 

(<-') 

0.57 

0.38 

0.0054 

0.0050 

0.51 

0.36 

0.0048 

0.0046 

(<1) 

0.01 

-0.02 

-0.0007 

0.0006 

0.01 

-0.02 

-0.0010 

0.0006 

(e) 

0.27 

-0.03 

0.0045 

0.0037 

0.25 

-0.03 

0.0051 

0.0034 

(0 

0.18 

-0.24 

0.0000 

0.0030 

0.20 

-0.31 

0.0084 

0.0032 


First, comparing the predictions made using the planar 2D models for the sim¬ 
plified energy function (Eq. |3 2ll[ ) against the full energy function (Eq. 3.6( |, the 
model predicted beam shapes for each of the end point positions and orientations 


(specified in Table 3.11 showed no change based upon which energy function was 
used (see Fig. 


). This suggests that the inclusion of t he strain energy from 
direct shear is negligible for t his application. Further, making this simplification 
reduced the average computational time for these six test cases from 1.71 seconds 
to 0.15 seconds. 


A comparison of the predicted end loads and stored potential energy between the 
IFM using the full potential energy function (Eq. |3.6| and the simplified energy 
function (Eq. |3.20fr is shown in Table pT2| The results show some variability 
between each predicted end load. This is expected due to the differences in the 
formulations. However, the discrepancy is often not huge ami thus justifies the 
simplification of optimizing only bending energy in the calculation in order to 
realize the gains in computational efficiency. 

















09 



Figure 3.10: A set: of six plots showing the output of the IFM usiug the simpli¬ 
fied energy function (Eq. 3.21 )[ | compared against the IFM using the full energy 
function (Eq. 
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The beam shape predicted by the model using a sixth-order polynomial parame¬ 
terization was t hen compared against the experiment idly recorded beam shapes. 
Figure [iPT] shows comparisons between predicted and experimentally determined 
beam shapes for six unique sets of input position and orientation. For each set, 
the solid line represents the beam shape predicted by the model and the dots 
represent the corresponding experiment idly determined beam shape. 


A set of predicted beam shape solutions was also generated using the IFM with 
an eighth-order polynomial parameterization. These results were then compared 
against the experimental beam shapes and the sixth-order polynomial parameteri¬ 
zation as shown in Fig. T121 For this comparison, the same end conditions found 
in Table [T 1| were used. For each set, the soliil line represents the beam shape 
predicted by the sixth-order model, the dashed line represents the beam shape 
predicted by the eighth-order model, and the dots represent the corresponding 
experimentally determined beam shape. For each of the evaluated end coordi¬ 
nates, the eighth-order polynomial model produced nearly identical predictions to 
the sixth-order polynomial model. 


As a point of comparison, the use of a linear spline parameterization for the 
beam shape was also explored and compared against the experimental results. 
These model predictions are shown in Fig. |3.18| The experimental results are 
again shown by dots while the spline models are shown os solid, dashed, dash- 
dotted. and dotted lines for 6-segment. 8-segment. 10 -segment, and 28-segment 
splines, respectively. As can be seen, the lower segment number spline models had 
difficulty producing a reasonable beam shape for some position and orientation 
sets, most notably Sub-Figs, (b) and (f). 

A comparison between the sixt h-order polynomial parameterization and t he 28- 
seginent spline parameterization relative to the experimental beam shape is shown 
in Fig. |&l4| Mere, the polynomial solution is depicted with a solid line while the 
spline solution is shown with a dashed line. 
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Figure 3.11: A set of six plots showing the output of the IFM with 6th order poly¬ 
nomial parameterization compared against the experiment ally determined beam 
shape where [x r „^, y r:ni i.O eni i] correspond to the sub-figure references in Table 


3.1 








































72 



Figure 3.12: A set of six plots showing the output of the IFM using a 6th order 
polynomial, an 8th order polynomial, and the experimentally determined beam 
shape where [x rlu {, y fIl j. correspond to the sub-figure references in Table 


3.1 
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<C) 


(<l) 



<C) 


(0 


Figure 3.13: A set of six plots showing the output of the IFM using spline parain- 
eteriziit ions. 6-segmeut» 8-segment, 10-segment, iuul 28-segment, [j’,,,,/. y, n <i- 0emf] 
correspond to the sub-figure references in Table I 
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Figure 3.14: A set: of six plots comparing the output of the IFM using titli or¬ 
der polynomial parameterization versus the 28-segment spline parameterization. 
[r r „rt, i/rnrf-^rrirfi correspond to the sul>-hgure references in Table 
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A comparison of the spline predicted end loads with both the polynomial pre¬ 


dicted and experimentally determined end loads Is also shown in Table |3.3| The 
load applied to the distal end of the beam was recorded using the six-axis load 
cell. As this experiment was conducted using the planar (two-dimensional) case, 
only F : . F r . and M u loads are shown. The load cell was oriented such that the 
y-eomponent of the force and the x- and /.-components of the torque were zero. 
Thus only three of the possible six measurements from the load cell were used. 
These results are represented graphically in Fig. |3.15| The results show quite a 
bit of variability between the parameterization methods. The error between the 
experimental and model predicted end loads for both the sixth order polynomial 
and 28-segment spline parameterizations are plotted against the strain energy cal¬ 
culated from the polynomial parameterization for each of the tested experimental 
conditions. The strain energy used here is the optimized strain energy as calcu¬ 
lated using the method demonstrated in Section 13.3.1 


There is some observed 


dependence between the absolute force calculation errors and the calculated strain 
energy. However, on a percent error basis. this trend is not as dear indicating 
that additional sources of erroi such as experimental measurement or calibration 
error are likely to be significant. Further, qualitatively it was observed that lit¬ 
tle variation was present between parameterization order with the exception of 
the lower segment number spline methods which occasionally produced largely 
different force and moment predictions. 


To provide a means of comparing the relat ive quality of the model predictions, it: 
was necessary to develop a metric for this measure based upon the coefficient of 
determination. The coefficient of determination itself is a measure of the goodness 
of fit between a dataset and a prediction of those values [12-1| . Typically, the 
coefficient of determination (/?-) is calculated using the ratio of the residual error 
between t he dataset and the model (unexplained variation) to the tot ill variation of 
the dataset using one dependent variable and one independent variable. Similarly, 
the coefficient of multiple determination calculates the same value but now with 
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Tabic 3.3: Experimental versus model predicted end loads for a sixt h-order poly¬ 
nomial model and a 28-segment spline model. Forces are given in N and moments 
in N • xn. 


Sub- 

Fig* 

Polynomial Predicted 

F : F r M y 

Spline Predicted 

F : F r M„ 

Experimental 

F ; F, M y 

(a) 

0.13 

-0.11 

0.0039 

0.17 

-0.04 

0.0031 

0.10 

-0.02 

0.007 

(b) 

0.02 

0.07 

-0.0020 

0.04 

0.09 

-0.0017 

0.05 

0.06 

-0.05 

(c) 

0.57 

0.38 

0.0054 

0.50 

0.32 

0.0040 

0.28 

0.15 

0.20 

(d) 

0.01 

-0.02 

-0.0007 

0.05 

-0.02 

-0.0015 

0.04 

-0.01 

0.00 

(e) 

0.27 

-0.03 

0.0045 

0.21 

-0.04 

0.0006 

0.15 

0.00 

0.06 

(0 

0.18 

-0.24 

0.0066 

0.13 

-0.12 

0.0047 
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Figure 3.15: Error between experimental and model predicted end loads for both 
6th order polynomial and 28-segment spline parometerizations plotted against 
strain energy calculated from the 6th order polynomial parameterization for each 
of the tested experimental conditions. 
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two independent variables and one dependent variable. However, as a means of 
calculating the goodness of fit between the IFM prediction for the beam shape 
and the experimental data, it was desired to parameterize the calculation on the 
basis of the arc position, a. As such, the standard /?-’ calculation was modified to 
account for variation in both coordinate directions. 


Under these conditions, the H~ value was calculated as 


/?-' = !- 


ss r 

SS, 


(ill 


where SS ret is the residual sum of squares (E<j. 
of squares (Eq. |3.44ft . 

SS r „ = Y, 

i 

ss* = Y. Pf 


(3-42) 


and SS/at is t he total sum 


(3-43) 

(3.44) 


Here, t, represents the residuals (Euclidean distance) between a data point and 
the corresponding point on the beam curve predicted by the model as shown in 
Fig. 3.1(> The value for e, was found by first determining the position, s. along 
the beam matched to each data point. Then the location of the corresponding 
position, s. on the model predicted beam was found. The distance between these? 
points was computed and used to represent the residual for each data point. 
The value P, represents the Euclidean distance between each data point and the 
average of the location for all the points within the dataset, ;dso shown in Fig. 

cn<3 

Using this calculation for the R- value, the sixth order polynomial and the 28- 
segment spline parauietcrizatious were compared. The results are shown in Table 
|3.4[ along with values for the total sum of the residuals for all of the data points 
as well as the maximum calculated residual. The output of the model for each 
parameterization is plot t ed in Fig. |3.141 for each of the corresponding experimental 
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Figure 3.16: Graphic illustrating the variables used in the calculation of the coef¬ 
ficient of determination. 


conditions, it can he seen that the model achieved R~ values greater than 0.98 
with either parameterization for each of the experimental conditions. 


le 3.4: Coefficient of determination, total residual error, and maximum resid- 
error calculated for each sub-figure. Values presented for both the sixth or- 
polvnomial and 28-segment spline paraiueterizatious. Parameterization with 
higher /?*' value for each dataset is indicated by boldface type. 


Suls 

Fig- 

6th Order Polynomial 

R~ tint f -mas 

28-Segment Spline 

R~ e lal fmar 

(a) 

0.99(539 

(5.85 

1.23 

0.99741 

5.33 

1.15 

(b) 

0.99799 

5.07 

1.10 

0.99676 

6.97 

1.19 

(c) 

0.99865 

3.90 

0.(53 

0.99525 

7.51 

1.26 

(d) 

0.99886 

3.(52 

0.93 

0.99642 

6.98 

1.28 

(e) 

0.99533 

8.10 

1.37 

0.99759 

5.15 

1.19 

(0 

0.99504 

7.76 

1.64 

0.98401 

— 

14.65 

3.03 

— 
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While there is no clear consensus on what constitutes an acceptable R ! value, it 
Ls possible to use it as a means of comparision between the different experimental 
conditions and the parametcrizations. From Table '1.4[ it is oiiserved that each 
parameterization achieved a higher It- value than the other on some of the ex¬ 
perimental conditions as indicated by the boldface type. This suggests that each 
parameterization is roughly equivalent in terms of prediction accuracy. There¬ 
fore. there Ls no clear benefit in terms of predicting the beam shape to using 
the 28-segment spline over the sixth order polynomial since the former requires 
significantly more elfort to compute than the latter. 


3.5 Summary of Manipulator Backbone Model 

The alternative beam model presented in this section has shown potential as a 
relatively simple method for predicting the shape of a continuum beam manipu¬ 
lator when driven to a prescribed tip location and orientation. The model uses 
inverse kinematics to calculate the potential energy stored within the beam and 
then selects the beam shape with the lowest potential energy ns the most likely 
solution among the set of possible beam shape solutions. These calculations are 
performed under the assumption that all externally applied loads occur at the 
distal end of the manipulator and that strain due to bending and torsion of the 
beam are the primary contributors to the total distortion energy stored within 
the beam. Additionally, the model neglects the effects of the beam mass itself 
and tints care should he taken if the length to diameter ratio of the manipula¬ 
tor backbone is particularly large. As the backbone diameter is decreased, the 
manipulator length would also need to be reduced in order to prevent large tip 
deflections due to the beam mass. 

The results shown in the previous section demonstrate that the IFM does a rea¬ 
sonable job at predicting beam shape for all of the beam tip configurations tested 
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with a maximum observed residual error of 3.03 mm in the most extreme ease. 
The model exhibited some difficulty in predicting beam shape for some of the 
tested conditions. However, us shown in Fig. 3.10| these errors did not change 
with or without the inclusion of axial and transverse strains. Thus it is likely 
that these errors can be attributed to experimental measurement error and errors 
in the estimation of the material properties. The performance of the IFM when 
using an eight degree polynomial was also evaluated as shown in Fig. 3.12 Here 
it was found that the eight degree parameterization did not perform significantly 
better than the six degree parameterization when it came to predicting beam 
shape. This would suggest that there is little benefit to using the higher order 
parameterization since it comes at the cost of computational effort and produces 
no consistent gains in accuracy. 


Next, a spline parameterization was explored as an alternative to the polynomial 
solution method. The overall observation, as expected, is that the accuracy of 
the spline model when predicted the experimental beam shape improves as the 
number of segments increases as shown in Fig. |3.13[ For each of the tested end 
locations corresponding to Sub-Figs, (a)-(f). the spline model was found to con¬ 
verge with 14. 18. 19, (i. 11, and 20 segments, respectively. Further, with the 
exception of the positions in Figs. |3.14e and 3.14l| the spline model converged to 
the beam shape predicted bv the polynomial model. It is idso noted that except 
in the higher curvature positions, the 0-segment spline model did a reasonable 
job at predicting the beam shape. The 28-segment spline model produced beam 
shape predictions that are better than the polynomial model for some positions 
and w r orse for other positions, its shown in Table 3.4| No clear trend indicating 
which model better predicted the solution was apparent from these results. Over¬ 
all, the beam shape predictions produced by the IFM demonstrated reasonable 
accuracy for the majority of tested boundary conditions by either of the tested 
parameterization methods. While the spline parameterization method may be ca¬ 
pable of producing better results under some circumstances, this method typically 
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required significantly more lice variables to produce a viable solution. Possible 
sources of error in this investigation include measurement error for the experimen¬ 
tally determined beam position and orientation as well as errors stemming from 
the tolerance built into the optimization algorithm (set to 1mm and 0.057°, for 
position and orientation respectively). These tolerance settings were appropriate 
as the IFM predictions were found to be relatively insensitive to small changes 
in the position and orientation boundary conditions; meaning that changing the 
boundary conditions did not drastically change the shape of the resulting beam 
prediction. 

The data presented in Table |3.3| shows the predicted end loads versus the ex¬ 
perimentally measured end loads. Predicted end load data was calculated using 
both the sixth-order polynomial parameterization and the 28-segment spline pa¬ 
rameterization. It was found that the IFM predicted results for the beam forces 
at the tip for a given beam shape did not correlate well with the experimental 
measurements when using either met hod of parameterization. The spline model 
produced slightly better predictions for some end loads, but this result was not 
observed consistently. Better end load predictions from the spline model are ex¬ 
pected since the spline model uses more variables to produce a solution and thus 
has more flexibility. It should also he noted that limitations of the experimental 
measurement hardware likely contributed to errors in t he force and torque mea¬ 
surements since the loads being measured were near the lower end of the limits 
for accurate measurement for the six-axis load cell (estimated resolution of 0.063 
N and 0.0008 .V - m for force and moment measurements, respectively). Future 
efforts should work towards producing a more robust validation of the end load 
predictions using measurement equipment with greater precision. 


Chapter 4 


Artificial Muscle Actuation Unit 


A robot for any purpose requires some method of introducing power into the 
system, often in the form of an actuator (e.g. motor driven, piston and cylinder, 
or artificial muscle based). Thus, this chapter first provides an overview of the 
available actuation options and rationale for the selection of the artificial muscle 
actuator in Section |P| Section |Q| provides the mathematical description of the 
muscle actuator and is followed by a discussion of the specification process for 
the actuator parameters and the analysis that contributed to it in Section 
Experimental results are presented in Section |4.4| to validate the appropriateness 
of the models and justify t he selection of the actuator. 


TT 


4.1 Motivation and Related Work 

There are a variety of methods available to achieve the controllable linear actua¬ 
tion necessary' to articulate the segments that allow for positioning of the distal 
end of the device. Some of the possibilities include electromechanical linear actu¬ 
ators. pneumatic actuators, or hydraulic actuators. A comparison of these three 
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choices reveals that each option has both advantages and disadvantages. Elec¬ 
tromechanical actuators are typically considered to be more efficient than their 
pneumatic and hydraulic counterparts, while hydraulic actuators have the highest 
power to weight ratio of the three choices presented here, assuming that the hy¬ 
draulic power source Ls located remotely |75) . By prioritizing the design coils! mints 
of the present device, it was determined that the size restriction due to the nature 
of natural orifice surgery was the most critical condit ion to be met. followed by 
the force output from the device that is necessary to achieve desired performance. 
Therefore, hydraulic actuation was chosen for its higher power output despite any 
possible lack in efficiency, which is less critical in surgical applications. 

Although hydraulic actuation offers high power density, mechanical rigidity and 
high dynamic response, the force capability, F. for a given pressure, P. of a 
conventional linear hydraulic actuator Ls limited by the piston area, .4, given 
by: F = P - A. In view of the stringent diametric requirement of the present 
application, the ilsc of hydraulic artificial muscle actuators (AMA) (see Fig. |4.1[ i, 
which will be shown to have better peak force capability than conventional linear 
nctuators for the some pressure and diameter, will he explored |I25| . While the 
concept of artificial muscle actuators and their use in robotics Ls not new. the 
majority of applications use pneumatic power to drive them p>71 108111201 1127) . 
However, through the iuse of a hydraulic medium it may he passible to mitigate 
some of the issues with responsiveness and rigidity that have been encountered 
previously. There are two important advantages to a hydraulic approach: 1) a 
higher pressure (by n factor of 10 compared with pneumatic) can be applied so 
that: force and power density can be further increased and the actuator diameter 
can be decreased: 2) liquid has a much lower compressibility and therefore better 
rigidity than compressed gas. 

Artificial muscle actuators consist of a contained internal bladder, surrounded 
by a flexible, braided outer sheath. It Ls the geometry of this outer sheath that 
transmits radial expansion of the internal bladder due to an applied pressure to 
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contractile force along the longitudinal axis of the muscle actuator 1 I2S| . As the 
radius of the bladder, and thus the outer mesh, increases, the individual strands of 
the inesh which are woven in an over-under crossing pattern rotate relative to each 
other and to the long axis of the actuator and shorten the longitudinal distance 
from one end of the strand to the other. The load capacity of the artificial muscle 
actuator is then a function of the geometry and orientation of the outer sheath 
and the pressure applied to the internal bladder )12!)| . As discussed in Section 
is the contraction of the AMA which provides input positioning load to the 
manipulator through t he application of tension to the attached tendon. 



4.2 Analytical Description of Artificial Muscle 
Actuators 

There are two methods present oil in the literature for modeling the transmission of 
internal pressure to contractile force of an AMA. The first is a theoretical approach 
based upon energy conservation [1‘25| . while the second is on examination of the 
force prolile of the surface pressure |13fl| . The first approach is based on the 
principle t hat energy supplied to the actuator by t he pressurized fluid must leave 
the actuator through the application of a load over some distance. The second 
approach is based on an examination of the distortion of the internal bladder 
under isobaric conditions. However, ultimately each of these methods arrives at 
the base model [127) . A summary of the first method is presented here. It should 
be noted that this model does not account for possible effects of compressibility 
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in the bladder, interactions between the bladder and the braided outer sheath, or 
any contribution from the ends of the muscle actuator. 

As shown in Fig. |4.2| adapted from C'hou and Hannaford |12t)| . it is first neees- 
snxv to dehue t he geometric parameters of the braided sheath. These parameters 
include the length of the individual braided strands, b. the number of turns each 
strand makes over the lengt h of the actuator, n, and the angle between the strands 
and the longitudinal axis of the actuator, 7 (f), which is assumed to be uniform 
for all strands within the braided sheath. 




Figure 4.2: Artificial muscle act uator geometry (left) and magnified view of braid 
geometry relationship for a single unwound strand (right). 


The overall length of the actuator. /-„(/), and the actuator diameter, D\t), can 
then be represented in terms of the constants, ri and b, and as functions of the 
braid angle, 7 (f). ns seen in Eq. ( )l.l) and Eq. l JX 2 | i. 

L n (0=6*cos7(<) (4.1) 

b - sin 7(0 

717T 



<4.2) 
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Then, calculating the volume of a cylinder and substituting in the functions for 
L a (t) and D{t), 

V(t) - jWtfUt) = -^-si«.M/)cas 7 (0- (4.3) 

The first derivatives of L a [t) and l’</) with respect to ~y(t) are calculated as 


<IU0 
•hi 0 


—6 • sin 7 (f) 


(4-4) 


dV(t) //sin 7 (f) 


(cos 2 7(f)--shr7(0) 


(4-5) 


d7(0 2rm 2 M "' 2 

From Eq. ( fT~i| ) and Eq. <JTToJ. the first derivative of V(/) with respect to L„(t) is 
given as 


dV{t) dV(t)/drf(t) -IP 


d.L„(t) dL a (t)/d^(t) 47m 2 
From the principle of virtual work we have: 


*(3 cos 2 7(0-1)- 


(4.6) 


F a (t)-‘IUt) = P a (t)-dV(t) 


(4.7) 


and solving Eq. ( |4.7| for the force output with Eq. ( )l.l| andresults in 

b 2 


hM) - PM) 


4-n- 

Jr D? 


.>M0 2 , 

I 


= P a (l) — «-[3cos 2 h(/))-l] 


(4-8) 

(4.9) 


where F„(/) is the contractile force. /’„(/) is the pressure differential across the 


bladder wall, and D a (obtained by evaluating Eq. 4.2 at q - 90°)is the t he¬ 


oretical maximum muscle diameter. Note that the strand number does not factor 
into the calculation of the contraction force. This is because this relationship is 
derived purely from the geometry of n single strand during n muscle contraction 


mid thus does not change with strand number. 
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Examining the term within the brackets in Eq. i ]4.9[ l anti observing that it can 
be greater than 1 (lim-.( f )_^o = ->• the force capability of an artificial muscle can 
be greater than that of a hydraulic piston actuator for the same diameter. Do, 
and pressure. However, this comparison is based on a maximum AMA diameter 


occurring at 7 = 90°. If F a (t) is set to zero, we see from Eq. 4.9 that ->(/) = 
54.7°. This represents the maximum achievable contraction of the muscle meaning 
that larger braid iuigles are not achievable. Thus this maximum value for 7 is 
independent of all other parameters and holds true under all conditions. This 
represents the true point of maximum contraction and corresponding maximum 
AMA diameter. D, IIIL[ . . Substituting this diameter into Eq. 4.9| it is found that 
the AMA force can be up to three times that of a similar hydraulic piston actuator 
as shown in Eq. |4.i0[ 


F a (t) = P a (t) 


-D; nar 3 cas J (',(/)) - 1 


4 ’ siir’(54.7° 


(4.10) 


It is assumed here that the force from the hydraulic piston actuator is F = P - 
.-1 where the diameter is equivalent to the maximum diameter of the AMA. In 
reality the available area may be reduced by the piston rod cross-sectional area. 
Notice, however, that t his advantage comes at the cost of having t he force/pressure 
relationship van' with 7(f) 01 actuator length f. a (f), illustrated in Fig. [4.3[ At 
the moment when the braid angle is 41.8°. the two actuators have equivalent load 
capacities. 


Refinements to this basic model have been presented in the literature which in¬ 
clude considerations for bladder end cap geometry, non-linear bladder material 
response, and frictional effects |l2(il 11271 1129) . However, the basic model has 
been experimentally shown to provide sufficient accuracy to gain an understand¬ 
ing of the theoretical system performance without introducing unnecessary model 
complexity |131| . 














88 


>* 

.*=: 

u 

10 

Q. 

ID 

□ 

o 

D 

3 

% 

33 

ID 



0 10 20 20 40 W GO 70 80 !>D 


Braid Angle (degrees) 


Figure 4.3: Ratio of force capacity for artificial muscle actuator to force capacity 
of a piston and cylinder when the maximum cross-sectional area of each actuator 
is equivalent {F = P- A nius ) as a function of braid angle. For the artificial muscle 
actuator, the maximum diameter is the diameter of the muscle at *>•(/) = 54.7°. 


The application of the AMA for use in driving a manipulator often involves arrang¬ 
ing multiple actuators in an antagonistic fashion |12S |. The result is that when 
one actuator is activated, the corresponding antagonistic actuator must stretch 
beyond its nominal length to accommodate. Under these conditions t he antag¬ 
onistic actuator not only produces a contraction load in response to its internal 
pressure according to Eq. |4.8| lmt also a force in response to the elastic stretching 
of the inner bladder. In order to account for this, the elastic stretch of the internal 
bladder is modeled as a linear spring force such that F, = k- A L a , where k is the 
spring constant and A L<, is the change in length of the internal bladder beyond 
its undeformed state. Thus. Eq. |4.i)| is modified to account for elastic stretch of 
the inner bladder when the actuator is elongated beyond its nominal length as 
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shown in Eq. |4.ilj 


F (/J , * [3 cos 2 ( 7 (f)) - 1] if A L a < 0 

' P„(t)^ • [3cos 2 (7(*)) - 1] + A- ■ A L a if A L a > 0 


(4.11) 


4.3 Braid Parameter Selection 


When designing an artificial muscle actuator for a given application, it is likely 
that the available supply pressure is known as well as the desired length, L m „, 
and maximum diameter. D maI , of the actuator which occurs when the braid angle, 
7 = 54.7°. Therefore, the optimal design is the one that maximizes the perfor¬ 
mance of the actuator in terms of the range of the muscle contraction. Due to 
the geometric constraints of the AM A as discussed in Section 4.2[ the contraction 
range is directly tied to another important geometric property, the range of pos¬ 
sible braid angles throughout the stroke of the AMA. The minimum braid angle 
(full muscle extension) is dependent on the strand diameter, /?,. and the number 
of strands in the braid, V, and occurs when adjacent strands make contact |131| . 
This relationship is presented as 


1 . 

= -sm 


-1 


D,N 

ttD 0 


= — Sill 

2 


-1 


L\N- 


(4.12) 


Recalling Fig. 4.3| which shows a plot of the relat ive force capacity of the artificial 
muscle actuator versus the braid angle of the mesh, it is possible to see that the 
force goes to zero when the braid angle reaches 54.7° as previously mentioned. 
The upper limit of the force, which is dictated by the minimum achievable braid 
angle, is shown by the left vertical line. This line can shift left or right based upon 
the characteristics of the mesh, U, and N. 


Therefore, it is necessary to consider all possible combinations of strand diameter, 
D„, and number. N , in order to optimally design the AMA. This can be done by 










minimizing the bniid angle. 7 min . Additionally, the stress within each individual 
strand of the braided mesh must also not exceed the limits of the mesh material 


and thus acts as a constraint on the determination of The tensile stress 
within an individual strand within the braid is presented in the literature )131[ as 


_ Pa(t)D(t)L(t) _ 2PD% sin 2 ~f(t) cos 7(1 > 
U 2«A*.4„.„ " -V / )■ 


(4-13) 


which is formulated by analyzing the fractional component of the hoop stress 
realized wit hin an individual strand of the braided mesh. Each of the N strands 
encircles the bladder ;< times. This result can then be compared against the 
tensile limit of the mesh material, which is typically either a nylon polymer (ay = 
, r >0 MPa) or stainless steel (a y = 500 MPa). Looking at Eqs. |4.12| and |j.l3[ it 
is seen that the effect of increasing the strand number, A", is that the minimum 
braid angle, y mirl increases (smaller contraction range) and the stress within each 
strand is reduced. 


4.3.1 Finite Element Analysis of Internal Bladder 

The use of the artificial muscle actuator as a means of power input to the sys¬ 
tem carries with it several advantages previously discussed. However, it is also 
necessary to understand the failure limits for each component. A 11 expression 
for determining stress in the braided mesh was discussed in the previous section; 
however, it is necessary to understand the failure limits of the internal bladder 
due to applied pressure as well as the bladder is likely to he the weaker of the 
two components that make up the actuator. The likely mode of failure for the 
bladder is by the wall expanding through the gaps between strands, as shown 
in Fig. |4.-1[ of the braided mesh and rupturing JUS)- Further, it Is desirable to 
avoid exceeding the elastic limit of the bladder material, which varies based on 
material composition and strain rate, estimated as 3 MPa for general latex rubber 
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1132) . Towards this end. a finite eleinenf analysis of stress within the bladder was 
performed for a variety of possible combinations of actuator design parameters. 
The relevant parameters of each diamond shaped unit segment of the bladder are 
then the edge length. EL, and the braid angle, 7 , as shown in Fig. |4.5[ The 
edge length is a constant, measure of the mesh density and is determined in the 
manufacturing by setting the distance between subsequent strands. The value for 
the edge length can be calculated its 


el = Jl = ^i 

tiN N 


(4-14) 


where b . n, and A' are the strand length, number of turns per strand, and strand 
number, respectively, as before. The braid angle is a measure of AMA contraction 
and tints changes as the AMA is pressurized. This relationship is described in 
Section |4.2| In addition to the parameters EL and 7 . it is also necessary to 
consider the wall thickness. /, of the bladder material when determining a proper 
design to prevent failure. Further, the pressure within the bladder. P a is an 
important consideration for evaluating the stress within the bladder wall. All of 
these possible variations make determination of the failure criteria for the internal 
bladder of the artificial muscle a difficult problem to evaluate. The complexity 
of this problem is further increased due to the intricacies of shell mechanics for 
an irregular shape. To simplify the analysis, the bladder was first divided into 
the smallest repeated unit, which takes the form of a diamond where the sides are 
defined by the strands of the braided inesh. as shown in Fig. |4.5| The model used 
here assumes that the edges of the bladder segment are fixed in place by the braid 
strands. This assumption is reasonable as each side of a given segment is shared 
by adjacent segments which are each subject to the same conditions. Further, 
this model does not account for any stress relief that may be present due to the 
diameter of the strands restraining the bladder. This results in the overestimation 
of the stress within the bladder and thus means that the selection of parameters 
based on this model will be conservative. 
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Figure 4.4: Illustration of the bladder mesh bulging through a gap in the strands 
of t he braided sheat h. 

A model of this individual unit was developed in the finite element analysis soft¬ 
ware package, Abaqus Unified FEA (Dassault Systemes). As it is possible that 
each of these four parameters will vary independently, it is necessary to evaluate 
all possible combinations. To accomplish this, the input file was scripted using 
the Python progmmmiug language to allow for automation of the job submission 
process. This script is shown in Appendix[C] Briefly, the script was set up to 
evaluate for six possible edge lengths (€ 0.25 mm. 1.5 nun]), six braid angles 
(e [5°, 55°|). five bladder thicknesses (e [0.05 mm. 0.25 nun]), and five internal 
pressures (e [0.25 MPa. 1.5 MPa]). This results in a total of 900 possible combi¬ 
nations of EL. 7 , t. and P„. The range for each parameter was selected to cover 
the likely range encountered bv the AM A when used for this application. For each 
combination, the script defines the geometry of the bladder segment, applies the 
boundary conditions and pressure load, and submits the job to the solver. 



Figure 4.5: Graphic showing relevant planar dimensions describing an individual 
bladder segment. 
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Am example of the results output by the Abaqus simulation is shown in Figs. 


0.15 mm, 700 kPaj. Figure fl~G| shows the deformation of n bladder segment for 
n given set of parameters. As expected, the greatest deformation occurs at the 
center of the segment while the deformation at the edges is zero ns specified by the 
boundary conditions. Here it is shown that the maximum calculated deformation 
is on the order of 0.5 nun for this example. Figure f-T7| shows the magnitude of the 
Von Mises stress within the bladder segment for the example set of parameters. 
Here it is shown that the highest stress concentration is calculated at the edge of 
the segment near the obtuse angle. The maximum stress at this location is on the 
order of 2.8 MPa. Moving away from the boundary, the next highest stress occurs 
at the same location as the maximum deformation within the segment (geometric 
center) and is on the order of 2.0 MPa. 


4.6 mid 4.7 where the input piuameters were set to [i?Z.,7,f.P a ]=[1.0 mm. 15° 
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Figure 4.G: Example output of Ahnqu.s simulation results for deformation of an 
internal bladder section. 


The results from each of the 900 simulations (each combination of EL. 7 . t. and 
/*„) were then stored in individual output tiles. The result of interest in this 
ease is the maximum stress that occurs in the bladder segment. To extract this 
information from each of the results hies, a post processing script was also written 
in Python (see Appendix [Tj to read the Von Mises stress data for each element in 
the mesh and then find the maximum |133| . This maximum Von Mises stress was 
then written to a text lile along with the corresponding parameter values for each 
of the 900 simulations [13-lj . This text hie represents a four dimensional array 
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Figure 4.7: Example output of Abaqus simulation results for Von Mises stress of 
an internal bladder section. 


pTun* 
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that can he used as a lookup table to Hud the stress corresponding to a given set 
of input parameters. 

The array contains the values for maximum stress calculated under the given range 
of selected input parameters which thus necessitates the use of an interpolation 
to extract the relevant information for any one specific set of input conditions. 
For this, a cubic 4-dimensional interpolation function was used. As shown in 
Fig. [4.8j it was found that for any arbitrary set of input parameters where only 
the braid angle was varied the maximum stress occurred when the braid angle 
was 46° for all combinations of the inputs. P„, A", and t. Note that the strand 
number. A', was used in place of EL for presentation of these results ns it is more 
relevant ns a design parameter. For a given value of A'. EL is calculated using 
Eq. |4.11| for submission to the interpolation function. The value for L\ was set 
to 6 mm for this application. This result for the braid angle is true for any set 
of input parameters and is not unexpected since the surface area of the diamond 
shaped unit of the internal bladder is largest when the braid angle is 45°. Based 
on this result, it would be appropriate to use a braid angle of 4o 5 in subsequent 
calculations in order to account for the worst case scenario. 

As discussed in Section [PI the strand diiuneter. D,. and number, A\ arc the 
two critical parameters for desigu of the braided mesh. In terms of evaluating 
the bladder stress, strand diameter affects only the stress concentration at the 
boundaries of each bladder segment (larger strand diameter equates to lower stress 
concentration at the boundary). Thus assuming the most conservative condition, 
strand diameter is not factored into the calculation of the bladder stress. Instead, 
strand diameter is relevant in determination of failure within the braided mesh 
itself. An evaluation of the bladder stress was performed with respect to the effects 
of the various independent parameters. The braid angle was fixed at 45° (worst 
case) while the strand number, A\ actuator pressure. P a , and bladder thickness, 
t. were varied individually. This analysis Is shown in Fig. |4.9[ First, it is shown in 
Fig. |4.!)a| that there is a direct linear relationship between actuator pressure and 
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Braid Angle (degrees) 


Figure 4.8: Calculation of maximum bladder stress for variable braid angle. Values 
indicated in the legend correspond to /’, (MPa), A\ I (mm)]. 


stress in the bladder wall. Thus a change in actuator pressure yields a proportional 
change in bladder stress with a constant of proportionality that ranged from 0.05 
(A'=80. i=0.25 mm) to 31.5 (A*=10. /=0.1 mm). 

Next, to look at the effects of strand number and bladder thickness, the actuator 
pressure was set to 700 kPa while the strand number. N , and the bladder thickness, 
f, were varied. As shown in Fig. 4.9b varying the strand number. A*, shows that 
the bladder stress increases quickly as the strand number approaches 15 where the 
spacing between strands becomes large. Thus, for this application, the bladder 
stress can be used to set the lower limit for the strand number. It is desirable 
to approach this lower limit as fewer strands means that the performance of the 
actuator is better due to having a longer contraction range. 


Finally, the effects of bladder wall thickness were examined. Figure 4.9c shows 
that the wall stress increases quickly when the wall thickness is smaller than 
approximately 0.15 mm and changes more gradually above that value. This result: 
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Figure 4.9: Plot of the bladder .stress calculated using FEA versus (a) actuator 
pressure, (b) strand number, and (c) bladder wall thickness. Braid angle is set 
to 45° for all coses. Horizontal line drawn at 3A/Pa to indicate estimated elastic 
limit of a latex rubber bladder material 
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is important to consider ns the bladder wall thickness plays a role in determining 
the minimum achievable diameter of the AMA. 


4.3.2 Parameter Selection Methodology and Fabrication 


The previous text provided a menus of designing the AMA by analyzing the 
constraints in terms of strand stress and bladder stress for a given set of input 
conditions. Equation |4.13| when combined with the results presented in Section 
|4.3.1[ allows for a determination of the optimal combination of strand diameter. 
D„ and strand number, N. within the braid for which the strand stress does not 
exceed t he tensile strength of the braid material and the stress in the bladder wall 
does not exceed the limits of the bladder material. The thickness of the bladder 
Wiill is also an independent design parameter that can be minimized in order to 
allow the actuator to reach full elongation and therefore the wall thickness can be 
included as n design variable. 


When designing a given muscle actuator, it is assumed that the desired maximum 
elongated length. L„ uir , and the desired maximum diameter, D, na ,, are known. 
Prom these two input parameters, the constants representing the strand length, 
b, and number of turns per strand, n, can he calculated using the optimized value 
for % mi (Eq. |- 1 .12[ as 

6 = Lmat (4.15) 


cos-y miri 


b _ 6 sin 54.7° 
uD q ~D, naJ 


(1.16) 


For the braiding machines that manufacture this sort of braided mesh, the neces¬ 
sary input parameter is the pick count or the number of times the braid is repeated 
per unit length |135 |. This input setting can be calculated from the optimized 
strand number, A\ and the number of turns that each strand performs, n, over 












the nctuator length ns 
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Pick = 


nN 

2r 


(4.17) 


where the length. L c , Is the actuator length at the diameter of the core, D c that 
the mesh is being braided onto such that Eqs. -hi and |4.2| become 


L, = b • cos 7 C 


(4.1S) 



(4.19) 


Therefore, for a given combination of wire diameter, D„ and strand number, .V, 
the braiding machine can be configured using the appropriate pick per unit length 
setting calculated from Eq. |4.17| 


4.4 Testing of AM A Load Capacity 


An evaluation of the accuracy of the predicted beam load capacity as a function 
of AM A extension and pressure, as formulated in Section |4.2| was carried out. 
It was not possible to produce a muscle actuator at a scale appropriate for this 
application as an inner bladder material with the correct diameter was not found 
to be available. Thus a larger version of the muscle actuator, with an S.S mm 
maximum outer diameter, was produced (Fig. 4.10| l using latex surgical tubing 
(OD 3.2 nun, ID 1.6 nun) as the internal bladder and nylon expandable mesh 
(C)D 4.4 mm. ID 3.2 mm) as the outer sleeve, for which the minimum achievable 


braid angle was approximately 16°. The actuator from Fig. [-1.1 ()| was connected 
to a rigid support at one end and to a calibrated spring scale (OHAUS, 4 kg 


capacity) on the other end (Fig, 4.1 lb. The internal bladder was inflated using mi 


instrumented syringe (BARD Caliber Inflation Device) which provided a measure 
of the inflation pressure. 
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Figure 4.10: Photograph of a prototype art ifinal muscle actuator. 



Figure 4.11: Illustration of test setup for measuring actuator load as a function 
of braid angle. 


The procedure was performed under two static pressure conditions. (>89 kPn and 
413 lcF'a. The pressure within the bladder was increased to the static set point and 
held constant while the load on the actuator was increased. For each data point, 
the contraction force exerted by the actuator on the spring scale was measured as 
well as the length of the actuator and corresponding braid angle. These results 
are shown in Fig. |4.12| where tlie solid line represents the predicted contraction 
force calculated using Eq. |l.8| nnd the dots represent the experimental data. The 
results show reasonable agreement. It- values of 0.9578 for the <>89 kPn case and 
0.9079 for the 413 kPa case, between the predicted and experimental forces. Better 
agreement was seen for the high-force, elongated actuator condition and greater 
discrepancy for the low-force, contracted muscle condition with the largest error 
occurring at a braid angle of 50° where the error is 50% of the predicted value 
(G89 kPa case). These finding are similar to results reported in the literature 
11291113f| . Deviations between the experimental and theoretical actuator forces 
are likely due to frictional effects which become more observable at higher braid 







angles. 
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Figure 4.12: Plot of force output from the artificial muscle actuator predicted 
analytically (solid line) and determined experimentally (dots). 


4.5 Summary of Actuation Unit Design 

The use of artificial muscle actuators for robotic applications lues been expanding 
as they provide several advantages including dexterous mobility and compliance 
when coming into contact with other surfaces or objects |13ti) . This is particu¬ 
larly important in the application of minimally invasive surgery where the robot 
maneuvers in an unpredictable and sensitive environment. Further, the use of a 
hydraulic artificial muscle actuator for this purpose provides the opportunity for 
greater force output for the given size constraints, as shown in Fig. |4.4| where it is 
seen that the theoretical load capacity of the AM A is up to three times that of a 
conventional hydraulic actuator for the some diameter. Noting, however, that this 
advantage is lower in practice based on the real value for the minimum achievable 
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braid angle. The real force capacity advantage is likely in the range of 2.5 to 
2.85 times the traditional hydraulic actuator based on an achievable minimum 
braid angle of between 20° and ltr. In Section [PI a method for modeling the 
AM A and its load capacity was given showing that actuator force is a function 
of internal pressure and actuator length. The procedure for delining this method 
includes several simplifying assumptions such as the absence of considerations for 
frictional effects; however, it was demonstrated that at the prototype scale for 
this application, the predicted output compared well with the experimental re¬ 
sults. If the actuators were to be decreased in size by increasing the operating 
pressure, the effects of these errors would become more pronounced |128j . The 
determination of these errors requires experimental analysis of appropriate small 
scale AMAs. As was shown, the maximum contraction of the AMA is set by the 
braid geometry at a braid angle of 5-1.7°. However, the maximum elongation of 
the AMA is something that con be designed and thus allows for an optimization 
of the AMA characteristics in order to achieve the greatest stroke length while 
avoiding failure. 

The methods presented here make it possible to identify the appropriate braid 
characteristics to achieve maximum AMA performance. Maximum performance 
is obtained by minimizing the achievable braid angle. The minimum braid angle 
is dependent on the number of strands within the braided mesh and the diameter 
of those strands |131) . The limiting conditions placed on these two quantities are 
the yield stress of the strands mid the stress limit of the inner bladder. The yield 
stress of the strand can be calculated directly knowing the strand number, strand 
diameter, pressure, and braid angle of the AMA (see Eq. |4.13[ i. The stress wit hin 
the bladder requires more careful consideration as it is also parameter dependent 
but is less well defined in terms of a solution method. Therefore, to determine 
bladder stress, a finite element analysis was used. Here it was shown that under 
all tested conditions the maximum bladder stress occurred at a braid angle of 45°. 
Thus iu performing an optimization calculation, the braid angle can beset at this 
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value to eliminate one of the variables. Further, since it is desired to minimize the 
achievable braid angle, the bladder wall thickness can be set at a low value for this 
calculation as well while having the option of increasing the thickness to permit 
higher operating pressures. Finally, under the appropriate stress constraints, it 
was shown that it is possible to calculate the smallest allowable strand diameter 
and the smallest number of strands that can be used to produce a usable AMA. 

Using a prototype AMA, it was demonstrated that a load capacity of 40 N was 
possible with a bladder pressure of 089 kPa. However, this value does not repre¬ 
sent the maximum achievable value since the muscle was not in its fully extended 
condition when the limit of the spring scale was reached. The theoretical pre¬ 
diction was found to be accurate at the tested braid angles with the exception 
of observed deviations as the braid angle approached its maximum limit. If the 
theoretical calculation for AMA load capacity is extended towards smaller braid 
angles, then the load capacity of the prototype AMA would approach 80 N as the 
braid angle approached 10°. If we then extend this model to the design scale with 
a maximum outer actuator diameter, Do = 6mm, the predicted load capacity of 
the AMA would be 37.2 N as the braid angle approached 10° for the same supply 
pressure. 



Chapter 5 


Control Valve 


This chapter describes the design and testing of n novel hydraulic control valve 
for tine in the minimally invasive surgical manipulator. Section [f77T| provides moti¬ 
vation for the development of this design followed by a description of the various 
considerations relevnnt to the valve and its application, such as the method of 
manipulation for the valve. In Sections and |5,3| a theoretical valve design is 
presented as well as a mathematical description of the operating principles. Next, 
Section |5.4| discusses two possible methods of valve activation, piezoelectric and 
electromagnetic. In Section |5.5[ two physical realizations of the valve design are 
presented as well as a discussion of how they operate relative to the theoretical 
design. Finally, Section [5d>] presents the results of experimental testing performed 
on valve prototypes which me used to evaluate various design options and help 
inform the selection of the final configuration. 
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5.1 Motivation and Related Work 
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The ability to control the flow of the fluid Is essential to any hydraulic system. 
Typically this is achieved through the use ol valves which can control whether the 
How is ON or OFF and possibly the rate at which the fluid is flowing. The method 
by which flow is stopped con include rotational motion such as in a rotary valve, 
linear motion such as in spool or poppet valves, or with n flapper valve that acts 
parallel to the flow direction. Additionally, proportioual control can be achieved 
through varying the size of the valve opening or by control of an ON/OFF valve. 
Either method carries advantages and disadvantages. A proportional valve may 
provide more straight forward control methods: however, it introduces additional 
difficulties such as power consumption which can be particularly troublesome for 
small scale applications. Alternatively, an ON/OFF valve may have a simpler 
design but considerations must be inode for inertial and transitional losses. An¬ 
other possible solution is the use of digital hydraulics which enables a near analog 
output using digital or ON/OFF components |137| . For all methods, precise flow 
control is an integral component of a successful hydraulic system. 

There hiis been a significant amount of research committed to the design and 
evaluation of flow control valves. However, the area particularly relevant for the 
present application is valve miniaturization due to the limited space availability 
within the surgical tool. As previously mentioned, the held of fluid power at the 
micro-scale is referred to as microfluidics, a field with origins beginning in the 
1990s and has quickly expanded to suit a variety of potential applications [7N| . 
Much of the research in microfluidic valves makes use of any one of a short list 
of enabling technologies including electrostatic, magnetic, piezoelectric, thermal, 
chemical, or pneumatic methods [1331113911140111411114211143] , Of t hese options, 
possible valve operation methods such as piezoelectric or electrostat ic can provide 
interesting solutions to the problem of flow control and have been explored exten¬ 
sively for control of pneumatic systems [l44M14o| . A MEMS based design solution 
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was explored as a way to achieve How control for this application. However, sev¬ 
eral difficulties were encountered in terms of valve robustness and fabrication. A 
MEMS based design is presented in Appendix [A]as a possible design alternative 
for future exploration. 

As previously discussed, the availability of a valve for this application which meets 
the needs of the device in terms of controllability and size was not found to be 
available commercially. However, there are several approaches being pursued aca¬ 
demically to bridge the gap between traditional- and micro-fluidics. Two meth¬ 
ods presented in the literature include the use of magnetics for valve actuation 
[1-1(11 IH7| . Fu et: al. present a solution designed to use 3 nun diameter iron balls 
to facilitate clostue of the valve orifice |l-l6| . The position of the ball Is driven by 
a coil with a length of 10 mm and a diameter of S mm. The valve was demon¬ 
strated to operate with supply pressures of 50 kPa to 200 kPa and could operate 
in proportional mode up to 112.5 kPa before the ball was pulled into the valve 
seat. This design is promising but has too huge of a footprint and operates at 
too low of pressures for use in this application. Another solution presented by 
Pernod et nl. has a similar pressure capacity and is operated by a 1 nun long by 3 
mm diameter permanent magnet combined with a 500 winding coil to manipulate 
the valve membrane [1-I7| . This valve has a smaller footprint than the design 
presented by Fu et id. but does not vet meet the requirements for this application 
in terms of size and cont rollability. Finally. Peirs et al. investigated the use of 
both electromagnetic and piezoelectric methods for valve manipulation |s0]. In 
that study it was found that either approach provided significant challenges for 
producing t he desired performance safely and consistent ly. In part icular they on- 
countered challenges with integrating the valves into their design under similar 
size constraints to the present application. Therefore, given the lack of available 
solutions either commereiiilly or academically, the design of a new valve to lit this 
application and meet the necessary design considerations w'ns undertaken. 




5.2 Control Valve Design 
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In order to achieve precise manipulation, it is necessary to obtain a control valve 
that will provide a controllable actuator pressure while occupying as little space 
as possible within the device. A review of the available commercial valves and 
common micro-fluidic valve solutions was conducted and 110 existing valve design 
was found suitable for the preseut application due to the limitations of high supply 
pressure, proportional control, and met hod of activation. Therefore, development 
of a novel control valve with a diameter of (> mm, capable of manipulating high 
pressure flows, which could be on the order of 0.7 to 7 MPa for hydraulics, on a 
small scale was performed. 


For this application, the outlet of the valve will be connected to the AM A de¬ 
scribed in the previous chapter. To provide useful manipulation, it is the feed 
pressure, P a , to the actuator that must he regulated by the valve such that the 
contractile force produced by the actuator (Eq. |-1.8[ l, which is proportional to 
/’,. can he manipulated. Each actuator provides one directional contraction only 
and thus an arrangement of three actuators at a minimum operated different ially 
and antagonistically such that the manipulator can operate in three-dimensional 
space (discussed in greater detail in Section 0.2.2[ l. For each of the actuators, a 
simple hydraulic circuit such as is shown in Fig. |5.1| can be used. This circuit 
consists of an upstream, high-pressure supply, P,. and a fixed orifice How restric¬ 
tion. /?[. us well as a downstream variable orifice flow restriction, /?.»(/). and a 
return to atmosphere, P a(in . Between the two How restrictions is the feed point 
for the actuator, P a {t). By adjusting the downstream restriction relative to the 
fixed, upstream restriction it would then he possible to manipulate the pressure 
at the actuator feed point as a pressure divider. For this purpose, a flapper-style 
valve can he used as a variable orifice |14S{ . With appropriately sized orifices, 
the flapper would then act as the input for manipulating the downstream orifice 
within its available stroke range. 
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Figure 5.1: Hydraulic circuit concept for control of a single actuator. 

As each manipulator segment would require three actuators, the hydraulic circuit 
could be modified to the form shown in Fig. |5.21 which shows three sets of 
upstream orifices, three actuator feed points, and a single central flapper that 
is used to regulate the three actuators in a differential/imtagouistic fashion as 
required. As such, when the flapper moves nearer to one of the downstream exit: 
orifices (increasing the pressure at the associated actuator feed point), it is moving 
away from the other two exit orifices at the same time (decreasing the pressure at 
the associated actuator feed points). 

This style of valve has the advantages or low complexity and insensitivity to 
contaminants in the fluid [1-18| . Further, the manipulation of the valve can he 
intuitively mapped to the output of the manipulator since each have three inputs 
and are operated different ially/antagonistically. One of the basic principles of a 
flapper style Vidve is that there is a constant leakage flow. Thus the supply must 
be capable of accomodating the leakage. This presents little concern us the flow 
rate required to inflate the actuators is significantly smaller than the available 
flow rate of the supply, five liters per hour based on actuator fill time (assuming 
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Figure 5.2: Three actuator hydraulic circuit concept for the flapper style valve 
design. 

689 kPa supply pressure, 0.5 seconds fill time for a single actuator) for a design 
scale actuator versus 55 liters per hour based on available flow rate for endoscopic 
irrigation pumps |14KI 1141/] . Further, the capabilities of existing irrigation pumps 
do not. constrain the design as the supply pump could be modified to meet the 
necessary requirements. 
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5.3 Analytical Description of Flapper Based Valve 
Design 


The valve design is modeled as n flapper style valve mid thus the flow through 
each upstream orifice (Q^j) and each downstream valve orifice (Q 2< i)> shown in 
Fig. PTTH| is calctilated in Eqs. fTT| anti fTTZj 


QiAt) = A 0 C m 




(5-1) 


CkAt) = AfM?* 


-n, m 

p 


= *Da9i{t)C4 


-Pa At) 

p 


(5.2) 


where P, is the supply pressure and P 0i , Ls the pressure at each actuator feed point. 
Here it Ls shown that the upstream flow rate Ls consistent with the standard orifice 
equation, where To and Coo are the area and discharge coefficient of the upstream 
orifice, respectively, while the downstream flow exiting the downstream orifice is 
restricted by the flapper. In this case the curtain area (.-1/.,(/)) rather than the 
orifice area is used in the How rate calculation. The curtain area Ls 
where D\ Ls the nozzle diameter and g,(t) Ls the gap between the orifice and the 
flapper. Further, the discharge coefficient for the flapper orifice (C,//) is typically 
on the order of 0.8 • C | 1-18| . In both of the previous equations, /> represents the 
density of the fluid used. The gap. </,(/), acts os the input variable for control 
of the valve and varies based on the spatial position of the flapper. If the three- 
dimensional position of the flapper is described in Cartesian coordinates where 
the XY-plane describes planar motion of the flapper and the Z-axis describes 
out-of-plane motion, then the change in the size of the gap, </,-(*), can be found as 


A 9vit)‘ 


1 0 - sin ft 


A x(t) 

A fhU) 
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A y(t) 
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112 

assuming that the gap, y i, is aligned with the X-axis and ft is the taper angle 
of the flapper as discussed further in Section |5.5| Equation |5.3| assumes that the 
downstream orifices are equally spaced around the circumference of I he valve and 
that a positive change in y,{l) corresponds to closure of the gap. 

Additionally, the How rate to or from the actuator must satisfy the difference 
between the upstream and downstream flow rates for each actuator as in Eq. |5.4[ 

QM = Quit) - Quit) (5.4) 


Since Q a ,i{t) is equal to the rate of change of actuator volume, \^(/), under the 
incompressibility assumption, actuator flow rate, Q a .,(f), can be found with the 
first time derivative of Eq. |i.4[ which results in 

<M‘) = v,(t) = tL i# - ■ LA*) (5-5) 


where the actuator length. L aj {l), is dependent on beam geometry (will be de¬ 
scribed in Chapter[tjj|. This allows for a calculation of through a combina¬ 
tion of Eqs. [Oj [iT?] and [O] as 


run = 


pQl, + aw. - 




+ A j.iCar 


f 2P.AJC1 + 2 - P Ql, 


(2ific4 + 2.4 


-i 


(5-G) 


It Ls convenient to establish the working range of the control pressure by examin¬ 
ing the blocked-load condition, Q a A0 = l'* this steady-state condition when 
actuator volume Is constant. the upstream and downstream flow rates become 










equal nnd the equal ions simplify to 
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(5.7) 


This relationship is shown in Fig. 5.3 lllMj . This plot shows that there is a large 
region of approximately linear response between 0.4 and 1.6 for the orifice ratio. 
Further, it can be seen that the pressure response flattens out at an orifice ratio 
of approximately 3.0, meaning that this is a reasonable upper limit placed on the 
size of the flapper valve exit area relative to the upstream flow restriction. 



Figure 5.3: Plot of pressure ratio versus orifice ratio for a flapper valve. 


5.4 Valve Actuation 

Under the conditions described in the previous section, it: is necessary to consider 
the options available to provide operation of the valve. The two most promising 
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solutions are piezoelectric and electromagnetic 114-11 I14< r >| . Each of these valve 
actuation methods must he able to lit: within the constraints of the device and 
the available space within the valve, manipulate the three downstream orifices 
either with one component or three single-acting components, and proride analog 
manipulation over the range of motion specified in Sect ion Discussion of the 

piezoelectric and electromagnetic options is provided here with recommendations 
of the most viable solution for this application. 


5.4.1 Piezoelectric Actuation 

One possible option for valve actuation is to use a piezoelectric (PZ) actuator as 
the input method for the valve, which has the advantage of being MRI compatible. 
Piezoelectric act uators have characteristics of high force output and fast response 
time. While large forces are possible with PZ actuators, large displacements are 
difficult to produce, typically in the range of 0.1% strain as a maximum |l~>l)| . 
To overcome the strain limitations it may be necessary to consider methods of 
force amplification. One such mechanism has been developed with a capability 
of increasing the strain output to 1.4% using a single amplification or 20% strain 
when using a double mechanism [ 15 H . This type of strain amplification comes at 
the cost of force output and overall size of the act uator. An example comparison 
of the use of a PZ actuator for valve manipulation is shown in Table |5.1| Here it 
is assumed that the force output of a PZ stacit is 30-^— and the strain without 
amplification Is 0.1%. These numbers relay the fact that when an amplification 
layer Is added to the PZ stat ic, the force requirement increases and thus the cross- 
sectional area (W • 11) of the stack must get larger. It is evident from these 
dimensions that for the application to the previously presented valve, it would be 
necessary to give significant consideration to the size of the PZ stack. Another 
approach to this would be to use a hi-morph style PZ actuator. However, similar 
size limitations are present and must be overcome. 
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Table 5.1: Comparison ol PZ stark with and without amplification. 


Actuator Type 

Stack Dimensions (L x \\ x II) in nun 

No Amplification 

30 x 0.4 x 0.4 

Single Layer 

2.14 x 1.54 x 1.54 

Two Layer 

0.15 x 5.S x 5.8 


Some work exploring the use of PZ actuators for valve control has been presented 
in the literature where it was suggested that a PZ stack can be used effectively 
for this application but care must be taken to limit the effects of moisture |80| . 
Exposure to excessive humidity was found to limit the life of the actuator. Other 
limitations of a PZ based design that must be taken into account include the driv¬ 
ing voltage and operating frequency. Piezoelectric devices require high voltage to 
operate: typically in the range of 100 volts or greater. This would present consid¬ 
erable risk to the patient mid thus care would need to be talcen to mitigate this 
risk. Finally, piezoelectric devices typically exhibit poor performance at low oper¬ 
ating frequencies due to charge drain within the device which results in diminished 
actuator performance |1"2[ . Manipulation of the valve would have to account for 
this which may limit the usage scenarios of the valve under proportional control. 


5.4.2 Electromagnetic Actuation 

Electromagnetic actuation is another option for manipulation of the valve dis¬ 
cussed in Sect ion [572| The use of this type of actuation has the disadvantage of 
not being MRI compatible but has the advantages of being low complexity, inex¬ 
pensive, and able to operate in an aqueous environment. The basic realization of 
an electromagnetic actuator is such that a magnetic lield is generated by wound 
coils of small diameter magnet wire. An example of a coil produced using -1-1 AWG 
(American Wire Gauge) wire which has a diameter of 0.05 mm is shown in Fig. 
|5.4[ The key design parameter for the driving coil is the wire diameter which is 
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available commercially in a variety of discrete sizes. The resistance of the wire 
is then directly a function of the wire diameter and length necessary to produce 
the coil. Further, the current applied to the wire is variable such that the coil 
produces the desired electromagnetic force to manipulate the valve. The driving 
coil has to lie capable of closing the gap between the coil and the flapper against 
the force of the saline exiting the orifice. The force that the solenoid is capable 
of producing can lie calculated by first finding the magnetic flux density, B. from 
Ampere's Law 


B = 


Nlfi 


(5.8) 


where / is the coil current, .V is the number of coil turns, g is the gap between 
the solenoid and the flapper, and /i is magnetic permeability of the core material. 
Force generated can then be expressed as 


F = 


B-A 

17 


(Nl)-/iA 


(5.9) 


where .-1 is the cross sectional area of the solenoid. By rearranging Eq. |5.9| the 
current requiiement for closure of an electromagnetic relay can be shown as 


I = 


2g 2 F g j2F 
V-'p.4v N\ JiA 


(5-10) 


It is apparent that the closing force, F oc I~ A' 2 . Further, it is necessary to consider 
that the possible uiunbei of coil turns. N x -i- based on circle packing density 
[lf>j[ and / oc D~ for a given force requirement based on wire resistance being 
inversely proportional to cross-sectional area, where D„ Is the wire diameter. 
Putting these relationships together demonstrates that the closing force capacity 
of a given coil Is independent of wire diameter, Another consideration is the 
power consumed by the coil during operation. Knowing that the power carried 
by the wire is related to current and resistance as P — I 2 Ii. it is seen that 
P oc Di its well since wire resistance, if «A- This implies that for larger wire 
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diameter, more power is introduced to the system and thus more heat has to he 
dissipated from the driving coil in order to prevent a short by melting the wire 
insulation. This suggests that the smallest possible wire would be the most ideal 
for the coil in terms of heat dissipation from the system. One additional benefit 
to selecting the smallest possible diameter for the coil wire is that jis the diameter 
decreases, the achievable packing density for a given available space increases 
[154| . By optimizing the packing density, it Is possible to achieve the necessary 
force capacity from the solenoid for a smaller outer diameter. For problems of 
this sort. it. Is found that the packing density converges towards a volumetric 
efficiency of 90.09% as the wire diameter decreases [lofS) . For this application, 
greater than 90% efficiency is achieved for diameters smaller than 0.21 mm, based 
on the analysis of Feikert et ah. 1992 |l.il| . Therefore, any wire diameter smaller 
than this limit provides reasonable packing density. Given this analysis and the 
available sizes of magnet wire, wire with a diameter of 0.05 mm was selected for 
prototyping and testing purposes. 



Figure 5.4: Photograph of a 44 gauge wire coil used to test heat dissipation. 


Revisiting Eq. 
Sect ion ITH 



is possible to apply the design constraints set forth in 




(5.11) 




and (\jf — 0.8 • Caa together with the equation for the flapper gap 


US 


.4/ _ 3.75.4 0 
-D s -D N 


( 5 - 12 ) 


If the force requirement for the flapper is equal to the pressure drop across the 
valve multiplied by the nozzle area. As. then Eq. 15.101 becomes 



(5.13) 


Finally, the number of coils of wire, .V. that can fit within can be approximated 
using the optimal packing density and the available area. L r • D r j 2. where L c and 
D c are the length and effective diameter of the coil, as 

, r _ Lr.- D, -0.9009 

n nj 
T^u 



Putting it all together. 


3.75.4 0 - x Dl 1 2P„ 
'rD,\ L c -D c - 0.9009 y fi 


(5-15) 


Substituting in the appropriate values for this design and assuming 242 turns of 
wire are used to form the toil around a steel core, it is found that the necessary 
driving current is 41.0 mA. Therefore, based on the wire diameter previously 
selected, it is calculated that the power being dissipated into the system is on the 
order of 7.8 m\V where the resistance of the wire Is calculated using t he resistivity 
of copper. 


Since the energy entering the system through the coil is primarily dissipated 
through heat, it Is necessary to determine how much of that heat can be de¬ 
livered into the fluid that is flowing past the coil. To evaluate this, it is necessary 
to develop an understanding of the capacity of the fluid to remove heat from the 
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coil as it passes through the valve. The heat transfer coefficient for this flow is 
calculated as 


h = 


D 


■Xu 


(o.l(i) 


// 


where k w is the thermal conductivity of the fluid. Du is the hydraulic diameter 
of the flow path, auil Xu is the Nusselt number. As the How within the valve is 
likely to he turbulent based upon the How rate of the fluid and the diameter of 
the orifices, the Dittus-Boelter correlation [li>(i| was used to calculate Xu as 


Xu = 0.02:i/?e f "'/Y" 


(5.17) 


where Re is the Revnolils number (10,298) calculated on the bnsLs of the hydraulic 
diameter, Fr Ls the Prandtl number (13) estimated for 0.9% water salinity |157| . 
and u is equal to 0.4 since the coil temperature is greater than the fluid temper¬ 
ature (heating condition) |l.jS| . Finally, the heat flow in watts Ls calculated as 
shown in E<j. 5.IN where A is the surface area of the fluid interaction and AT is 
the temperature difference between the coil and the bulk fluid. 


Q = hAAT —XiiAAT 
Du 


(5.18) 


Substituting in the design parameters for tills valve (Section |(i.3[ i and the proper¬ 
ties of saline (A u . = 578.8^ [T591 ITfiOj ) it: is found that the available heat flow 
into the bulk fluid is on the order of 16G in\V per unit kelvin of temperature dif¬ 
ference. Thus as the temperature difference increases, the bent flow into the fluid 
increases as well. Since the total capacity for heat flow into the fluid is significantly 
greater than the maximum power entering the coil, the coil temperature should 
stabilize below the melt temperature of the wire insulation even when operating 
at maximum current. 





5.5 Description of Valve Realizations 
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Two possible realizations of this valve design were explored. The first consists of 
three flapper nozzle valves |1-I8[ controlled by one single tetrahedral valve plug 
that plays the role of the flapper for each of the orifices (see Figs. |5.5f . A second 
realization (Fig. |5.6[ ) of the flapper style valve design was also developed, which 
differs from the first design in that the supply enters the valve at the center and 
exits radially as shown in Fig. |5.A[ opposite the orientation of the first design. 
With this design the orifices are located on the outer diameter of the valve body 
mid are opened or closed using a ring which Is manipulated radially. 



Figure 5.5: Exploded view of the plug-based valve model showing the valve body 
and internal components. 
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Figure 5.6: Model of ring based flapper valve design. 

This valve design is normally open such that the majority of the flow bypasses 
the actuator and dumps to the return line. When the manipulated port (plug 
or ring) is activated in a particular in-plane direction, the orifices are activated 
differentially: the orifice(s) that the plug gets closer to become more restrictive, 
and the corresponding actuator pressures increase: the orifiee(s) that the manipu¬ 
lated part gets further away from becomes less restrictive, and the corresponding 
actuator pressures decreases. This in turn produces a directional actuation of the 
manipulator to bend in the direction where the orifices become more restrictive. 
Out of plane movement of the manipulated part pressurizes or de-pressurizes the 
actuators simultaneously. This has the effect of increasing or decreasing the ma¬ 
nipulator stiffness and also provides the ability to independently manipulate the 
pressure at the three actuator feed points. These orifices are sized based upon the 
necessary flow rate and the valve characteristics described by Eq. |5.7| ns detailed 
in Section [(CTl 


5.5.1 Plug-Based Valve Design 


The plug-based valve consists of a single supply pressure input which feeds three 
upstream (fixed) orifices (Fig. |5-5) . The central plug Is manipulated electromag¬ 
netic-ally (Fig. 5.7) by three independent control signals. The tetrahedral shaped 
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vilive plug can move both in plane and out of plane and can be manipulated to 
block each downstream (variable) orifice individually, in pairs of two, or all three 
depending on how the control signals are applied. When an orifice is blocked, 
pressure increases at the fe<xl point for the corresponding actuator. Conversely, 
when the orifice is unblocked. How passes to the return line and pressure at the 
actuator feed point decreases. Figure [5T»| shows one third of the How path when 
the variable orifice is open (top) and when the orifice is blocked (bottom). In the 
open position, the majority of the How passes through the valve following the least 
resistance path and exits to the return line. When the variable orifice is blocked, 
the path to the return line has a high How resistance and thus more of the fluid 
How feeds into the actuator. 



Figure 5.7: Eularged view of plug-based valve showing downstream orifices and 
location of a single driving coil. 


5.5.2 Ring-Based Valve Design 

A second realization of the flapper style valve desigu was also developed and 
prototyped (sec Fig. |5.6j l. For this version of the valve, the operation principles 
are the same as for the previous version. The primary difference for t his desigu is 
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Figure 5.8: Model of tlie control valve showing part of the How path for both open 
(top) and blocked (bottom) variable orifice conditions. 


that the supply enters the valve at the center jukI exits radially .as shown in Fig. 
5.9j opposite the orientation of the other valve version. The flapper nozzles are 
located on the outer diameter of the main valve body. These orifices are opened or 
closed differentially using a ring placer! over the orifices mid manipulated radially 
(Fig. |5.t>{ |. Further, both the valve body and the ring are tapered to allow for 
stiffness control of the three actuators together. With a greater taper angle, the 
travel distance along the long axis of the waive that is necessary for the ring to 
open or dose the downstream orifice is reduced. This in turn results in a ring 
that is smaller and thus easier to manipulate. However, the smaller ring with a 
higher taper angle becomes more difficult to manufacture mid is more sensitive 






12-1 



thus requiring more precise manipulation. 



Figure 5.9: Model of the ring-based control valve showing part of the flow path 
for hath blocked (top) and open (bottom) variable orifice conditions. 

One of the primary advantages of the ring-based valve design over the plug-based 
design is that this design is more easily produced using conventional machining 
techniques. For example, it is possible to produce this design using a standard 
lathe. 


5.5.3 Valve Connection for Supply and Return Lines 

The design of each valve geometry Is such that both the high-pressure supply 
and the low-pressure return lines can connect to the valve body in a direction 
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Figure 5.10: Illustration of the ring-based valve showing the taper angle. 


parallel to the longitudinal axis of the valve. This method of connection keeps the 
overall size requirements of the valve in the radial direction iis small as possible. 
One realization of this method of connection would employ a hydraulic line with 
the cross-section shown in Fig. 5.11[ where the black areas represent the tubing 
material and the inner white areas represent the multiple lumens of the tubing. 
This tubing is then connected to the base of the valve body such that the outer 
two lumens pass over the lower portion of the valve body and are sealed to the 
valve at the distal outer diameter as shown in Figs. |5.12| The inner lumen passes 
inside the inner diameter of the base of the valve body and is sealed to this inner 
wall or the barbed fitting. The hydraulic supply and return then behave as shown 
in Figs. |5.S| aiul [531 The use of this style of tubing permits the option of supplying 
the valve by running the tubing down the length of the endoscope working channel 
and connecting to an external pump located in the surgical suite. However, when 
more than one valve is placed in series within the working channel such as would 
be necessary with a multi-segment manipulator, providing space for the tubing 
from t he distal valve to bypass the proximal valve may prove challenging. 


Another option Is to utilize the existing irrigation and suction channels common 
to most current endoscopes. These channels exist to facilitate surgical procedures 
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Figure 5.11: Example cross-section of a multi-lumen tubing used to provide both 
feed and return lines for the control valve. 

by Hashing the .surgical site of debris and removing fluids from the site, respec¬ 
tively. Milking use of these existing fluid channels would reduce the number of 
lines running the length of the endoscope and eliminate the need for a secondary 
external fluid pump as existing endoscopic irrigation pumps are typically capa¬ 
ble of meeting the design requirements or could with minor modification such as 
increasing the available pressure [!-!*)) . 


127 



Figure 5.12: Illustration of valve connection to the hydraulic supply and return 
lumens of the described tubing for the (a) plug-based valve design and (b) the 
ring-based valve design. 
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5.6 Valve Prototyping and Experimental Results 

5.6.1 Testing of Driving Coils 

Several driving coils were fabricated for testing the dissipation of heat from the 
coil under load at a scale of 2.5 times the design scale. An example of one such 
coil is shown in Fig. [5.4[ The first coil tested had no material at the center 
of the coil (air core). The coil was placed near a fabricated ping made out of 
low-carbon steel. The coil was then energized and the gap between the coil and 
the ping at which the plug was pulled into contact with the coil was observed. 
It wns found that the gap for which pull-in occurred was les.s than 0.25 mm (the 
designed gap size) for coil currents up to 150 mA. At currents above 150 in A. the 
coil temperature quickly exceeded allowable limits resulting in coil failure due to 
melting of the wire insulation. Therefore, it was concluded that n coil with an air 
core was not suitable for this design. 

Next, a coil with a steel core material was tested under the same conditions. For 
the solid core coil, it was observed that the pull-in gap was approximately 1 nun 
when the coil current was 80 mA. Further, it was observed that the temperature 
of the coil remained at a reasonable level (warm but not hot to the touch) even 
when the coil current was applied for more than 10 seconds at a time. This 
result suggests a much greater likelihood that plug pull-in can be achieved at coil 
currents which allow for proper heat dissipation. It is important to note that these* 
tests were performed under conditions such that the primary mode of cooling was 
passive dissipation to the surrounding air. In practice the coil will be subject to 
forced How of saline which has a much higher heat capacity. 

Finally, testing of the heat dissipation of the coil was tested when the coil was 
submerged in room temperature water. When submerged it was found that the 
current in the coil could be increased to approximately 300 mA for several seconds 
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before 1 he coil temperature approached allowable levels. This observation confirms 
the notion that the current capacity of the coil is significantly higher when the 
surrounding fluid is water/saline rather than air. 


5.6.2 Valve Prototype Testing 


Plug-Based Valve Testing 


Both versions of the valve were also prototyped and tested. The plug-based version 
of the valve was printed using a Strntasys Dimension 1200es rapid prototyping 
machine at three times geometric scale such that the outer diameter was 18 mm. 
The orifices were scaled according to the flow characteristics (Section |5.3[ | such 
that the downstream orifice was 1.3 mm. upstream orifice was 0.7-1 nun, and the 
gap was 0.4 mm. A scaled prototype was selected due to difficulties associated 
with fabrication at the design scale (Fig. [571 


for which the outer diameter is 
0 mm, downstream orifice was 0.7 nun, upstream orifice was 0.4 mm. and gap 
was 0.25 mm. Shown in Fig. |5.14| one version of this valve prototype was made 
to be manually manipulated (Fig. 5.14a| >. To achieve this, a hole was drilled in 
the top cap of the valve to provide access to the inner plug. This hole is located 
downstream of the primary outlets at the same location as the intended return 
line such that the inclusion of this hole did not affect valve operation. A second 
prototype with a sealed cap was also produced (Fig. |5.14b| |. Included within this 
prototype was three electromagnetic coils intended to provide manipulative force 
to the inner plug. Each of these prototypes were fitted with a cylindrical ring to 
seal the outer supply channel and tubing to provide the necessary connections to 
the supply, return, and actuators. Each of these connections was rigidly sealed 
using steel reinforced epoxy (.1-13 Weld). 


Pressure to the valves was supplied using a rotary vane pump (Procon) driven 
by an electric motor (Fig. |o.ir>||. A relief valve was used to regulate the supply 
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Figure 5.13: Photograph of the design scale prototype of tlie plug-based valve 
design. 

pressure at a steady (>89 kPa. The How rate from the pump was sufficient to 
provide constant pressure to the valves. 

Manual manipulation of the valve plug demonstrated that it was possible to vary 
the pressure at the feed-points to each of t lie outlets that feed the actuators. How¬ 
ever. the increase in pressure at. the feed-point was not as large as was expected. 
Direct measurement of the pressure was not possible; however, estimations based 
on actuator inflation suggest a pressure of approximately 17(1 kPa. or 25% of the 
supply pressure, was achieved. From Fig. |5.3| it should be possible to achieve 
at least 90% of the supply pressure using the flapper valve. Possible causes for 
this may be poor geometric tolerances in the printing of the valve body or in the 
valve plug resulting in poor surface contact. In particular, the imprecision in the 
inside corners of the plug seat may have restricted the plug from closing the gap. 
Further, accessing the plug through the hole in the cap proved difficult and pro¬ 
vided a limited range of motion. Using the electromagneticallv manipulated valve 
in the some configurat ion produced similar results. A small change in feed-point 
pressure was observed but the change was not significant enough to produce any 
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Figure 5.14: Photographs showing the scaled plug-based valve prototype in both 
(a) manually manipulated and (b) elect romaguetically manipulated configura¬ 
tions. 






132 



Figure 5.15: Photograph of hydraulic supply setup. 

significant output at the actuator; an estimated pressure of 1(H) kPa or 15% of the 
supply pressure. This result may be due in part to higher flow forces at the valve 
orifice lor the three times scale version of the valve. Increasing the diameter of the 
valve orifice significantly increases the How forces at the valve plug. It is difficult 
to predict any additional reasons for the poor performance of the valve since it 
was not possible to observe the internal valve operation during testing. Future 
exploration of this valve design should include developing a means of investigating 
the internal operation of the valve to locate any flaws in the design oi fabrication 
of the prototype. 


Ring-Based Valve Testing 


The ring-based valve design was also prototyped and tested. However, to accom¬ 
modate the need for manual manipulation and to reduce manufacturing costs, the 
scale of the valve was increased. At prototype scale, the outer dimension was 20 
mm and the overall length was 32.5 mm. An upstream orifice diameter of 0.74 


mm was used to provide the desired flow rate and then E<|. |5.7| resulted in a 
downstream orifice diameter of 1.3 mm in the fully open condition. Using this 
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diameter, the necessary gap. y, between I lie orifice and the flapper under the fully 
open condition is 0.4 mm. Also, for this prototype, a taper angle of 3° was used in 
order to make manual manipulation of the valve ring possible. A prototype of this 
design is shown in Fig. 5.1ti The prototype was produced using CNC machining 
techniques (Island Tool k Die. Pine Island, MX) and connected to the hydraulic 
lines using steel reinforced epoxy (.1-13 Weld) and a hose clomp, as shown in Fig. 

Em 



Figure 5.1G: Machined prototype of the second valve design. 

The ring-based valve prototype was tested using the same setup shown in Fig. 
5.1ol Testing of this prototype involved manual manipulation of the outer control 
ring to differentially open and close the flapper nozzles. It was found that the 
prototype of the ring-based design performed significantly better than the first 
design. It was possible to inflate artificial muscle actuators connected at the 
actuator feed points to their fully contracted states and thus the pressure Increases 
at these locations were demonstrably superior than was observed for the lirst 
valve desigu, on the order of 650 kPa or 95''! of the supply pressure. Further, 
it was possible to simultaneously increase the actuator outlet pressures in pairs 
by positioning the ring such that two orifices were sealed. Finally, it was found 
that the relative stiffness of the three outlets together could be manipulated by 
positioning the ring up or down the taper. The effectiveness of this valve design is 
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also explored in Section [(T.4.31 where it Ls combined with the rest of the manipulator 
components to evaluate overall performance. 



<»>) 


Figure 5.17: Photographs showing the scaled valve prototype for the ring-based 
design with (a) outlet and supply lines attached and (b) assembled within the 
return line. 





5.7 Summary of Control Valve Design 
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After review of the available valve technologies for control of a hydraulic system, 
no options were found at the scale necessary for this application. The design pre¬ 
sented here is based on a single-jet flapper valve that facilitates control over the 
ratio of actuator pressure to supply pressure. The flapper valve design provides 
benefits of a relatively linear performance curve on the basis of outlet pressure 


versus downstream orifice size as shown in Fig. 5.3 [l-lS| . The t heoretical behav¬ 
ior of the valve is based on the principles of pressure drop across a flow restriction. 
As such, these predictions are dependent upon the quality and accuracy of the 
valve fabrication processes which can introduce unexpected errors due to poor 
surface finish or dimensional tolerances. For smaller valve orifices, these errors 
can become more pronounced. The valve design presented here consists of three 
such flapper nozzles for the control of three actuators. Two possible realizations of 
this design were discussed Each realization contained a manipulated component 
(plug or ring) that plays the role of the flapper for each of the actuators, differ¬ 
entially opening and dosing the three nozzles through planar motion. Further, 
axial motion of the plug or ring controLs stilfness of the three actuators together. 


Testing of the electromagnetic driving coils showed that the coil was able to ma¬ 
nipulate the plug, when a steel core was placed at the center of the winding, 
without exceeding the fusion temperature of the wire. Further, when tested while 
submerged in water, it was possible to increase the coil current nearly fourfold 
without shorting, greatly exceeding the rated current capacity of the wire. This 
suggests that the coil actuator provides a suitable method for plug manipulation 
in a submerged condition. 


Finally, it was found during repeated attempts at fabricating a workable proto¬ 
type of the plug-based valve design presented here that manual manipulation of 
the valve produced a muted effect on the pressure nt the actuator feed-point. 
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Findings were similar for the powered valve. The likely cause of this result is poor 
surface finish in the fabricated valve due to the rapid prototyping process used for 
fabrication. The ability to properly seal the nozzle exit is important for flapper 
valve operation. Testing of the ring-based valve design revealed that this design 
performed significantly better than the first under similar operating conditions. 
It was possible to increase the pressure at the actuator feed points individually, in 
pairs, or all at the same time. This result can likely be attributed to the fact that 
the prototype of this version could be fabricated using conventional techniques. 
This allowed for a better surface finish and thus a superior seal at the flapper noz¬ 
zle exit. These results indicate that this valve design represents a feasible solution 
for the control valve at design scale. 



Chapter 6 


System Assembly 


As the final chapter describing the design of the surgical manipulator, this c hapter 
focuses on t he assembly of t he various components that form the surgical manipu¬ 
lator. First, in Sect ion fT2] the models of each of the components are collected and 
assembled into a system model. This system model can be used to provide some 
insight into predicted lmuiipulator performance and outcomes. Next, in Section 


tifv what the manipulator would look like at design scale in such a way that it 
would be capable of meeting the design requirements. In Sect ion [7T-H the system 
is analyzed using the analytical models. Further, a system prototype is then pre¬ 
sented which has been scaled up from the intended design scale. This prototype 
is used to evaluate some of the design metrics and provide motivation for future 
miniaturization towards a design scale prototype. 


the design specifications for each of the components are provided to quan- 
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6.1 Motivation for System Level Analysis 
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As discussed in Sect ion [1.1.51 there are a number of surgical robot concepts being 
explored, each with strengths and weaknesses. As a means of producing a relevant 
comparison against these other solutions, it was necessary to develop an assembled 
model of the manipulator as a system. By coupling the system sub-components 
together as designed, it is possible to analyze the theoretical system performance 
and identify areas of possible improvement . These results can idso lie compared 
against the design criteria set forth in Chapter [3] in terms of system performance. 
Further, the system model can enable future iteration of the design process for 
each of the sub-components whenever there is a change in design requirements or 
application. As a whole, the system model permits virtual assessment of system 
performance for control system design. 


6.2 Analytical Description 

To develop a model of the assembled system, each of the components must be 
assembled together. The system components me coupled together through both 
hydraulic and geometric relationships. The control valve supplies pressure directly 
to the actuators. Each actuator mounted in its elongated configuration to allow 
for .antagonistic stretch and is coupled to a tendon which in turn is routed along 
the manipulator backbone. Backbone shape directly determines relative length of 
each actuator. In both the two-dimensional and three-dimensional cases, only a 
single manipulator segment Is modeled. Future work should build upon this and 
include the ability to model mult iple manipulator segments. 

For both the two-dimensional and three-dimensional cases, it is assumed that 
the number of disks along the actuator tendons are sufficiently large such that 
the angle ,i(i) between the tendon and each intermediate disk i, as shown in 
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rig. O is sufficiently small so that each tendon can he considered lo be normal 
to the disk. Further, it is also assumed that any frictional interaction between 
the tendon and the disks is negligible and thus tension within each tendon is 
efficiently transmitted through each intermediate disk. These two assumptions 
dictate that the interaction between the tendons and the intermediate disks have 
negligible effects on the beam shape (JJ7J. As such, the contractile force generated 
by each AM A is transmitted to the terminal disk by each associated tendon and 
the tendons present loading on the beam only via the terminal disk, in the form 
of a force directed along the line of action of the tendon and perpendicular to the 
terminal disk as shown in Fig. |t>.2[ but not via the intermediate disks. Because the 
tendons are attached to the terminal disk at an offset distance from the disk center, 
the net result is an applied moment and an applied force which are counteracted 
by the moment and force produced by the beam through deformation. 



Figure 6.1: Graphic showing angle change of tendon due to rout ing through in¬ 
termediate disks. 


6.2.1 System Model in Two Dimensions 


In the two-dimensional case (Fig. 6.3fr , the individual system component mod¬ 
els can be coupled together as shown in Fig. |6.4( The two-dimensional system 
includes one valve consisting of two active orifices and two passive orifices, two 
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Figure 6.2: Graphic showing resulting free body diagram for assumption of many 
intermediate disks. 


actuators, and a single segment beam that can exhibit planar motion. The two- 
dimensional beam model assumes planar bending and does not take into account 
any deformation due to beam mass or the effects of shearing or axial extension. 
For the overall system model, the control input is the size of the gaps for the 
valve's variable orifices. Since the two-dimensional model incorporates two an¬ 
tagonistic actuators, only two active valve orifices are included, which operate 
differentially, in addition to the two upstream fixed orifices. The first block iu 
Fig. (T3 is the valve hydraulics block. This block accepts the control input, g, and 
the time derivative of the actuator lengths, L a . From these inputs, the pressure 
at the actuator feed points, P a . are calculated using Eq. |5.6| (repeated below in 
Eq. 16.1(1, where Qaj is calculated as shown in Eq. 6.2 The actuator pressures 


then act as the input to the actuators block along with the actuator lengths. L„. 

for each of the actuators, the load from each 


Using Eq. |4.111 (repeated in Eq. 
of the actuators is calculated as F a> i and F a ;>. Recalling that when the actuators 
iirc stretched from their nominal length, the force load due to the actuator stretch 
is modeled as a spring force (F, = k-AL a ) which results from stretching the inner 
bladder, as discussed in Section |1.2| Under nominal conditions when the valve 
plug is centered, the actuator pressure Is GdAkPa as defined by the ratio of the 










downstream and upstream orifice diameters. 
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Posit) = 


pQl, + 2-4’C^P. - I -Vanft,. 






(2^ti + 2-4?,f'i) 


-1 


( 6 - 1 ) 


g-,<o = v;<o = (& 2 - 3i.,(o 2 ) • t,,<o 

f/^/)^.[3cos 2 h(/))-l] ifAL.<0 

I Poit)^ - [3cos 2 (M0) - 1 ] + A- - A L a if A L a > 0 


( 6 . 2 ) 

(6-3) 



Figure 6.3: Skeleton model showing the basic components that make up the two- 
dimensional system. 


The quantities relating to the actuator length. L„ and L„. which act as inputs to 
the actuator* and valve hydraulics blocks, respectively, originate from the beam 
mechanics block. The beam mechanics block takes as input the position and 
orientation of the distal end of the manipulator, [j. r,0], which is initialized in 
the undeformed configuration. From this information, the actuator length can be 
described as a function of the orientation of the distal end of the manipulator. This 
relationship can be found by examining a differential segment of the manipulator 
(Fig. |6.5] | |fifl| . If ds represents the arc length of the manipulator backbone 
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Figure 6.4: Simplified block diagram of system model. 


over tills segment and the corresponding arc length of the location of tendon 
attachment is defined as ds*. then the relationship between the two arc lengths 
can be found by first writing 


ds' = {p- r)dO (6.4) 

where r is the distance from the center of the backbone to the compressive side 
tendon, 0 Is the angular orientation of the backbone, and /> is the radius of cur¬ 
vature of the manipulator segment which can be defined as (dO/ds)~ l . Rewriting 
this expression to get 

rf.s* = ds — r - dO (6.5) 

and then rearranging and integrating both sides over the length of the beam. 

e 0 ,L\ results in an expression for the total length of the tendon, .s*(L), as 
a function of the beam length. .s(L), and the orientation of the beam at .s = L, 
0{L), .as shown in Eq. 

s*(I) = s(L) - rO(L) (6.6) 


6.6 


Tims, the change in length of a given tendon located at the outer edge of the 
manipulator Is equal to r0(L). Therefore, by knowing the orientation of the 
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terminal dusk, the corresponding tendon position and thus actuator length can he 
calculated. Note here that the tendon is assumed to be inextensible. The change 
in tendon length that is being calculated is really the portion of the tendon within 
the flexible section of the manipulator. Because the actuator is connected to the 
terminal disk of the backbone via the tendon, the length of the actuator, L a (t) 
can be determined by the change in the tendon length as the manipulator bends. 
The actuator length is given by: 

ut) = L«.o - r • 0(t. L) (6.7) 

where L an is the length of the actuator when the manipulator is straight (i.e. 
0(s) = 0 for all s € [0, L)), and r is the radius of the routing disk at the point 
of the tendon connection. For the other actuator the same relationship exists 
but instead the change in length, rt)(L). is added to the initial length for positive 
values of 0(t,L). 



Figure 6.5: A differential segment of the manipulator with the backbone and 
tendon shown. 


As it is assumed that an end-effector is placed at the distal end of the manipulator, 
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translational and rotational inertial dynamics allow for calculation of the location 
and orientation of the distal end of the manipulator due to the applied control. 
A force and moment balance Ls applied to the end-effector mass, treated as a free 
body ns shown in Fig. <>7t* Loads applied to the end mass in the global coordinate 
system originate from t he actuators, the reaction load from the backbone, F& mid 
A/fc. and from any external loading. F& and M /., due to tissue manipulation or 
collision, as shown in Eqs. 


ITTTT 


Y. = fft-r + Flj + /v.l_r + F' o2j 

(6.8) 

y F : = mi = F h; + F L : + F„i.. + F t2 . : 

(6.9) 

y -W = 10 = -F ol • r + F a i - r + M, t + M L 

(6.10) 


The result, of these equations are values for the cartesian components of the trans¬ 
lational acceleration, x and if, and the rotational acceleration. 0. which can be 
integrated twice to find the position and orientation of the distal end of the ma¬ 
nipulator. 



Figure 6.6: Free body diagram for mass at the distal end of the manipulator. 


Further, using the model described in Section |3.3.l| for a given position and ori¬ 
entation. the beam mechanics block outputs the bemn forces and moment acting 






145 


on the end effector at the distal end of the beam can be described by Eq. ln.2-'t[ 
copied below. 



OPE' 



OPE' 

92a, d 



OPE' 

0H, lul 


(e-ii) 


6.2.2 System Model in Three Dimensions 


For the three dimensional model of the system, the components that make up the 
assembled tool are the same as those for the two dimensional system. However, 
in order to achieve non-planar motion of the manipulator, a total of three muscle 
actuators and three tendons are used. Since the tendons are connected to the 
terminal disk in such a way that the tendon loads are applied in a direction 
parallel to the manipulator backbone, the resulting load on the terminal disk due 
to the tendon loads is realized as a moment in two Cartesian directions (X and 
Y) and one force (Z-direction) as shown below 


F ; 


1 1 1 


F a A ' 

M r 

= 

11 0.866 t -0.866 -r 


F a .2 

My 


—r 0.5 • r 0.5 • r 


Fai 


( 0 - 12 ) 


where the actuator loads, F ai , which are transmitted by the tendons, are con¬ 
nected to the terminal disk as shown in Fig. |6.71 Tints for a given manipulator 
segment, three degrees of freedom can be manipulated by the three actuators. 
The addition of more actuators would not provide control over any additional 
degree's of freedom. Therefore, due to the limitat ion on available space, the use of 
three active inputs for a given manipulator segment was selected for this design. 
These three actuator/tendon combinations are arranged equidistant around the 
circumference of the manipulator as shown in Fig. 16.7 












Neutral 

Axis 
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Figure 6.7: Illustration of the terminal disk showing the tendon connection points 
ns indicated by the black squares. 


The system model in three-dimensions is arranged much the same as the two- 
dimensional case described in Section 6.2.1 The key difference being that there 
me six variables necessary to describe the position and orientation of the distal 
end of the manipulator rather than three. The modified system block diagram is 
shown m Fig. |6.8| The first block in Fig. (TS is again the valve hydraulics block. 
This block accepts as its inputs the control inputs, y,. and the time derivative of 
the actuator lengths. L aJ . From these inputs, the pressure at each of the three 
actuator feed points. are calculated using Eq. |6.l[ The actuator pressures 
then act as the input to the actuators block along with the actuator lengths. 
Using Eq. O for each of the actuators, the load from each of the actuators is 
calculated us F 0i >, F Ui2 . and Again, when the actuators are stretched from 
their nominal length, the force load from the actuator is modeled as aspring force 
which results from stretching the inner bladder. 


The quantities relating to the actuator lengths, and L a j, which act as inputs 









Figure 6.8: Simplified block diagram of the 3D system model. 


to the actuators and valve hydraulics blocks, respectively, originate from the beam 
mechanics block. The beam mechanics block takes as input the position and 
orientation of the distal end of the manipulator, (x. //. z.o, 0 . til'], which is initialized 
in the undeformed configuration. From this information, the actuator length can 
be described as a function of the orientation of the distal end of the manipulator. 


This relationship is similar to what was shown in Eqs. |6.-1| - 16.7| However, in the 
three-dimensional case, the orientation angle necessary to compute the change iu 
actuator length is not simply 0- Instead it is necessary to determine the change in 
orientation of a vector, P,. which describes the location of the tendon termination 
point for each tendon relative to the center of the backbone as shown in Fig. 16.10 


The arc length derivative of the H(s) rotation matrix (Eq. 3.2-lfr can then be used 
to operate on the reference position vector. r.,o. for each of the tendon termination 
points and determine the change in orientation of r, in the deformed state by 
substituting in the values which describe the orientation of the distal end of the 


manipulator, o(.s), 0(s), and for a differential segment as shown in Fig |6.9 
For each differential segment, the length of the tendon can be found using the 
change in centerline position from one end of the segment to the other, J(.s). the 
components of which are shown in the integrands of Eqs |3.32| - |3.3-l| along with 
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the change in position of the tendon termination points, r,. as shown in Eq. |ti. 13[ 

^i = r(.,)’ +(0.13) 

Integrating the magnitude of this vector over the entire beam then yields the 
resulting length, It., of each of the three tendons tinder the given beam shape. 
This tendon length con then be compared against the undeformed tendon length 
such that the change in length of the attached actuator. , can be calculated as 

L a At) = L afi + AL tj (6.14) 



where L „.o is the actuator length in the reference configuration. 


Figure 6.9: A differential three-dimensional segment of the manipulator with the 
backbone and tendon shown. 


Further, from the beam mechanic* block, using the three-dimensional model de¬ 
scribed m Section for a given position and orientation, the reaction loads 
at the distal end of the beam can be described by Eqs. |3.38| and |3.39| (repeated 
below in Eqs. |6.15| and |6.1<1[ >. These results are then added in to the net load 














acting at the distal end of the manipulator. 
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Then using; the relationships defined in Eqs. |6.17| and |6.1<)| in the system model 
(Fig. iG.Sh, the three dimciudonal position and orientation of the end moss can be 


computed from the translational and rotational inertial dynamics similar to the 
two dimensional case. A force and moment balance is applied to the distal end of 
the manipulator, treated as a free body as shown in Fig. |0.10| Loads applied to 
the end mass originate from the actuators, the reaction loads from the backbone, 
F, and Mb, and from any external loading. Fl and Ml. as shown in Eqs. |6.17| nnd 


(i.!8| Force biilances are with respect to the global reference frame while moment 
balances are with respect to the body frame described by the angles, 6{L). 0(L), 
and 

~ ’ - - — (6.17) 


T = m j = Ft + F l + y P 0 ,> 


i=l 


y; M = /ii + 0x lltpj = J7 b + Ml + y (i t l x Fn,, 


(6.1S) 


where / contains the moments of inertia of t he end effector about the local prin¬ 
ciple axes. The quantity r, x F aj is the cross product between the position vector 
from t he center of rotation to the point where the tendon terminates for each of 


the three tendons attached to the end mass as shown in Fig. (i.!0| and the force 
vector applied at that point in the direction of the tendon in the deformed state 
and gives the contribution to the total moment due to each applied actuator load, 
M aA . Here it is again necessary' to use the rotation matrix. Eq. |3._-1| to operate on 
the position vector. i r ( , and t he force vector. F a i . for each of the actuators. Since 
















the orientation of both of these vectors are known in tlie nndefomied configura¬ 
tion, the rotation matrix can be used to represent these vectors in the deformed 
configuration. This in turn makes it possible to calculate the contribution to the 
net moments acting on the end mass by each of the act uators. 



Figure 6.10: Free body diagram for the three-dimensional mass at the distal end 
of the manipulator. 


6.2.3 General Solution Method 


The general procedure for evaluating the system model as described in Sections 
|6.2.1 and 6.2.2 is the same whether considering the two-dimensional or three- 
dimensional case. Using the analytical descriptions of each of t he system compo¬ 
nents described in the preceding chapters, the state of the system at each time 
step can be evaluated as suggested in pseudocode below. The only steps winch 
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do not allow for chxsed form calculation are those involving the beam mechan¬ 
ics as the formulation for computing the beam loads and shape makes use of a 
computational minimization method 11 s discussed in Sect ions |3.3.11 and |3.3.2[ 

set constants 
initialize variables 
for step in time 

set control inputs 

solve for actuator pressures 

solve for actuator forces 

solve for position, orientation accelerations 
integrate to get position, orientation velocities 
integrate to get position, orientation 
get beam loads from boam.mochanics 
get actuator lengths from beam_mechanics 
end for 


6.3 Component Specifications 

The specific design of each of the system components requires consideration not 
only of overall constraints placed on the system described in Chapter [2] but also 
considerations unique to each individual component (Chapters^ [I] and|5j| and the 
interactions between each of the components that make up the assembly (Section 
(i.2[ |. These values presented here describe the design scale characteristics for the 
minimally invasive surgical manipulator intended to meet the design requirements 
set forth in Chapter [5] While the requirements of any given surgical procedure 
may require modification to these specifications, the numbers set forth here should 
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provide general guidelines for a capable robot that can be employed for a wide 
variety of abdominal procedures using tranggnstric access. 


6.3.1 Manipulator Backbone 


For the manipulator backbone, the size of the manipulator must be constrained 
to lit within the available tool space and thus the diameter of the routing disks 
located along its axis must have a diameter no greater than 6 mm. Further, these? 
disks should be spaced such that t he manipulator can achieve maximum deflection. 
This is calculated using the ratio of disk radius to distance between disks where 
the ideal ratio occurs in the range of 0.2 to 0.5 as demonstrated by Li and Rahn. 
2002 [HU. Thus, for a required maximum disk radius of 2.5 mm at the tendon 
routing location, the spacing between t he disks should he € [5mm. 12.5mm| and 
most optimally at 6.25 inm [87) . 


Recall the two workspace specifications, an overall positional volume of 50 mm 
cubed and an angular range of ±30°, described in Section [?72| Figure [57TT] shows 
the working envelope for a two segment manipulator where each segment lwvs a 
length of 50 nun. The shaded region in each Sub-Fig. represents the achievable 
positions for an unloaded manipulator where the augular range of each segment 
is limited to 30°, 60°, 90°. 120°. and 180° for each Sul>-Fig., respectively. This 
working envelope can be revolved around the longitudinal axis to project it into 
the third dimension. Assinning a two segment manipulator with a bend range 
of ISO", a manipulator segment length of approximately 50 nun would allow the 
device to reach the majority of the likely workspace, ns shown in Fig. 


rrn 


If the number of segments were increased, the length of eaeli segment could he 
shortened. 


Finally, the selection of a backbone material with a high capacity for strain would 
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Figure 6.11: A set of five plots showing the planar working envelope of a two 
segment manipulator for a per segment bend range of (a) 30°, (b) 60°. (e) 90°, 
(d) 120°, and (e) 1S() C . 












aid in making the manipulator maneuverable without introducing permanent de¬ 
formation. One attractive option is nickel titanium alloy, nitinol, wliich exhibits 
such characteristics. Also, it is desirable for the beam to be flexible in bending 
due to an applied load yet aide to support itself in the unloaded state. For the 
scale of the manipulators discussed here, a diameter of approximately 0.5 mm is 
appropriate based on qualitative ol»servfttions. For Martensite phase? nitinol, the 
elastic modulus is on the order of 28 GPa |123| compared to a similar material 
such as stainless steel for which the elastic modulus is on the order of 196 GPa 
[Hil| . As such, t he bending stiffness of the nitinol option is much lower than the 
stainless steel option. 


G.3.2 Actuation Unit 


For the specification of the artificial muscle actuator, it is assumed that tIre the¬ 
oretical diameter limit, D, t . is 6 mm. Further, the braid material is assumed to 
be nylon <<7*- == 50 MPa) and the bladder is assumed to be made of latex rubber 
[ay =s 3A I Pa) |132) with n wall thickness of 0.15 mm in order to maximize the 
available contraction range. Equation |6.3| gives the relationship between actuator 
length, L„. actuator pressure. P„, and contractile force, F„. Thus two of these 
values must be specified to solve for the third. First, the desired condition of 
producing manipulator retroflexion will he used to design the actuator length. 


To design the length of the actuator it is necessary to know the required muscle 
contraction, for which there is an achievable limit of approximately 32.5% of the 
elongated muscle length |l2.i| . As shown in Eq. 6.7| the muscle contraction can 
be described as being equal to the value r • 9(L) in the planar ease. This means 
that for a single manipulator segment, manipulator bending of 180° would require 
a length change during actuator contraction of r ■ 0[L) = 2.5 mm • tt = 7.85mm. 
Due to the necessity for a given actuator to stretch by the same amount when 
the antagonistic actuator is activated, the total change in length bet wren the 
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extended and contracted .states must be doubled. Thus assuming a contraction 
limit of 52.5%. an actuator of length of 48.2 min would meet this requirement. 

Next., for the artificial muscle act uator, it is necessary to consider both the param¬ 
eters of the braid and the length of the muscle. Applying the braid optimization 
methodology’ described in Section |4.3| at design scale and assuming a supply pres¬ 
sure of 700 kPa suggests that the use of at least 30 strands and a minimum strand 
diameter of 0.14 mm yields the smallest braided mesh size while meeting the con¬ 
straints for strand stress and bladder stress, assuming nylon wire is used for the 
mesh and latex rubber for the bladder. This corresponds to a theoretical mini¬ 
mum braid angle of 5.2° and minimum diameter of I).(i7 mm. The real achievable 
minimum diameter is limited based on the deflated size of the inner bladder. From 
this value for 7„„„. the necessary strand length, b. can be calculated from Eq. |4.15| 
ns 48.4 mm. 


Using the specified design parameters yields a pick count of approximately 17 per 
inch (English units used to match braiding machine input). A sample mesh was 
produced using a braiding machine (Steeger HS 80/8VEA). A photograph of the 
machine braiding apparatus is shown in Fig- El 


The mesh sample, shown in Fig. 6.13| demonstrates the braiding machine output 
with the braid parameters specifications suggested previously. The sample 1ms 
a pick count of 15 which corresponds to an approximate edge length of 1 mm 
between each strand that crosses the long axis of the actuator. 


Finally, it is necessary to consider the resistance of the antagonistic actuators and 
the bending resistance of the beam in the retroflexed state. The beam resistance 
is estimated using Eq. |3.3| and the antagonistic actuator resistance is calculated 
from the nominal pressure applied to the actuator in the open valve condition 
as well as the contribution from the elastic stretch using a spring constant of 
0.07 N/rnm, estimated from mechanical properties of general latex rubber and 
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Figure 6.12: Photograph of a muscle actuator braiding process setup. 



Figure 6.13: Photograph of braid sample at 15 pick per inch. 


the geometry of the inner bladder |162[ . Putting this all together results in mi 
estimated actuator force requirement of approximately 5.68 N in order to achieve 
retroflexion of the planar manipulator. doing back to Eq. |6.3| it is possible to 
solve for the necessary actuator pressure under these conditions of approximately 
575 kPa. It should be noted that these calculations were performed in the ahsence 
of load in order to establish a baseline for manipulator positioning. Either the 
pressure would need to be increased or the actuator would need to be longer in 
order to accommodate any load at retroflexion due to the relationship between 
actuator length and force capacity shown in Fig. 
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In terms of the load requirement of 5 X placed on the manipulator, the corre¬ 
sponding force requirement from the actuators is dependent upon manipulator 
tip position and orientation. For example, if the external load is applied in a 
direction parallel to the axis of the manipulator backbone at the distal end, then 
the force requirement from the actuator to begin moving the load would be the 
same magnitude as the load. Actuator forces greater than the load would result 
in a change in manipulator tip position and orientation until a new equilibrium 
position is found. Thus a 5 N load applied parallel to the manipulator backbone 
would require an actuator pressure of at least 100 kPa to begin moving the load 
under the condition of full actuator extension. However, as the orientation of 
the manipulator changes, the pressure requirement would increase quickly due to 
contraction of the actuator and direction of the applied load with respect to the 
tip of the manipulator. 


6.3.3 Control Valve 


FYom the required actuator si/.e. it is possible to determine the necessary design 
parameters for the control valve. For the actuator size given in the previous 
section, the volume change of the actuator would be on the order of 0.64 mL. 
assuming the deflated volume is zero and the actuator inflates to approximately 
5 mm in diameter in the contracted state. With the design goal of at most 0.5 
seconds for full range of motion, a valve flow rate of 1.4 — would be desired. 
Using Eq. 57T described in Section |5.3| and assuming a supply pressure of 700 
kPa together with the necessary design average flow rate of 1.4 — results in an 
upstream (fixed) orifice diameter of 0.4 mm. To simulate the effects of reduced 
How rate with increasing actuator pressure, the pressure was ramped to 350 kPa 
(conservatively selected as an approximate tipper hound on actuator pressure dur¬ 
ing high speed position changes) wit hin a 0.5 second span and the total volume 
set; to the design value above. Next, applying the principle of Equation 5.7 with 


5.7 
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an orifice ratio of 3 as discussed in Section 15.31 then results in a downstream orifice 
diameter 0.7 mm in the fully open condition (flapper nozzle behaves as though 
the standard orifice area is used in place of .-1/). Using this diameter, the gap, g , 
between the orifice and the flapper is sized such that its nominal value falls under 
the fully open condition, here set to be 0.25 mm. 

For each viilve design it is also necessary to consider the taper angle for the flapper. 
For the plug-based design, the taper angle of the plug is dictated by the geometric 
limitations of the valve body. Based on inspection of these geometric constraints 
the taper angle was set to 15° as shown in Fig. |5.5| For the ring-based valve 
design, it is also necessary to consider the taper angle for the middle valve section 
as shown in Fig. |5.1tl( As previously mentioned, the higher the taper angle, the 
less the travel distance that is necessary for the ring to vary the manipulator 
stiffness. Based on the limits of the geometry a taper angle of approximately 10° 
is the highest practical and results in a ring width requirement of at least 3.25 
mm. This geometry would be a challenge to fabricate. However, the taper angle 
could he reduced and the ring size increased to enable cheaper manufacturing. 


6.4 System Analysis 

6.4.1 Theoretical System Analysis 

A model of the system was developed in MATLAB/S1MULINK ®(The Math- 
Works Inc.) as a combination of the models described in the previous sections. 
From this model it was possible to theoretically evaluate the output of the system 
assembly under various input conditions. 

First, the system response was modeled using the method of Section |ti.2,l| for a 


two dimensional single segment manipulator in planar bending operated with a 
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single active actual or by rainping the control input (variable orifice gap) from 
the open position to the closed position. The antagonistic valve orifice remains a 
constant flow restriction based on the relative orifice sizes since the orifice gap is 
large enough to not aifect the flow restriction. Material and geometric properties 
are as stated in the previous section. In thus case, the antagonistic actuator is 
in its deflated state and provides only resistance against the action of the active 
actuator due to elastic stretch of the internal bladder material. The result is shown 
in Fig. [6.14| where actuator pressure denotes the pressure at the fet'd point to 
the activated actuator. P„. and tip orientation is the angular orientation of the 
distal tip of the manipulator. The tip orientation is zero when the manipulator 
is straight and has no deformation. At this nominal position, the actuators are 
in their undeformed state with lengtlis of -10.35 min. The plot shows that in the 
absence of external loading, the tip orientation approached 180 °. corresponding 
to full retroflexion ns specified in the design requirements, as the pressure at 
the actuator feed point approaches 570 lcPn, consistent with the discussion in 
Section [0-3.21 h* the fully retroflexed state, the active actuator has contracted 
to a length of 32.5 mm with an applied contraction force of 5.07 N. while the 
antagonistic actuator has elongated to a length of -18.2 mm corresponding to mi 
applied load of 4.26 N. This result shows that the desired beam deflection can he 
achieved with approximately 83% of the supply pressure when no load is applied to 
the manipulator. However, the pressure to achieve retroflexion could be reduced 
by increasing the nominal length of the actuator. Further, since the primary 
contributor to the active actuator pressure requirement is the force applied by 
the antagonistic actuator, this requirement could be reduced by increasing the 
downstream orifice diameter such that the nominal pressure in the antagonistic 
actuator can reach a lower value. The simulation shown here was stopped once 
the angular orientation reached the design goal of 180° as this point represents the 
design limit of motion based upon the elongated length of the actuators. It was 
found that the dynamics of this model was limited by the flow rate of the valve 
rather than the inertia of the manipulator due to low mass and low moment of 
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inertin at the distal end. assumed to be on the order of 10.2 mg mid 10.2 rng'mnr . 
respectively. The response shown here For the entire motion of the manipulator 
was completed in less than half of one second through the model prediction. This 
time was determined to be sufficient for most applications, especially true when 
considering that most applications will not require full retroflexion. 



Figure 6.14: Plot of t ool t ip orientation versus actuator pressure for planar bend¬ 
ing using the two dimensional model. 


The range of mot ion for a single segment manipulator with no external loading 
Is then defined by the sweep between the undeformed configuration and the fully 
deformed configuration as shown in Fig. 6.15 This sweep range from the planar 
approximation model can be rotated around the longitudinal axis of the unde¬ 
formed manipulator in order to project the workspace into the third dimension. 


Using the same analysis for the three-dimensional case where a single actuator 
Is activated and two antagonistic actuators provide backside resistance by means 
of elastic stretch, it was found that the tip orientation (relative to bending in 
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Figure 6.15: Plot of planar beam deflection for increasing actuator pressure with 
no external load for the two dimensional model. 


the direction of the activated actuator) followed a similar curve as for the two- 
dimensional case as shown in Fig. but Instead of making it to 180° at 570 
kPa the pressure required was only 455 kPa. This result is reasonable despite the 
fact that there iire two antagonistic actuators since the termination points of the 
passive tendons are half the distance from the neutral avis when compared with 
the two-dimensional case as shown in Fig. |6.7| This means that the elongation 
of the passive actuators due to bending of the beam is reduced and thus the 
proportional antagonistic force would be slightly lower as well. The two passive 
actuators each elongate to a length of approximately 44.2 mm with applied loads of 
3.11 N. Similarly, the active actuator contracts the same as in the two-dimensional 
model but with a contractile force of 4.52 N. 


It. is also possible to increase the stiffness of the three actuators simultaneously 
through the use of out-of-plane positioning of the valve flapper ns discussed in 
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Figure G. 16: Plot of tool tip orientation versus actuator pressure for planar bend¬ 
ing using the three dimensional model. 


Section |5.5| This has the effect of increasing the nominal pressure in the actuators 
when the flapper is centered between the three downstream orifices. The flapper 
was initially positioned such that the nominal gap for each actuator was 0.18 
mm. The flapper was then moved towards the downstream orifice for a single 
actuator, as shown in Fig. |6.l7| while simultaneously moving away from the 
downstream orifices corresponding to the two antagonistic actuators (see Fig. 


G.18a( |. The rate of change of the flapper position was such that the tool tip 


orientation reached 180° within one second as shown in Fig. 0.18l> This change 
in flapper position increases the pressure for one actuator while simultaneously 
decreasing the pressure for the two antagonistic actuators (Fig. |6.1ba^ which 
results in a contraction of the pressurized actuator and lengthening of each of 
the antagonistic actuators (Fig. 6. Further, the pressure in the antagonistic 
actuators decreases until the corresponding gap size is large enough to have no 
effect on the How restriction. At tills point, the pressure applied to the actuator 
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is dictated by tin* fixed diameters of the upstream and downstream orifices. 


\ 



Figure 6.17: Illustration of valve flapper in-plane position change where dashed 
lines correspond to orientation of each downstream orifice. 




Figure 6.18: Plot of (a) gap size for both the active and antagonistic actuators 
mid (b) tool tip orientation versus time. 


Looking at the bending range in the three-dimensional case by individually acti¬ 
vating each of the three act uators produces a similar result to Fig. |6.1o 


m 


that the 


beam deforms approximately along a circular arc in the direction of the activated 


actuator ns shown in Fig. 6.20 The umbrella like shape of the curve defines the 


workspace of a manipulator segment in the unloaded case. 
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Figure 6.19: Plot of (a) actuator pressure and (b) actuator length versus tool tip 
orientation for simulation with increased actuator stiffness. 



Figure 6.20: Plot of beam deflection for varying actuator pressures applied to each 
of the three actuators with no external loath 
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Activating on addit ional actuator has t lie effect of revolving I lie manipulator about 
its long axis by changing the resultant direction of the applied moment. An 


example of this is shown in Fig. |6.21| Here the pressure applied to each actuator 
was ramped towards .'100 kPa with a small time delay between the activation of 
the first and second actuators, as shown in Fig. |6.22| The result shows that 
the manipulator deformed in planar bending as before prior to the activation of 
the second actuator. Upon activation of the second actuator, the manipulator 
exhibited out of plane deformation about the long axis with increasing magnitude 
as the pressure in the second actuator increased. At the end of the simulation, 
the orientation of the manipulator tip measured relative to the positive 2 -axis was 
approximately 187°, slightly greater than the retroflexion condition. This was 
achieved at a much lower pressure t han was calculated under previously simulated 
conditions using a single active actuator due to the use of two actuators to achieve 
manipulator bending. 



Figure 6.21: Simulation of three dimensional beam deflection through activation 
of two actuators. 
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Figure 6.22: Actuator pressures for each of the activated actuators. 

Tlie application of external load on the end of the manipulator can produce in¬ 
teresting results as the beam shape and the load on the end effector can change 
depending on current beam configuration and direction of the applied load. An 
example analysis of the effects of external load is shown in Fig. |ti.2.i| First, the 
manipulator was deformed to the configuration shown by the dashed line through 
the application of 2711 IcF'a to a single actuator. Next, a 0.5 N load was applied in 
the positive longitudinal direction as indicated in the figure. It is interesting to 
note that the beam tip not only deflects in the longitudinal direction due to the 
applied load blit also shifts in the negative transverse direction as the beam ad¬ 
justs to a new, stable configuration. This is a result of the unactuoted manipulator 
degrees of freedom which passively determine the beam shape when the manipu¬ 
lator is subject to a load. A second similar analysis was performed in which a load 
was first applied to the beam with no actuator inflation followed by an increase in 
actuator pressure, the results of which are shown in Fig. |6.21| An external load of 
5 N was applied in the negative transverse direction causing the beam to deflect 
to a fully extended configuration as shown by the dashed line. Then the pressure 
in the actuator which contracts the right side of the beam was increased to 700 
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kPa. This resulted in a recovery of some of tlie initial deflection experienced by 
the beam, moving the load approximately 20e mm back in the positive transverse 
direction, indicating the ability to move a load of this magnitude. 



Figure 6.23: Calculated beam shape for a beam deflected by a single actuator 
with anil without an end load applied. 


6.4.2 Computer Rendered System Assembly 

Using the dimensions of the design scale system, a computer rendering was pro¬ 
duced to give a visual of how the components could be assembled together (lull 
model available at [TTjTJ]). Contained within the body of the endoscope are the 
valves {plug-based design) and actuators uecessary to operate the manipulators, 
as shown in Fig. |(>.2ri| The actuators for each manipulator segment are arranged 
in series t:o maximize the available space wit hin the endoscope working channel. 
Arranging the actuators in series requires that fluid transfer lines be used to dis¬ 
tribute the How of fluid from the valve to each actuator. Additionally, retainers 
are used to provide a fixed plate to which the actuators can be anchored such that 
AMA contraction created through pressurization of t he actuator is applied to the 










Figure G.24: Calculated beam shape for a beam deflected by an end load with and 
without n pressure applied to the actuator. 


connected tendon. Each retainer includes slots to allow the anchored end of the 
actuator to change position as it inflates. This change in position also requires 
that the fluid transfer lines be flexible such that they can move with the actua¬ 
tors. The retainer slots also provide passage for the fluid transfer lines for (list ally 
located actuators. An example of the working channel cross-section is shown in 
Fig. |<>.2(>| where a siugle actuator is shown in its fully inflated state alongside two 
deflated actuators. The fact that the fully inflated actuator does not occupy the 
full working channel diameter allows enough remaining space for the fluid transfer 
lines as well as channels for tendon routing without interference. 


The valves are connected to the existing suction and irrigation lines within the en¬ 
doscope (see Fig. |ti.27^ to eliminate the space requirement for the extra hydraulic 
lines. Putting it idl together, two, two-segment manipulators might maneuver 
within the defined workspace as shown in Fig. (>.2<S In this figure, an endo¬ 
scope is shown with two tool channels; each containing the components for two. 
two-segment manipulators. The components for each manipulator segment are 
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Figure (i.25: Computer rendering showing the full manipulator assembly in (a) 
perspective view and (b) side view. (c> the internal components used to operate a 
single segment of one manipulator, and (d) t he connection between the valve and 
a single actuator shown in its inflated state. 
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Figure 6.2G: Illustration of the cross sect ion of the working channel of an endoscope 
showing three actuators arranged in parallel, one of which is fully inflated and the 
routing channels for the manipulator tendons. 


arranged serially within the working channel and the tendons for the proximal 
segment are routed aJong the outer radius of the working channel as shown in Fig. 
G.2(i|in order to avoid interference with the components lor the distal segment. 


0.4.3 System Prototyping and Testing 

Limitations due to the valve not performing ns intended as discussed in Chapter[5] 
prevented the hydraulic surgical robot system from beiug assembled in a controlled 
fashion. However, it was possible to test the interaction between the artificial 
muscle actuator and the continuum beam manipulator by manually inflating the 
actuator usiug a pressure instrumented syringe (HARD Caliber Inflation Device) 
us shown in Fig. |ti.-!)| or through manual manipulation of the ring-based valve. 
This made it possible to get a better understanding of the behavior of t he actuator 
and manipulator in response to changing pressure. For this testing a prototype 
of a manipulator backbone with attached actuators was produced from available 
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Figure G.27: Illustration of the valve within the working channel of an endoscope 
with channels providing connection to the suction and irrigation lines. 

materials. The backbone was made from a length of nitinol wire with a diameter 
of 1 mm (£7 = 1.37 • 10 ~ a N ■ in 2 ) and a length of 185 mm. The disks used to 
route the tendons had a radius of (i mm at the location of the tendon. 

As it was not possible to produce a muscle actuator at design scale, a version 
of the muscle actuator was produced using latex surgical tubing as the internal 
bladder and pie-manufactured nylon mesh, with an inflated outer diameter of 8.8 
mm, as shown in Fig. |4.lti| and described in Section |4.4| This muscle actuator 
was found to perform ;is predicted analytically in terms of contraction range as 
demonstrated in Section [I.-l| For tins to be possible, care was taken to coordinate 
the elongation of the outer mesh with the deflated length of the surgical tubing. 
In particular, it is important to ensure that minimal clearance exists bctw r een 
the outer mesh and the bladder in order to prevent the bladder from crossing an 
elastic instability point (aneurysm formation) |1G4| . To do this, it is necessary to 
stretch the inner bladder and outer sheath together. It. is also necessary to ensure 
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Figure 6.28: Computer model of design scale, two segment per manipulator system 
with workspace shown. Shown with example end effector tools including a grasper 
nnd a scissors (8-l| . 



Figure 6.29: Photograph of a pressure regulated syringe and manipulator testing 
apparatus. 
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that the actuator is not taut when connecting it to the tendons thus allowing the 
actuator to elongate when an antagonistic actuator Is activated. In practice this 
has the effect of slightly increasing the pressure necessary to inflate the bladder 
due to an increased level of wall stress that must lie overcome. To minimize this 
effect it is ideal to select a bladder material with high elasticity (compliance >0.05 



la the planar case with a pressure up to 1.03 MPa applied to only a single actuator 
using the instrumented syringe, the actuator and manipulator were found to per¬ 
form as predicted by the model. The manipulator behaves as in the planar ease 
presented in Section 3.3| exhibiting no out-of-plane deformation and reproducing 
a circular arc for which the radius of curvature decreases with increasing actua¬ 
tor pressure. This result is expected ns the inactive actuators provide symmetric 
opposing forces on the antagonistic side of the manipulator. The forces from 
these actuators reduce the output from the active actuator (but not significantlv) 
and actually improved the lateral stability of the manipulator. A photograph of 
the manipulator under these conditions Is shown in Fig. |G.30[ The pressure re¬ 
quired to achieve full retroflexion was recorded as approximately 281) kPa This 
pressure is subject to variability according to the nominal actuator length, actu¬ 
ator construction parameters such as amount of pretension, and geometry of the 
manipulator. When jui additional actuator was activated, the curvature of the 
manipulator revolved around the longitudinal axis to the new equilibrium orien¬ 
tation. In the absence of external loading, the manipulator shape remains planar 
as predicted by the model. 


It was also possible to evaluate system performance through manual manipulation 
of the ring-based valve prototype. Because the actuator used here was larger than 
the design scale, it was necessary to increase the size of the ring-based valve 
orifices to accommodate the higher flow rate requirement of 5.0 —. Using the 
same procedures described in Section [5.3[ an upstream orifice diameter of 0.74 mm 


was necessary. Equation |5,7| then results in a downstream orifice of diameter 1.3 
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Figure 6.31): Photograph of manipulator backbone under load. 


mm in the fully open condition. 1 'sing this diameter t he necessary gap. g. between 
the orifice and the flapper under the fully open condition is 0.4 nun. These valve 
orifices for both the upstream and downstream orifices are sized to achieve both 
sufficient flow rate ami appropriate valve control for this valve prototype. Also, for 
this prototype, a taper angle of 3° was used in order to make manual manipulation 
of the valve ring possible. 

The outlets of the ring based valve design were connected to the some muscle 
actuators and pressure was supplied using the same pump setup as shown in Fig. 
5.1o| with a constant supply pressure of 689 kFo. To test the interaction between 
the ring-based valve design and the actuators, the ring flapper was manually 
manipulated to differentially open and close the flapper nozzles. When one of 
the flapper nozzles was fully blocked by the ring, the actuator was inflated in 0.35 
seconds on average. This confirmed the design goal of 0.5 seconds for full actuator 
contraction which was used to size the orifices. Further, it was found that two 
actuators could be inflated by positioning the near side of the ring between the 
two corresponding nozzles. Finally, by adjusting the out-of-plane position of the 
ring it was possible to adjust the relative neutral state between the three actuators 
(stilfness control). 
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Being able !o vary tlie stiffness of the actuators ns a unit provides some measure of 
control over the baseline rigidity of the manipulator jus deflection in any direction 
requires more load to overcome the actuator loads. To get. an idea of what the 
effect of this variable stiffness is, the actuators were inflated simultaneously while 
using a pressure gage. A lateral load was then applied to the manipulator using 
a spring scale and the deflection of the manipulator was measured as shown in 
Fig. |6.3l| The test was carried out at two positions along the length of the beam. 
Position 1 was located GO mm from the fixed base and Position 2 was located 185 
mm from the fixed base. The results of this test are shown in Table 16.II 

Table 6.1: Results from manipulator stiffness measurements. 


Inflation 

Loarl Position i 

Load Position 2 

Pressure 

Stiffness ( 7 ^ 777 ) 

Stiffness (—) 

0 kPa 

0.08 

0.014 

345 kPa 

0.12 

0.018 

689 kPa 

0.19 

0.024 


Supply 



Figure 6.31: Illustration of test setup for measurement of manipulator stiffness. 











6.5 Summary of System Assembly 
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The system model described here incorporates each of the individual components 
that make up the hydraulically powered surgical robot design. The model has 
been formulated in both two-dimensional and three-dimensional versions. The 
two-dimensional version allows for simulation of the assembled surgical manipu¬ 
lator subject to planar deformation and is suitable to simulate manipulator mo¬ 
tion when no out of plane external loading is applied. The model includes the 
analytical descriptions of each of the robotic components described in previous 
sections. Where necessary, these components are coupled together through the 
relevant hydraulic and geometric relationships. The two-dimensional model was 
able to demonstrate that the design goal of 0.5 seconds to achieve full retroflexion 
is possible with the included valve design. The simulation suggests the ability 
to manipulate a load of 5 N. However, the magnitude of the distance that this 
load could he moved is not large. This result: may indicate the need to explore 
higher actuator pressures to achieve greater capability depending on the specific 
intended task of the manipulator. Further, at higher operating pressures it would 
be possible to reduce the size of the actuators thus making parallel arrangement 
of the actuators feasible. This would in turn reduce the overall length of the ma¬ 
nipulator components within the endoscope as shown in Fig. |6.32 For example, 
increasing the operating pressure to 3.5 MPa would allow the actuator diameter 
to be reduced to a theoretical value of 1 mm to achieve the same load capacity. 
However, the increased pressure would necessitate a denser braid which would re¬ 
duce the contraction range and thus increase the required length of the actuator. 
Further, the use of higher operating pressures and higher act uator loads could be 
combined with a st ilfer nnuiipulator backbone in order to reduce transverse deflec¬ 
tion due to external load. Fortunately, the loads common to many surgical tasks 
are significantly smaller and thus more achievable. As compared against similar 
surgical robots presented in the literature, the system models suggest that appro¬ 
priate criteria for robotic manipulation speed and load capacity can be achieved 
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without the use of external actuators or motors. 

One challenge with the system model is computational effort required to compute 
the output. For both the two-dimensional and three-dimensional system models, 
the use of an optimization algorithm for modeling the beam response prevented 
real-time operation of the model. Pre-computation of the beam output for a 
variety of end loads was found to improve computation time but not enough to 
permit real-time simulation. Further enhancement of the system model will allow 
this simulation to be used as a plant model for system controller design. 

Testing of the system using manual inflation with an instrumented syringe con¬ 
firmed the model predictions for beam shape achieved during actuator contraction. 
Further, in the absence of external load, the manipulator exhibited planar defor¬ 
mation which is also consistent with the model prediction. When the actuator 
and backbone prototype was connected to the ring based flapper valve, it was 
possible to measure the contraction time of an actuator. This result confirmed 
the performance of the valve prototype which was designed to achieve full actu¬ 
ator contraction within a maximum of 0.5 seconds. Further, by simultaneously 
inflating each of the three actuators which are used to operate the prototype seg¬ 
ment. it was possible to get an idea of the effect of baseline actuator pressure on 
manipulator stiffness. These results indicated that greater pressures do indeed 
result in greater stiffness. This situation is useful in the case of increasing ma¬ 
nipulator stability during a given task. However, it comes at. the expense of load 
manipulation ability as the actuator acting in the direction of intended motion 
would have to act against the antagonistic actuators in addition to the load itself. 
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Endoscope 


Manipulator Segment 
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Figure 6.32: Computer rendering showing internal system components for a single 
manipulator section in (a) an exploded view, (li) assembled in plate, and (c) 
assembled within the working channel of au endoscope. 









Chapter 7 

Conclusions and Future Work 


7.1 Review 

The hypothesis of the research presented here was that hydraulic power can he 
scaled and adapted for use in embedded surgical robots for minimally invasive 
surgery. To test this hypothesis, it was first necessary to identify the requirements 
for novel minimally invasive surgical robots ami the needs that exist in the field 
through experimentation and a review of existing and developing devices. A 
vision for the proposed solution was then laid out and the major considerations 
specified. Under this vision, it was necessary to consider the specifics of each of 
the relevant hydraulic components that must be present to assemble a surgical 
robot. Finally, through analysis and testing of these components several things 
were learned about the realization of a hydraulic surgical robot and provided much 
insight to further developments in this field. 

While the future of NOTES is in question by many, the progression of gastroin¬ 
testinal surgjcal intervention methodologies over the years traces a clear path 
based around the continuous reduction in incision number and size. Following 
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this progression through current developments makes a compelling argument that 
single-incision surgery and eventually natural orifice surgery are inevitable if the 
surgical practice is to continue to progress. However, for this to be possible, the 
tools that enable these procedures must be available and fully qualified. The 
transgastric approach for NOTES Is often considered the most challenging from a 
technological viewpoint and tints is the subject of much research and development. 

The first chapter began by describing the history* of surgical intervention and 
the development of minimally invasive surgical techniques. The progression of 
the surgical practice in relation to the abdominal cavity from traditional open 
access techniques to emerging natural orifice procedures was presented. The tools 
that will enable these novel procedures arc developing rapidly. A review of these 
technologies was conducted and the strengths and weaknesses described as a means 
of developing the design path for this research. The rest of the chapter summarized 
the motivation for this research, stated the problem to be overcome, and laid out 
a vision for the design of a hydraulic surgical robot for minimally invasive surgical 
int erventions. 

The next chapter described the design requirements for the surgical robot designed 
for use in natural orifice transgastric endoscopic surgeries. First, the requirements 
for load capacity at the tool tip were investigated. As it was found that little infor¬ 
mation existed in the literature to describe this requirement, a set of experiments 
was conducted to quantify this value. The results of these experiments were pre¬ 
sented. Following this, a review of the literature was used to evaluate design 
requirements for the surgical workspace as well as considerations for safety and 
biocompatibility for the design. These considerations helped to guide decisions 
made during the design process. 

Chapter [3] provided a detailed description and analysis of the selected backbone 
structure. A model using the principle of minimum potential energy and inverse 
kinematics was presented to predict the behavior of the backbone. The model was 


181 

then experimentally validated to demonstrate that it is capable of producing accu¬ 
rate beam shape predictions under most circumstances. Recommendations were 
provided as to how to best utilize the model under the prescribed assumptions. 

In Chapter [T] the selection of an actuator to drive the manipulator was presented. 
The performance of the actuator was predicted analytically and confirmed exper¬ 
imentally. Further, an optimization algorithm was used to determine the ideal 
parameters necessary for fabrication of an actuator which maximizes range of mo¬ 
tion within specified const raints placed on actuator diameter and force capability. 
This was accomplished through modeling of the stresses present within the actua¬ 
tor using finite element analysis. It was then demonstrated that the requirements 
placed upon the actuator can he met using the given design. 

Chapter [^described the method of providing control input to the system by way of 
a control valve. To meet the ueeds of the presented device within the constraints 
of the application, a novel valve design was conceived. Design parameters that are 
necessary to achieve desired valve performance were calculated. Two realizations 
of the design were fabricated at a larger scale to provide a means of evaluating 
the design performance and confirming the hydraulic principles used to develop 
the design. The results presented suggest that the design can be scaled to a size 
appropriate for application to the field of minimally invasive surgery. 

In Chapter [C} a description of the assembly of the novel manipulator was provided. 
A model is shown which combines each of the manipulator components to allow 
for prediction of manipulator performance. Design specifications were given for 
each of the components that make up the assembly such that the design Is theo¬ 
retically capable of meeting the requirements set: forth in Chapter[2j Additionally, 
results from testing of a manually operated manipulator assembly were presented 
to demonstrate the efficacy of the design. It was shown that it was possible to 
achieve the desired performance objectives using t his desigu. 


7.2 Conclusions 
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The Held of surgical robotics is rapidly moving and t lie development of new surgical 
robots is taking place continuously from both academia and industry. These 
developments are motivated by both the expansion of laparoscopy as a method of 
surgical intervention and the exploration of new surgical methods. The only robot 
performing abdominal procedures on a wide scale is the da Vinc i Surgical System 
bv Intuitive Surgical. While the da Vinci robot has helped gain acceptance for 
robot surgery, it has several limitations such as cost, size, and mobility which have 
made it inaccessible to smaller, rural hospitals and for use in remote or battlefield 
situations. Further, the da Vinci system is not currently well adapted for newer 
natural orifice access techniques. However, recent developments have suggested 
that the company is working towards expanding the platform for these types of 
surgical developments (ilia) . Beyond the cla Vinci, a wide variety of novel surgical 
robots have been described in the literature. A number of these were presented in 
Section [1.1.5[ While each of these devices approach the problem from a different 
direction, they all face the same challenges in balancing between scale and robotic 
capability. It is from this balance that the hypothesis of adapting hydraulic power 
for use in embedded surgical robotics was born. 

After a review of the literature, it was found that the load requirements for per¬ 
forming procedures in the abdominal cavity are not well documented and much of 
the information available was in the form of estimation and speculation. There¬ 
fore, a series of experiments was conducted to help quantify these requirements. 
The common surgical tasks that were identified as likely having the highest tool 
force requirements included organ or tissue manipulation, suture knot tightening, 
and suture breaking. The results of this investigation contribute to the relatively 
small pool of available quantifications of tool force measurements for abdominal 
surgical tasks that con be found m the literature. This hel|xs to lend more rigor 
to mi area that has largely been taken for granted when designing novel surgical 
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robots. Of particular note, when reviewing the literature available on robots in 
development, it is evident that many of these robots do not achieve this load 
capacity specification. This may limit the practical applicability of these novel 
robots to lighter duty activities. The results for tool force load combined with 
considerations for workspace, safety, and biocompatibility assembled from avail¬ 
able resources helps to define the requirements for surgical robots and guide the 
design process. 

A continuous backbone was selected as the structural member for the manipula¬ 
tor. This selection was made because a continuous backbone provides enhanced 
usability in complex environments or for fine motor manipulation. However, a 
downside to the use of a continuous beam manipulator is that its behavior under 
load cannot be easily predicted due to the complexities of spatial large deflec¬ 
tions. While a variety of options for the analytical prediction of beam behavior 
have been demonstrated in the literature, the majority of them require a pro¬ 
hibitive level of computational effort to arrive at a solution. As an alternative to 
available solution methods, an inverse kinematics model based on the assumption 
that the beam would assume the configuration with the lowest potential energy 
was developed. The two-dimensional version of this model was validated exper¬ 
imentally and found to be reasonably accurate for most conditions. This model 
makes it possible to predict the beam shape and end forces/moments for a given 
set of end position and orientation parameters thus making it possible to use these 
predictions for further design. 

The continuous backbone is driven by hydraulic actuators connected to thin wires 
which are routed through evenly spaced disks along the backbone. The disk at 
which these wires terminate determines the location where the majority of the 
actuator load is applied to the backbone. For this application, artificial muscle 
actuators were selected rather than traditional hydraulic actuators due to their 
ability to achieve greater force output for the same pressure and size at the expense 
of contraction range and a force output which varies with actuator length. This 
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realization is important for a severely .size constrained application. Further, the 
response of the muscle actuators is similar to that of human muscle tissue in 
tonus of compliance, aiding to the safety of the manipulator during an operation. 
Art ificial muscle actuators are typically used in pneumatic applications; however, 
the use of water or saline results in similar actuator performance and less risk in 
the case of a burst. It was not possible to prototype a muscle actuator at the 
design scale, and thus huger actuator prototypes were produced and confirmed 
the predicted actuator performance. For use in the surgical manipulator, the 
artificial muscle actuator works well since the limitations in diameter can be met. 
Further, the surgical manipulator hns less stringent limitations for the length of 
the actuator and this helped to make the muscle actuator a feasible solution. 

Flow within the system to each of the actuators is manipulated by a novel flow 
control valve. To arrive at a control solution for the surgical manipulator, it: 
was necessary' to find a control valve or combination of control valves that could 
regulate t he pressure supplied to each of the three actuators that articulate a given 
manipulator segment. For this to be feasible, the valve must lit within the available 
diameter. Hydraulic valves of this size were not found to be available commercially. 
Therefore, the design of a novel control valve was presented. The requirements for 
the valve were that it be capable of proportionally regulating the pressure supplied 
to each of the actuators while meeting the constraints of size. Ideally the valve 
would provide for both differential operation of the actuators and stiffness control 
of the three actuators together. Several designs for this actuator were considered. 
The two designs presented in Section [o] make use of three dimensional motion of 
a manipulated part (plug or ring) that acts as a '‘flapper". Planar motion of the 
flapper opens and closes the flapper nozzle for each of the actuators differentially. 
The flapper can be positioned such that pressure in any one actuator is increased 
individually or in pairs of two. This method of control is appropriate as the 
actuators would typically be controlled in this way for manipulation of the robot. 
Axial (out-of-plane) motion of the flapper can be used for stiffness control of the 
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three actuators together. A scaled prototype of each of these valve designs was 
produced. Testing of these prototypes revealed that the actuator pressure could 
he manipulated as expected through manual operation of the valve. However, 
for the plug based design, the available range of this pressure control was smaller 
than expected most likely due to imperfections in the rapid prototyping based 
fabrication method. While it was not possible to test the use of the driving coils 
in valve operation, test ing of t he coils independently found that the small diameter 
coils were capable of manipulating the plug without overheating when a steel core 
was used. Further, when submerged in water, the current capacity of the coil 
was increased significantly. Thus, in practice the coil should provide sufficient 
input for valve manipulation. The ring based design proved much more capable 
at increasing the pressure at the actuator feed points. This is likely due in part 
to the larger manipulated component which is more easily handled and the fact 
that the valve could be fabricated out of aluminum through traditional machining 
techniques. With the ring valve connected to the artificial muscle actuator it was 
found that, on average, the actuator could be inflated in 0.35 seconds when the 
flapper nozzle is fully blocked by the ring. This suggests that the orifice size for 
the desired flow rate is appropriate and thus this valve design shows feasibility of 
the valve concept. 

The design of a hydraulic surgical robot has proved to be a difficult endeavor with 
many challenges. The application of endoscopic robots to transgastric natural 
orifice surgery introduces severe constraints in terms of the size of the device, 
load capacity, and safety and biocompatibility considerations. While the design 
of many devices have been proposed in the literature, the majority of them do 
not meet several of these requirements and do not represent a viable solution for 
this type of surgical procedure. In theory, the application of hydraulic power to 
this field presents several desirable characteristics such as improved power den¬ 
sity. greater stiffness, and faster response than comparable electromechanical or 
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pneumatic designs. To map these principles to the held of surgical robotics re¬ 
quired significant effort as hydraulic components of an appropriate size are not 
available commercially. Therefore it was necessary to consider alternative solu¬ 
tions and novel designs for the various robotic components to assemble a design 
that was capable of meeting the design specifications while taking advantage of 
hydraulic power as a system input. Through these efforts, much new information 
was produced and developments resulted from the research presented here. The 
considerations of hydraulic power for surgical robots were evaluated and are appli¬ 
cable to other surgical tasks where hydraulic power may be used advantageously. 
A quantification of the load requirements for surgical robots performing abdominal 
procedures was produced which will provide a guide for other researchers develop¬ 
ing surgical robots. These values have typically been reported anecdotally in the 
literature and thus are a valuable resource for the field. An alternative model for 
predicting the behavior of continuum beams under load was developed to provide 
a inverse formulation for computing beam shape and end loads. This is useful 
as continuum beams are widely used for minimally invasive surgical manipulators 
as well as in a wide variety of other applications. Finally, a novel valve concept 
and two possible designs realizing this concept were developed. These valve de¬ 
signs facilitate control over the three actuators in an antagonistic arrangement. 
Further, the valve designs enable proportional control of the three actuators at a 
size scale not commercially available. In summary, the design of a novel hydraulic 
surgical manipulator as a summation of its parts has been performed and it has 
been demonstrated that fluid power can be used advantageously for the design of 
embedded surgical robotic tools for minimally invasive surgery. The pursuit of this 
research has provided many unique challenges and the work presented here has 
addressed many of them, as well as laid the foundation for future developments 
in the application of hydraulic power to the growing field of surgical robotics for 
minimally invasive surgery. 



7.3 Future Work 
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This dissertation has presented the design of a novel hydraulic surgical manip¬ 
ulator as a contribution to the rapidly growing body of knowledge surrounding 
the hold of minimally invasive surgery and the tools necessary to perform such 
procedures. The following sections outline several opportunities for which future 
work could proceed, building upon the foundation of work presented here. 


7.3.1 Prototype Iteration 

As previously discussed, difficulties with the fabrication of the valve prototype 
prevented testing of the manipulator in a controlled manner. The next steps would 
then bo to refine the valve prototype fabrication to achieve a working version. For 
the plug-based version of the valve, the primary limitations of the prototype were 
that the complicated geometry of the valve body required the use of a rapid 
prototyping process. Therefore, a fabrication of the valve at design scale should 
be explored. To do this it will likely be necessary to look at ways to simplify the 
geometry for manufacturability. Further, the use of alternative manufacturing 
processes may help to realize the design. In particular, the tetrahedral seat for 
the plug is the most difficult feature to produce. The use of electrical discharge 
machining may make the tetrahedral shape possible. With these considerations 
it is believed that producing a working prototype of the plug-based valve design 
will be possible. 

The ring-based valve prototype demonstrated that the valve design was feasible 
and worked as expected when manually manipulated. The future goal for this 
valve prototype would be to achieve greater miniaturization through precision 
manufacturing. The primary benefit of miniaturization, aside from meeting the 
size requirements, would be to make the ring small enough in mass that it can be 
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manipulated with either an electromagnetic or a piezoelectric actuator. As the 
scale of this valve design decreases, the tolerances of the machining operations 
become increasingly tight. Therefore, the use of precision micro-machining will 
likely be necessary to make t he design scale prototype possible. 

Once the valve design prototype has been produced, it will be passible to integrate 
it with the other components of the manipulator. Fabrication of a custom internal 
bladder will permit miniaturization of the artificial muscle actuator. This can be 
accomplished by contracting with a custom extrusion company to produce tubing 
at an appropriate diameter ;ls no solut ions were found to he available off-the-shelf. 
This will permit further testing of the system prototype to ensure that the design 
requirements can be met at design scale. Further, this will make it possible to 
test system response and find an appropriate mapping from user input to desired 
manipulator action. 


7.3.2 Sensing and Control 

The \Lse of active sensing at the distal end of the surgical manipulator Is one area 
that garners much attention within the field of surgical robotics. However, there 
is not yet a consensus regarding the method in which it: is employed or even the 
necessity of its use. When it comes to force-feedback to provide a sense of touch 
for the surgeon who is manipulating the robot, the answer is unclear. Surgeons 
often proclaim that the presence of force-feedback is a necessity for robotic min¬ 
imally invasive surgery, yet surgeons regularly perform procedures without this 
information whether it. be with hand-held tooLs or with commercially available 
robots such as the da Vinci [80] . Some work has explored the use of visual or 
auditory feedback to relay force information or even to automatically avoid ob¬ 
jects all together [581 [Ilk'll 11(17] . However, the effectiveness of these methods seem 
to have significant variation among individual surgeons [1<18[ . More investigation 
into the methods by which t his information is relayed to the surgeon and the costs 


or benefits of proving it Ls necessary. 
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Another use for sensing on the distal end of the robot is to provide position 
feedback to the controller. This form of feedback can be provided through active 
measurement of the tool joints (not possible with a continuum based robot) or 
possibly through other means such as strain gauge based measurements on the 
manipulator backbone. It may also be possible to gather this information tltrough 
measurement of the actuator contraction or pressure. Further, some work has 
explored the use of algorithmic based optical tracking of the tool using either MR I 
or some other imaging modality [lt>9| . Current commercially available robots make 
use of the surgeon as a means for correcting absolute positional error produced 
by the robot, which has been measured to be as large its 25 mm [79j- In this way 
the surgeon is able to compensate for any realized error due to tool behavior or 
unexpected loads and increase or decrease the magnitude of his or her input iis 
appropriate. 

It is also necessary to consider the method by which the control valve will be op¬ 
erated. To achieve this, a cont rol algorithm must be developed. This development 
can likely build upon similar work done in areas such its human power amplifiers 
[170111711 . which are an appropriate solution for many applications due to their 
high power density and ability to accept direct lnunan manipulation. Addition¬ 
ally, considering the principles of teleoperatiou, it is then possible to produce an 
output that can be scaled in both force and displacement.. Tools such ns these 
enable a virtual physical connection between the human operator and the task 
being performed. Additionally, through the use of force-feedback it Ls possible to 
provide an intuitive, haptic interface for the operator where they can receive in¬ 
formation about the loads being applied to the tool and adjust the tool operation 
appropriately. 

Similar techniques have been applied to t ranslate lmm;ui manipulation to a smaller 
scale system using a human power attenuator. An attenuator allows the human 
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operator to manipulate the controls on a scale that is more natural while the tool 
performs the same motion at a more precise scale. Much work has been performed 
in this area along with its application to surgical robotic systems for minimally 
invasive surgery |72J IT721 1173] . 

The advent of minimally invasive surgery has thrust this Held into the forefront. 
As both surgeons and patients seek out new technologies that use smaller incisions 
and minimize the risk of injury, the necessity for increasingly small surgical instru¬ 
ments is growing. With decreasing size of surgical instrumentation, it becomes 
more difficult for the surgeons to manipulate these tools in a precise, controlled 
manner. Therefore, surgical robotic systems provide a means to achieve precise 
manipulation and force application while maintaining the surgeon's control over 
the surgical tools. 

The control objective for such a system is to accept the power input by the opera¬ 
tor at the control joyst ick and produce au attenuated power output at the surgical 
tool. Difficulties may arise in the synchronization of the input joystick and output 
hydraulic act uator due to model uncertainties inherent to hydraulic systems |174| . 
Uncertainties may he presented in the form of external disturbances, fluid leakage, 
and friction. 


7.3.3 User Interface 

A proper user interface is critical for the success of a surgical manipulator. As 
previously mentioned, the da Vinci surgical system makes use of a surgeon console 
that has gained widespread acceptance for the operation of a surgical robot. This 
type of console works well for remote operation of a surgical robot as it provides 
visual feedback to the surgeon, places the surgeon in an ergonomic position, ami 
allows the surgeon to provide input to the robot. A similar surgeon console would 
be necessary for the manipulator described in this work. Since each segment of the 
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manipulator contains two degrees of freedom, it is possible to measure user input 
with two sensors. These sensors could be of the form of rotary potentiometers. 
Using two potentiometers oriented perpendicular lo each other would permit the 
capture of both degrees of freedom at a single location similar to a joystick. This 
would also permit expansion of the user input to additional segments by adding 
addition pairs of potentiometers. By positioning the potentiometer pairs at nat¬ 
ural locations such ns t he elbow ;uid wrist, operation of the surgical robot: would 
he more intuitive. This method of user input would he consistent with the use of 
artificial muscle actuators in that the manipulator would move and behave in a 
familiar manner. This can lead to reduced training time for surgical students as 
well as reduced warm-up time before each procedure. An idealized user interface 
for control of a surgical robot would provide three-dimensional visual capability 
through a head mounted system, allow the surgeon to don and dolf the console 
at. will, and include all the necessary forms of information delivery to ensure the 
best possible surgical outcomes |fia]. 
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Appendix A 


Design Alternatives 


Manipulator Backbone 

Spherical Joint Structure 

One possible alternative solution to the continuum beam backbone that has been 
explored is a jointed struc ture as shown in Fig. [A.i| which would be driven by an 
array of tensile wires (three or four) that run along the outer edge of the joiued 
links. The jointed structure provides an option that is more easily modeled as 
the tip position can he predicted by knowing the orientation of a small number 
of segments. Further, the use of rigid links means that the resulting manipulator 
will be more precise than a similar continuum beam based manipulator since 
each there are a finite number of degrees of freedom that can each be actuated. 
The primary disadvantage of this approach is that the resulting structure is more 
complicated and thus more difficult to manufacture. 
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Figure A.l: Jointed support structure. 


Analytical Description 

Tills jointed structure Is comprised of multiple independent pieces connected 
through spherical joints, which were chosen in order to provide minimal resistance 
ngainst manipulation and rotational freedom in all directions. It was necessary to 
develop an analytical description of the proposed component in order to estimate 
how the device might perform more accurately than can be done by intuition. For 
the purposes of developing the analytical model for this structure, each piece was 
treated as an individual free body. First, it is necessary to define the geometric 
properties of the of each link and the relationship between t he links and the tensile 
wire as shown in Fig. |A.2| The geometry of each free body at the location of the 
spherical joint Is assumed to be a point for amplification. 
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Figure A.2: Sketch showing the labeled geometric references used in later equa¬ 
tions. 


As shown here, c is the angle between the edge and central axis of a link. 0 is the 
angle between a line perpendicular to the base and the link’s central axis. *y is the 
angle between a line perpendicular to the base and the force vector of the tensile 
wire in the proximal direction, and ,3 is the angle between a line parallel to the 
base and the force vector of the tensile wire in the proximal direction. Further, 
F, Ls the force applied to the tensile wire, h is the central axis height of the link, 
and .v Ls the side length of the link. Then, by geometry it can be shown that the 
angle 7 can be derived as. 



sin( 0 , + i") — sin t 
cos( 0 , + <!>) 


) 


(A.l) 


and aLso that 

A* Ji (A.2) 

where the subscript / denotes that the variable can be defined as such for any of 
the links within the joint. 


Using this notation, a free body diagram of an example three-link joint (Fig. |A.3fr 
can be used to derive the moment and force balance equations that allow the 
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joint to he fully defined in space for a given tensile force. For simplification of the 
model, the mass of the individual components .-vs well as any frictional losses in 
the spherical joints or tensile wire interactions is assumed negligible. Further, the 
interactions between the tensile wire and the free bodies are assumed to occur ns 
a point location and in only one direction. 



Figure A.3: Free body diagram of a sample three-link joint. 


For the distal link, labeled l in Fig. A.3j the balance equations are realized as, 
7: A/ = 0 =F, sin {$\8 sin(i’ + ) - F t cos Q x s cos {4< + 0 X ) 


- Fiuad sin 9\l sin 0 X - F llxlti coscost 


(A.3) 
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£ Fr.i = o = It,. 1 - F, sin ft + F w sin (9 t 
5^ F „.i = 0 = /? w .i - F, cos ft - F w cos 0 ,. (A.5) 

Additionally, for the remaining two elements, the balance equations are, 

y: A/, =0 =F, sin fts sin( ii< + 0 ,) - F, cos il,s eos( i • + 0,) 

- Ft cos( 7 i_i + sin(V + 0.) + F,sin< 7 ,-i + 0,)s cos(t- + 0.) 

— It v .,-i sin 0 ,/f sin 0 , — /? ViI _i cos H,h cos 0 , 

-+- F t i _i cos0 ,h sin 0, — ft ,.,-1 sin 0,h cos 0, 

(A. 6 ) 

F,j = 0 = /?,j-F,sinA + F,cos( 7 i _ l + 0 ,) + /? Jlii _i sin 0, - F r .,_, sin0, (A.7) 
X] i = 0 = R U' - F > ros A + F f Hin( 7 i-,+ 0 ,)-/? Vi i_, COS0,sin0„ (A.S) 
where i — 2.3 corresponding the the elements labeled 2 and .V in Fig. |A.3[ 


Conclusions 

Solving these e(|uations simultaneously allows the motion of the jointed structure 
to be predicted for a given load condition and control inputs. Ultimately, the 
model suggests that despite the advantages that this design presents in terms 
of freedom of motion, it may not: provide sufficient control over tip position for 
the application. In particular, the proximal link in the series actually deviates 
in the opposite direction from the applied tension until the following link reaches 
the limit of its rotational motion as demonstrated using the prototype in Fig. 
|.A.-I| This situation could be avoided by having low mass segments and applying 
a slight pre-tension to the wires. One additional difficult)' for the realization of 
this design would be find a robust method of fabricating the joints such that they 
could withstand repeated use reliably at the appropriate scale. 
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Figure A.4: Demonstration of problematic joint response. 


Control Valve 

MEMS Valve 

As ;ui alternative solution to the valve design presented in Section [5] a review 
of the conmion niicro-lliiidic valve solutions was conducted and no existing valve 
design was found suitable for the present application due to the limitations of high 
supply pressure, proportional control, and method of activation. Therefore, the 
feasibility of existing MEMS technologies to the application of proportional flow 
control was investigated. In order to successfully regulate flow pressures greater 
than a few thousand pascals, it is necessary to design for pressure balancing, an 
example of which Is shown in Fig. [A.~i| due to the low actuation forces that can 
be produced by most iuicroelectromechanical systems (MEMS) based applications 
[142| . Further, it is often difficult to achieve huge displacements, such as those 
that: would be necessary to achieve the desired How rates, using MEMS techniques. 
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VALVE CLOSED 


OUT 



Figure A.5: Example of a pressure balanced valve concept |H2j . 


()ne attractive solution in this space is the electrostatic comb-drive actuator, which 
lias been shown in the piLSt to be capable of large displacements relative to actuator 
size )175| . It Ls necessary’ to obtain huge scale actuator displacements in order to 
achieve a sufficient valve orifice size for the desired fluid How rate. For the desigued 
system, a maximum displacement of 60//m would be sufficient; however, greater 
displacement capability would permit a lower orifice aspect ratio which would 
yield better flow characteristics. 

A review of the literature reveals several electrostatic comb-drive actuator de¬ 
signs intended to extend the nuige over which the actuator is capable of operating 
j 176111771117811179| . Each of these solutions have inherent advantages and dis¬ 
advantages. However, of primary importance in the design of large displacement 
electrostatic actuators is the balance of displacement capability and lateral in¬ 
stability, sometimes referred to side pull-in |18tlj . Side pull-in occurs when the 
applied voltage difference across t he comb-drive electrodes is huge enough that 
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the electrostatic force between the interdigitated lingers causes them to pull to¬ 
gether, thus disabling the actuator |181) . To combat this effect, the comb-drive 
actuator must be designed such that the lateral stiffness of the suspension system 
is sufficiently high to avoid pull-in over the operating range of the actuator. 


Comb-Drive Principles 

One of the key components of the design for a novel microHuidic control valve is 
the method for achieving large displacements. In the present design, the method 
chosen for producing the mechanism displacement is use of an electrostatic comb- 
drive actuator. The comlr-drive actuator achieves motion through the application 
of a voltage difference across fixed and movable interdigitated comb fingers. This 
voltage difference produces an electrostatic held between the neighboring fingers 
and thus induces an attractive force that can be modeled us a parallel plate ar¬ 
rangement. set* Fig. |A.(»| 11751 118‘il . 



Figure A.6: Geometry of a single set of comb-drive fingers |1S‘J |. 
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Using tliis parallel plate model, with the geometry shown in Fig. [A.fit the capac¬ 
itance of parallel plates can be expressed as 


^, _ 2ue 0 h A l AQ 
so 


(A.9) 


where n is the number of fingers and is dielectric constant in free air. From 
this, the force. Fta, can lie described as presented in the literature |l7oj 1182 ) : 


f m = ~V* = 


1 sc. 

2 Sy 


-o 


I A-10) 


Concurrent to this, it Is necessary to model the effect of the spring used to return 
the actuator to its initial position when the voltage difference across the couih- 
drive lingers is removed. One design for the spring suspension system Is the folded 
flexure design shown in Fig. A.7 |177) . 
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Figure A.7: Folded flexure design commonly used for large displacement electro¬ 
static comb-drive actuators |I77| . 


The force that must be applied to the spring to achieve the desired deflection 
can be computed by modeling the individual springs as beam elements with fixed 
ends. Then the force of the spring Is modeled as presented in the literature 117l)| : 
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F ipr ,„ g = 6^=26^. (All) 

By setting the actuator force equal to the spring force, it is possible to design the 
electrostatic actuator such that the desired displacement can be achieved. The 
parameters that are able to be designed are the number of fingers, ». the spring 
length. L the spring width. Is- the structural thickness, Iia- the finger gap. 2 o, and 
the operating voltage. V. 

The scaling effects can be particularly significant and thus should be considered 
during the design process. For example, the force necessary to drive a fixed-fixed 
beam spring as proposed in the present project is greatly subjective to the effects 
of geometric scaling. As the overall size of the beam suspension system is scaled 
down by an order of magnitude from the original size, the force necessary to 
achieve the new scaled displacement reduced by two orders of magnitude. 

This effect is beneficial since it allows the electrostatic actuator requirement to 
be greatly reduced for even small changes in the suspension system size. One 
possible realization of this is in terms of the driving voltage for the electrostatic 
actuator. Assuming consistent geometry of the actuator, a smaller requirement in 
force would mean that the driving voltage could be reduced, thus resulting in a 
more efficient device. Additionally, it is interesting to note that the force output 
from the actuator is independent of geometry scaling assuming that the finger gap 
scales proportionally to structural thickness of the actuator fingers. Therefore, the 
limitation of geometric scaling of the actuator is only limited by the constraints 
of the fabrication process used to produce the device. 

Based upon these principles, a valve design was developed that would be operated 
by an electrostatic comb-drive actuator. The design included 100 pm long fingers 
in order to achieve a desired translation of at least GO /nn. Additionally, the 
orifice of the valve was sized at 50 /nn by 200 pm in order to achieve a flow rate 
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of 9 mL/min. The layout shown in Fig. A.8 represents the three layer masks 
necessary to manufacture this valve design. The green background represents the 
silicon substrate upon which the valve is const meted. The yellow layer locates the 
positions of the valve orifices and tints represents the area which will be etched 
away from the silicon wafer. The blue layer shows the st ructure of the valve itself, 
which is produced from polysilicon. Finally, the red layer represents the anchor 
which allows the polvsilicon to adhere to the wafer below. The details of this 
process are described in the following sect ion. 



Figure A.8: Mask layout for valve design. 


Processing 


A possible fabrication process is presented in the following flowchart. To ac¬ 
company each of the steps, a cross sectional model of the wafer throughout the 
fabrication process is presented alongside each description. Information regarding 
specific deposit ion/etch rates and recipes to be used have been adapted from the 
website of the Nanofabrication Center at the University of Minnesota - Twin < 'itios 
(183| . For convenience, a key correlating the different cross-hatching patterns used 
to the appropriate material is shown below in Fig. A.9 
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L Silicon Wafer | Silicon Niliide 
| Silicon Dioxide , Polysiicon 
| Gold 

Figure A.9: Key to he used to identify materials in the process design figures. 

1. Start with a bare (100) silicon wider. 

(a) Wafer should he cleaned using the RCA clean procedure. 

Si 

Figure A. 10: Bare silicon wafer. 

2. Deposit 2000 A of silicon nitride by LFC'VD. 

(a) Deposition for (i() min. using the default Low Stress Nitride recipe. 

Si 

Figure A.11: Wafer after silicon nitride deposition. 

3. Deposit 1 pm of silicon dioxide by LPCVD. 

(a) Deposition for 125 min. using the PSG recipe. 

(b) Anneal the wafer. 
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Si 

Figure A. 12: Wafer after silicon dioxide deposition. 

4. Perform photolithography on silicon dioxide. 

(a) Deposit, photoresist (PR-1813) by spinning. 

(b) Expose photoresist to UV using mask (dark field). 

(c) Etch photoresist using photoresist developer. 

(d) Etch unprotected Silicon Dioxide using HE etch. 

.. 

Si 

Figure A. 13: Wafer after silicon dioxide etching. 

5. Deposit polysilicon by LPC'VD. 

(a) Deposition for 7 hours using a combined doped/undoped Polysilicon 
recipe. 

(b) Anneal the wafer to reduce the chance of undesired deformation. 

—F—7—J? 7— 7 L, V 7— 

1/ / / / / / / / / / / 4 

i J.TTm.VL f. Ml. ■ ■■■>■ ... 

Si 


Figure A. 14: Wafer after polysilicon deposition. 
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(>. Perforin photolithography on polysilicon. 

(a) Deposit photoresist (PR 1813) by spinning. 

(b) Expose photoresist to UV using mask (light field). 

(c) Etch photoresist using photoresist developer. 

(d) Etch unprotected Polysilicon using R1E. 



Si 


Figure A.15: Wafer after polysilicon etching. 


7. Deposition of a reflective gold layer by electron-beam evaporation. 



Si 


Figure A.1G: Wafer aft or gold deposition. 

8. Remove sacrificial silicon dioxide by HE etch. 

(a) Wafer will need to go into the Critical Point Dryer to remove any 
remaining etchants. 
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Figure A. 17: Wafer after removal of sacrificial layer. 


9. 


Wafer diced into individual units and assembled with packaging to form final 
product. 


By following this procedure through step six. a prototype of this design was pro¬ 
duced and is shown in Fig. |A.li j 



Figure A.18: Valve design prototype shown at ox magnification. 














Valve Design Considerations 


224 


Experimentation with MEMS processing using the above procedure revealed that 
there arc some limitations of the methodology that need to be given greater con¬ 
sideration. Primarily, the use of a 1 pin sacrificial silicon dioxide layer was found to 
be insufficient and led to high sfictiou susceptibility. This issue may he remedied 
by increasing the sacrificial layer thickness to 2 pm. Further, by analyzing similar 
processing procedures described in the literature, it may be necessary to increase 
the thickness of the structural polysilicon layer iTTTTj : however, the method by 
which this can be achieved has yet to he determined. Other avenues of fabrica¬ 
tion have berm presented in the literature and shown to produce high-aspect-ratio 
comb-drive structures |18-l| . 


Appendix B 
IFM Code 


Two-Dimensional IFM Code 

The following code is written in Matlab M-code. This code accepts input beam 
end conditions and computes the most likely beam shape and end loads for the 
two-dimensional planar case. 

i function moment _L,forceZ,forceX,x] — energyMethod_standalone( 
thetaL.zL.xL.xO) 

.• % This junction returns the potential energy in the beam, the moment 
produced at the tip, and the force produced at the tip for some 
given input tip location and orientation . The solution is found by 
minimising the potential energy in the beam to find the optimal ( 
most likely ) beam s/ia/ie. 

* 

• % Input beam material and geometric properties. 

>, E = 28e9; 
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c; = 10c9; 

r I = 0.00045~4*pi/64; 
. A = 0.00045‘2*pi/4; 
u L = 0.05; 


>j % Convert thctoL to radians. 

k tliotaL = thctaL*pi/180; 


IS* 

«t. %%%%%%%% Optimisation 

IT 

i. Aeq = [L‘(i.L'5.L I.L 3.L'2.L|; %%% Theta end constraint specified. 
iu beq = tlietaL; 

XI 

.i % Optimization options set. and submitted to the fmincon function. 
options = opt unset (’ Algorithm'.* interior —point'.'Display', * oil ’, ‘ 
TolCon',0.001,' MaxFunEvnls’, 100000); 

ji [x, Penergy, exitflng .output .lambda — finincon(Bc*uergyFl ’N.x0.[j.[.Aeq. 
hni . 0. [1, *! (x )energyC’( )N(x,zL.xL),opt ions); 

21 

% \'irtual work principle to calculate end load. 
ju moment _L = lambda.eqlin: 

-•t foreeX = lambda.eqnonlin(l): 

-•* forceY = lambda.cqnoulin(2); 


%%%%%%%% energyFUN function 

12 

» function PE = energyF UN (x) 
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i4 % This junction returns the potential energy of the beam for a provided 
set of input parameters. 

it. % Input beam material and geometric properties. 
it E = 28e0; 

4. c; = 10ef); 

,, I = 0.00045‘4* pi/64: 

... A = 0.00045'2*pi/4; 

.. L = 0.05; 


u %%%%%%%% Parameterization of Theta 
IS. s = 0:L/100:L; 

l(i 

•r X Derivative of theta with respect to s. 

i* dtheta = x(l:6)*[6*s/5;5*s.’4:4*s/3;3*s.*2;2*s;onea(l,101)]; 


i % Potential energy function. 

-• integrand = E* I •'dtheta. '2; 

■' PE — 0.5*trapz(s.integrand); 


•t. %%%%%%%% energyCON function 

SI 

function c.ceq] - euergyC’ON(x,zL,xL) 


e .1 c 
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. L = 0.05; 

64 

s = 0:L/100:L; 

M 

,.t theta= x(l).*s.'6+x(2).*s.'5+x(3).*s.~4+x(4).*s.\3+x(5).#s.‘2+x{6).*s; 

C*JI 

-■ intx = cos(theta); 

7>i intv = sin(thcta); 

7j % Computation of end conditions given input beam shape. 

7 « xs = trapz(s.intx); 

74 ys = trnpzfs.iutv); 

ri. % Apply the constraints to the result . 

7 T ceq(l) = zs-zL; 

7 * ceq(2) = xs—xL: 


Three-Dimensional IFM Code 

The following code is written in Matlab M-code. This code accepts input beam 
end conditions and computes the most likely beam shape and end loads for the 
three-dimensional planar case. 

i function moment.phi.momentjheta,moment_psi,forceX.forceY.forceZ] 

= energy Method_standnlone3D(phiL,thetaL.psiL,xL.yL,zL) 
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.• % This Junction ieturns the potential energy in the beam, the moment 
produced at the tip , and the force produced at flic tip for some 
given input tip location and orientation. The solution is found by 
minimizing the potential energy in the beam to find the optimal ( 
most likely } beam shape . 

A 

> % Input beam material and geometric properties. 

E = 28e9; 

„ CJ = 10c9; 

r 1 = 0.00045'4-* pi/64; 

* A = 0.00045‘2*pi/4; 

u L = 0.05; 

aii 

x0 = [0.0.0.0.0.0,0.0.0.0.0.0.0.0.0.0,0,01; 

aa 

ii % Convert angles to radians . 

i< thetaL = thctaL*pi/180; 

i'. pliiL = phiL*pi/180; 

H. psiL = psiL'pi/180; 


IK 

ik %%%%%%%% Optimization 

an 

.i A<*j = [L 6.L 5.L 4.L 3,L 2,L .0.0,0,0.0.0.0.0,0.0.0,0; 

0,0.0,0,0,0, L‘(i.L‘5,L'4.L'3.L'2,L,0,0,0,0,0,0; 

11.0.0,0,0,0,0,0,0.0.0.0. L‘6.L‘5,L‘I.L‘3,L‘2.L]; %%% pin. 
theta, and psi end constraint specified . 
beq = jpliiL;thetaL;psiL]; 

... % Optimization options set. and submitted to the fmincon function. 
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.•r options = optimset(’AlgoritlunV interior — point VTolCon',0.001,’ 
MaxFuiiEvals',100000); 

“ [* .Penergy.oxitflng,output .lambda — finincon(tfieuergyFUN.xO. (].[ .Aeq. 
b «l'D*lI’'* 1 '( x ) ener ©’CON(x,xL.yL.zL 1 the t iiL,|iliiL,psiL).qptions); 
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in % l 'irtual work principle to calculate end load. 

ii mouient_phi = lambda.eqlm{ 1): 

u moment .theta = lnmbda.eqlin( 2 ); 

ii moment_psi = 1 ambda.eqlin(3>; 

i« forceX = lambda.eqnonlin(l): 

is forccV = lambda.eqnonlin(2): 

v. forceZ = lainbda.equonIiu<3); 

»T 

ill 

in %%%%%%%% energy FUN function 
II 

i3 function PE = energyFUN(x) 

ti % Tin* function returns the potential energy of the beam for a provided 
set of input parameters. 

44 

is % Input beam material and geometric properties. 

... E = 28e9; 

,r (! = 10e9; 

.* I = 0.00045 *4*pi/64: 

. 0 A = 0.00045‘2*pi/4; 

. L = 0.05: 


■i 

S3 


%%%%%%%% Parameterization of Phi. Theta, and Psi 



s = 0:L/100:L; 


% Parameterized variables. 

theta = x(7).*s. *6+x{8).»-s. * 5+x(9).*s.' 4+x( 10).*s. ‘ 3+x( 11 ).*s.' 2+x 

(12).*s; 

dphi = x(l:6)*|6*s/5;5*s/4;4*s.'3:3*s.'2;2*s;oues(l,101)]; 
dtlieta = x(7:12)#[()*s/5;5*s/4|4*s."3;3*B.'2;2*s;ones(l,101)]; 
dpsi = x(13:18)*’(>*s. "5;5“s. '4;4*s/3;3*s.’2;2*s;ones( 1,101) ]: 

% Pottnliul energy June lion. 

integrand = E*I*(dtheta/2+dphi/2.*cos(theta)/2)+2*G*l*(sin<theta) 
.*dphi + dpsi)/2; 

PE = 0.5*trapz(s. integrand); 


%%%%%%%% encrgyCON function 
function jc.eeq] — energy('ON(x,xL,yL.zL) 

' = 0 ; 

L = 0.05; 
s = 0:L/100:L; 

phi = x(l).*s. 'G+x(2).*s.'5+x(3).*s."4+x(4).»s/3+x(5).*s.'2+x((i).*s: 
t beta = x(7).«. *6+x(8).*s.' 5+x(9).*s.' 4+x( 10). *s.' 3+x( 11 ).*s.' 2+x 

(12).*s: 

psi = x(13).«/C+x( 14 ).«s/ 5 +x( 15 ). 4 s.' 4 +x< 16).*s.~3+x(17).*s/2+x 

(18).*s: 
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i intx = sin(tlieta); 

.• inly = —sin(phi).*cos(tliet.a); 
i intz — coa( theta). *co3(phi); 

. % Computation of end conditions given input beam shape. 
. x-s = trapz(s.mtx); 
r vs = trapz(s.inty); 

~ zs = trapz(s.intz); 

wii % Apply the constraints to the result. 

■ji ceq(l) = xs—xL: 

ceq(2) = ys-.vL: 

•« ceq(3) = zs—zL; 




Appendix C 


Abaqus Scripting Code 


Automation of Abaqus Simulations 


The following code is written in the Python programming language. This code 
consists of a series of nested FOR loops which iterate through the various model 
parameters, compute the shape of the bladder segment, mesh the segment, apply 
the load, and submit the simulation. The results of each run of the simulation are 
written to a unique .ODD lile. 

i # — *— coding: mbes — * — 

.• from part import * 

i from material import « 

* From section import » 

'. from assembly import * 

ii from step import » 

t from interaction import « 

• From load import * 
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from mesh import 

from optimization import * 

from job import * 

from sketch import « 

from visualization import * 

from connect or Behavior import 4 

import math 

^import time 

# SET VARIABLES 

EdgeLengths = [0.25.0.5,0.75,1.0.1.25.1.5] 

BraidAugles = [5.15,25,35,45,55] 

Thicknesses = [0.05,0.1,0.15,0.2,0.25] 

Pressures = [0.2,0.4.0.6.0.8.10] # Units art MPa 

E = G.5 
nu = 0.5 

i=l 

for EdgeLength in EdgeLcngtlis: 
for BraidAngle in BraidAugles: 
for Thickness in Thicknesses: 
for Pressure in Pressures: 

x — EdgeLeugth*matk.cos(BrnidAngle*mnth.pi/ ISO) 
y = EdgeLength*math.sin(BraidAngle«math.pi/180) 
xmldpoint — EdgeLength/2*inath.cos(BraidAngle*uiath. 
pi/180) 

ymid point — EdgeLeugth/2*math.sin(BraidAngle«>math. 
pi/180) 
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# CREATE THE MODEL AND DEFINE EDGE 
LENGTH. HR AID ANGLE. AND THICKNESS 

mdb.modelsj'Model—rj.ConstrainedSketcli(nnine=' 
__profile_\ shectSize=5.0) 

mdb.modelsj'Model— l'j.sketches|' ..profile.. ’ ]. Line( 
point 1=(—x, 0), point2={ 

0, .V)) 

mdb.modelsj'Model—1’ .sketches!' ..profile.. ’ ]. Line( 
point 1=(0, y). point2=( 

x, 0)) 

mdb.niodelsJ’Model— 1' .sketches!' ..profile.. ’ ]. Line( 
point l=(x, 11 ), point 2 =( 

", -y)) 

uidb.modelsJ'Model—l'j.sketches[’ ..profile.. ' ]. Line( 
pointl=(0, — y). poiut2=( 

-x. 0)) 


mdb.madels['Model-rj.Part(dhnensionality=THREE_D 
. n!une='Pnrt—1’. type= 

DEFORMABLE_B()DY) 


mdb. models 'Model— 1’ .parts['Part—1' . 

BaseSolidExt rude(depth=Thickness, sketcli= 
inclb.niodels[’Model—r].sketchesf' __profile_ ')) 
del mdb.modekf’Model-rl sketches!' ..profile.. ’] 


# SPECIFY THE ELASTIC MATERIAL 
PROPERTIES 

mdb.modelsj'Model—1’ .Material(name='Material— 1") 
mdb.modelsj’Model—rj.raateriaLs[’Materinl-r).ElnKtic( 
t!ible=((E, im), )) 



# SECTION AND ASSIGN TO PAST 
mdb.modelsi'Model—r].HomogeneousSoIidSection( 

material=’Mnterial—1\ name= 

’ Sect ion — 1'. t hickness=None) 
mdb.modelsi’Model— r].parts['Part— 1’. 
SeetionAssigiiment (offset =0.0. 

oiiset Field =", offset Type=MIDDLE_SURFACE, 

region=Region( 

cells =mdbjnodels[’Model—l’].parts 'Part—1’].cells. 
lindAt. (((().(). 

0 .0, 0.0), ). )). sectionName='Section-l’. 

t hickness Assignment= 

FROMJSECTION) 

# CREATE INSTANCE AND STATIC STEP 
indb.models Model— 1 ].rootAssembly. 

DatumCsysBy Default (CARTESIAN) 
mdb.models "Model— 1’ .rootAssembly.Instauce( 
dependent=OFF. mune=’Part —1 — 1’, 
part =indb.models[‘\Iodel— 1 ’ | .parts] 'Part — 1 )) 
mdb.models 'Model — 1'j.StaticStep(name= ‘St ep — 1'. 
previous=’Initial’) 

# SET THE PRESSURE LOAD AND BOUNDARY 

CONDITIONS 

mdb.models [‘Model-1"). Pressure! amplitude=UNSET. 
creat eSt epN ame=’St ep—1", 

distributionTy|je=UNIFORM, field=’\ magnitude^ 
Pressure, naine=‘Load—1‘, region= 

Region ( 



sidelFiiccs=mdb.inodels|'Model— l 1 ].root. Assembly. 

instancesf'Part -1 — r].faces.findAt( 

(( 0 - 0 , 0 . 0 , 0 . 0 ), ), ))) 

mdb. models} Model—T .PinncdBO(createStepName=' 
Step—1', locaICsys=None, name= 

BC’— 1’, region=Region ( 
faces=mdb.models[Model—l’j.root Assembly. 

instances!' P.art -1 — 1’].faces, find At ( 

((— xmidpoint. yinidpoint. Thickness/2). ). (( 
xmidpoint, yiuidpoint, Thickness/2), ), (( 
xmidpoint, —yinidpoint. Thickness/2), ). (( — 
xmidpoint. —yinidpoint. Thickness/2), ), 

))) 

# SEED THE INSTANCE AND MESH IT (MIGHT 

CHANGE MESH SIZE HERE) 
mdb.models]Model— T .raotA.sseinblv.seedPart Instance( 
deviat ionFact or=0.1, 

regions—(nidb. models[' Model — T j .root Assembly. 

instances!'Part—1—T], )• 

0.05) 

mdb.inodelsjModel— T .rootAssembly.generateMesh( 
regions=( 

indb.models[*Model—T].root Assembly.instances!' 
Part-1-1’], )) 

# CREATE THE JOB AND SUBMIT 
jobunme = Mob—%s’ %\ 

mdb..Job(ntTime=Noue. contactPrint=()FF, description 
=", echoPriut=OFF, 
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explicit Precision =SINGLE, 

get\leinoryFrouiAnaIysis=True, historyPrint— 
OFF, 


v. inemory=90, nieinoryUnits=PERCENTAGE, model 

=’Model—1\ modelPrint=OFF, 

•.» inultiprocessingMode=DEFAULT. name—jobimmc, 

nodalOutputPrecision=SINGLE, 

or immCpns=i, queue=None, scratch='\ tvpe= 

ANALYSIS, userSubroutine=”, 

« waitHouis=0. waitMinutes=0) 

oo nuib.jobs[jobname].suhuiit(consistencyCUccking=OFF) 

■ ni #time.sleep(30) 

iot i=i+l 

.... # END OF FOR LOOPS 


Post Processor Script 


After all of the simulations are complete, the following post processing script 
written in the Python programming language is used. The script opens each of 
the .ODB files generated previously, extracts the von Mises stress data, anti writes 
it to a text file for later use. 

. # READ OUTPUT FROM ODD RESULTS FILES 
j import odbAccess 

i maxStress=[] jf Initialize the variable of interest. 

for j in range(l): # Just need to edit, this with actual number of Job 
output files. 
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« jindex=j+l ff Adjust index to mutch job name convention used in 
jmiametric study. 

r jobuamc=’Job—9fs.odb’ %jindex 

* stressField =sessiou.openOdb(naine=jobunine).steps['Step—1’]. 

frmncs[-I].field()utputB[’S’] # Define path to stored data. 

•j stress Vector=[] # Initialize a temporary variable. 

for i in range(len(stressField .values)): if Loop through all the 
elements in the. model. 

von\Iises=stressField.vnlues[i ].inises ff Extract variable of 
interest (von .Wises stress in this case). 
u stressVector.insert(i.vonMises) ff Put that variable in a the 

temporary list. 

n maxS=innx(stress Vector) if Looking for max value in this case. 

i» maxStress. insert(jiiulex,mnxS) ff Store that max value in a list. 

,, ff end for loop 

St 

ir output =lile( '/sonie/path/Results.txt\'w') # Once all of the job files 
are processed, write max stress list to text file . 

For item in maxStress: 
in print»outpiit, item 
» if end output for loop 
-m output.eloseO if Close the text file . 



